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Chapter 1 



Introduction 



Sixty years after its invention by the Japanese mathematician Kiyosi Ito, stochastic ana- 
lysis has found a wide range of apphcations varying from pure mathematics to physics, 
enginering, biology and economics. Stochastic differential equations are an important tool 
for solving problems in pricing options at the stock market, finding solutions to boundary 
value problems, filtering signals in electrical enginering and for solving many other prob- 
lems. In 1984 Hudson and Parthasarathy published a paper [54] extending the definition 
of ItcVs stochastic integral and its subsequent calculus to the non-commutative world of 
quantum mechanics. The goal of this thesis is to apply this generalized stochastic calculus 
to problems in quantum optics. 



1.1 Some spectral theory and quantum probability 

The question of the mathematical foundations of quantum mechanics was settled by von 
Neumann in a series of papers written between 1927 and 1932, culminating in his book 
" Mathematische Grundlagen der Quantenmechanik" [87]. In this series of articles von 
Neumann developed spectral theory for normal operators on a Hilbcrt space, i.e. operators 
that commute with their adjoint. Avoiding tedious considerations regarding domains of 
operators we will now first focus on spectral theory for bounded normal operators on a 
separable Hilbert space H. 

Denote by B{M) the algebra of all bounded operators on H and let 5 be a subset of 6(EI). 
We call the set S' := {R G ^(H); SR = RS e S} the commutant of S in B(a). A 
von Neumann algebra .4 on H is a *-subalgebra of B{M.) such that A equals its double 
commutant, i.e. A" = A. It follows from von Neumann's double commutant theorem, cf. 
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[55], that a von Neumann algebra is a C*-subalgcbra, i.e. a *-subalgebra closed in the 
operator norm topology, of B{M.) that is closed even in the weak operator topology. 

It immediately follows from A" = A that the identity 1 G B{Il) is an element of the von 
Neumann algebra A. A state on ^ is a linear map p : A ^ C such that p is positive in 
the sense that p{A*A) > for all ^4 € .4 and such that p is normalised p{l) = 1. A state 
is called normal if it is weak operator continuous on the unit ball of A. The following 
theorem, see [56] sections 9.3, 9.4 and 9.5 for a proof, is at the heart of spectral theory. 

Theorem 1.1.1: Let C be a commutative von Neumann algebra and p a normal state on 
C. Then there is a probability space (Q, S, P) such that C is * -isomorphic to L°°{^1, S, P), 
the space of all bounded measurable functions on U. Furthermore, if we denote the *- 
isomorphism between C and L°°(r2, E, P) by C ^ C, we have 

p{C)= [ C^F{dio), CgC, 
Jn 

where C^ denotes the function C, evaluated at u, i.e. C^ = C,(a;). 

Since a normal operator, i.e. an operator that commutes with its adjoint, generates a 
commutative von Neumann algebra the above theorem can be applied. The simplest 
example of a normal operator is a Hermitian operator ^4 on a finite dimensional Hilbert 
space. The above theorem then states that A is equivalent with a function A, and that 
the algebra generated by A is isomorphic to the algebra generated by this function A,. 
Indeed, after diagonalisation, A is equivalent with the function that maps i, standing for 
diagonal entry number i, to its corresponding eigenvalue. Note that the above theorem 
also covers simultaneous diagonalisation of several commuting Hermitian operators. 

Given a probability space ($7, S,P), we can study the commutative von Neumann algebra 
B := L°°(f],S,P), acting on the Hilbert space L^(J7, S,P) by pointwise multiplication, 
equipped with the normal state p given by expectation with respect to the measure P. 
The pair (B, p) faithfully encodes the probability space (O, S, P) [70]. Indeed, the u-algebra 
S can be reconstructed (up to equivalence of sets with symmetric difference 0, a point on 
which we will not dwell here) as the set of projections in B, i.e. the set of characteristic 
functions of sets in S, and the probability measure is given by acting with the state p on 
this set of projections. We conclude that studying commutative von Neumann algebras 
equipped with normal states is the same as studying probability spaces. This motivates 
the definition of a non- commutative or quantum probability space as a von Neumann 
algebra equipped with a normal state. 

Prom the point of view of physics it is violation of Bell's inequalities [21] in for instance the 
Aspect experiment [6] , that motivates the study of non-commutative probability. See for 
instance [64] for a clear exposition of the point that this violation can not be accounted for 



1.1. SOME SPECTRAL THEORY AND QUANTUM PROBABILITY 



3 



by a local classical probabilistic model. The way out if one wants to preserve locality is to 
weaken Kolmogorov's axioms of probability ending up with a non-commutative probability 
theory, namely quantum mechanics. 

Let {B, p) be a quantum probability space. The events are given by the projections in i3, 
i.e. elements satisfying E"^ = E = E*. Two events E and F are called compatible if EF 
is again an event, which is equivalent to EF = FE. The interpretation is such that only 
events that are compatible can occur simultaneously. Note that in classical probability, 
i.e. B is commutative, all events are compatible. If E and F are compatible events then 
EF stands for the occurence of both E and F and E W F := E + F — EF stands for the 
occurence of E or F or both. When an event E has occurred we have to update the state 
p to pe with Pe{A) = p{EAE)/p{E), the state conditioned on E. 

In classical probability theory a random variable X is a measurable map from a probability 
space (Q, S,P) to some other measure space (Q', S'). The probability distribution of X is 
given by 

Px : S' ^ [0, 1]:S^ F{X-'^{S)) . 
Algebraically X is completely determined by the pull back 

jx: L°^(fl',S',Px) ^L°°(J7,E,P) : f^foX. 

This motivates the definition of a random variable on a quantum probability space {B, p) 
as a *-homomorphism j from some other von Neumann algebra A to B mapping 1_4 to 
Is- The probability distribution of j is then given by o" := p o j on A. A quantum 

stochastic process [I] is a family {jjjtgT of random variables indexed by time T which is a 
linearly ordered set such as N, Z, M or R_|_. See [61] for theory on non-commutative Markov 
processes. 

Let j : A ^ B he a random variable on {B,p). If ^ = L°°(M) then j is called a real 
valued random variable. The following brief exposition of spectral theory for unbounded 
selfadjoint operators shows that there is a one-to-one correspondence between real valued 
random variables and selfadjoint operators. Given j we can define a projection valued 
measure on the Borel sigma-algebra S of R by 

EiS) := jixs), -5GS, 

where xs stands for the characteristic function of S, i.e. the function that is 1 on S* and 
elsewhere. We can now define a selfadjoint operator A by 

A = [ XE{dX). (1.1) 

Conversely, since the spectrum of a selfadjoint operator A is real, we can define bounded 
operators r+ := (^4 -|- il)~^ and T_ := {A — The operators T_|_ and T_ generate a 

commutative von Neumann algebra A- Prom Theorem 1.1.1 it follows that there exist a 
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measurable space {Q,Ti') and a * -isomorphism i : L°°(i7,S') — >■ A- Define a measurable 
function A, from O to the extended reals M by 

Then wc can define a projection valued measure on the Borel sigma algebra of M, generated 
by the sets (—00, A] with A G M, by 

£^((-00, A]) := i{x[A,<x]), 

and in this way equation (1.1) reappears. Furthermore the spectral measure, in its turn, 
uniquely determines a real valued random variable j : L°°{M) B. The above ex- 
positions shows three equivalent ways of characterising real valued random variables or 
observables, i.e. by a selfadjoint operator A, by a projection valued measure E, and by a 
*-homomorphism j. A fourth way of looking at random variables is provided by Stone's 
theorem, see for instance [55]. 

Theorem 1.1.2: (Stone's theorem) There is a one-to-one correspondence between 
strongly continuous unitary representations t Ut of the abehan group M into a von 
Neumann algebra A and selfadjoint operators A affiliated to A, i.e. having all its spectral 
projections in A, such that 

Ut = exp{itA) := [ e'^^E{dX). 

JR 



Let A and B be von Neumann algebras. We denote their Banach space duals by A* and B*, 
respectively. A linear map T : A ^ B determines a dual map T* : B* ^ A* hy p ^ poT. 
An element A of a C*-algebra is called positive if it can be written as A = B*B for some B 
in that C*-algebra. Recall that von Neumann algebras are a special class of C*-algebras. 
The map T* maps states on B to states on A if and only if T is positive, in the sense that 
it maps positive elements of ^ into positive elements of ^S, and T{lj[) = Ig. The map T is 
said to be completely positive if for all n G N the map T Id„ : A M„(C) ^ B(Si M„(C) 
is positive. An operation T is a completely positive map such that T'(l_4) = Ig and its 
dual T* is also called an operation. An operation T* : B* ^ A* represent some physical 
procedure which takes as an input a state on a quantum system described by B and turns 
out a state on the quantum system A. The complete positivity ensures that operations 
on a quantum system A can always be extended to operations on A<Si W, the quantum 
system A with its environment W. 

The operational approach to quantum mechanics was pioneered by Davies and Lewis [34] , 
[33], [37]. In this approach all physical procedures performed on quantum systems are 
given by operations between their algebras of observables. An example of an operation is 
a random variable, i.e. a *-homomorphism i from an algebra A to B. Suppose B = A'S>W 
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and i = ld<E> lyVj then the dual of i is what is known in physics as a partial trace, i.e. 
restriction of the state to the smaller system A<Si lyv- A second example is adjoining with 
a unitary operator or a projection, i.e. time evolution or conditioning on a measurement 
result, respectively. A third example of an operation is a conditional expectation, the topic 
of the next paragraph. 

Let B be a von Neumann subalgebra of a von Neumann algebra A. A conditional expecta- 
tion of A onto H is a linear surjective map 6 : A ^ B such that = 6 and \\6\\ = 1. We 
are interested in a conditional expectation £p : A ^ B that leaves a certain state p on ^ 
invariant, i.e. p o £p = p. However, such a conditional expectation does not always exist, 
but if it exists it is unique [84] . The interpretation is as follows, if we start with a quantum 
system ^ in a state p and we gain (for example by observation) the additional information 
that S is in a state cr, then a o £p is the updated state of A [62]. It follows from [85] that 
£p is an operation and that it satisfies the module property: £p{BiAB2) = Bi£p{A)B2 for 
all Bi,B2 G B and A £ A. An example of a conditonal expectation is given by taking 
A = B ®W and £ = ldL®<j) with (j) a state on W. If .4 is a commutative algebra then we 
retrieve the conditional expectation of classical probability theory. 



1.2 Stochastic calculus on Fock spaces 



Let EI be a Hilbert space and think of its elements as the pure states of one particle in 
quantum mechanics. The particles we are interested in are photons. Since photons are 
bosons they have to be described by symmetrised wave functions. This motivates the 
definition of the symmetric Fock space over H as 

oo 
n=l 

It describes situations where the number of particles present is arbitrary. For every / € H 
we define the exponential vector e(/) G J^{M) by 




The linear span V of all exponential vectors is a dense subspace of jr(EI). On the dense 
domain T> we define for all / G M an operator W{f) by 

W{f)e{g) := exp ( - {f,g) - \\ft)e{f + g), ge H, (1.3) 

which is isometric and therefore uniquely extends to a unitary operator on J^{M). The 
operators W{f) : T{M) — J^{M.) are called Weyl operators and they satisfy the following 
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Weyl relations 

1. WifT = W{-f), / G M, 

2. W{f)W{g) = e^p{-ilm{f,g))W{f + g), f,9eM. 



For all / G H the family {W{tf)}tem forms a one-parameter group, continuous in the 
strong operator topology. Therefore it follows from Stone's theorem 1.1.2 that for all 
/ G H there exists a selfadjoint B{f ) such that 

W{tf) = exip{itB{f)). 

The operators B{f) are called field operators, see also Section 1.3. The domain of the 
operator B{fk)...B{fi) contains V for every fi,...fk G H and k E N (cf. [75]). For 
/, 5 G H and t G M it follows from the Weyl relations that on the domain V 

1. B{tf)=tB{f), 

2. B{f + g) = B{f) + B{g), (1.5) 

3. [B{f),B{g)]=2ilm{f,g). 



Let us take H = L^(M) and let us consider it in the canonical way as a real space, then 
the pair (H, Im(-, •)) forms a symplectic space. Denote by Hq the real subspace of EI given 
by {/ = (/i,/2) E H; /2 = 0}. From (1.5.3) it immediately follows that the family of 
operators {B{f); f G .ffo} is commutative. 

As in the previous section we can consider the bounded operators T'(/)+ := {B{f) + il)~^ 
and T{f)- := {B{f) — for / G Hq. They generate a commutative von Neumann 

algebra C. We denote by cj) the state on C given by 

</,(C) := ($,C$), with $ = leOeO... G J^(M). 

Using Theorem 1.1.1 we see that (C, cf)) is isomorphic to L°°(p., S,P) for some probability 
space (J7,S,P). We denote this isomorphism from C to L°°(i7,S,P) by C ^ C,. As in 
the previous section, see equation (1.2), we can let the operator B{f) (/ G Hq), affiliated 
with C, correspond to a measurable function B{f), on Q taking values in the extended 
reals R. 

The joint characteristic function of the random variables B{fi),, B{f2),, ... , B{fk), is for 
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xi, . . . ,Xfc G M given by 

]g|-gixiB(/i).gia;2B(/2). 

k 

1=1 

where we used Theorem 1.1.1 in the second step and equation (1.5.1) and (1.5.2) in the 
fourth step. For i G M we define a random variable Bt := B{x[o,t))» where X[o,t) stands 
for the function that is 1 on the interval [0, t) and elsewhere. For si < ti < S2 < t2 
it immediately follows from (1.6) that the joint characteristic function of the increments 
Bt, - Bs, = 5(x[5i,ti)). and Bt^ - Bs^ = B{x[s2,t2))» factorizes 

^^ix{B,,-B.,)^iy{B,,-B.,)^ = exp ( - ^{x\tl - Si) + y\t2 " S2))) = 

exp ( - exp ( - = E[e-(^*i-^n)]E[e'^(^*^-^-)] , 

i.e. {-Bt}t>o is a process with independent increments. Furthermore, from its characteristic 
function it follows that the increments Bf — Bs are normally distributed with mean zero 
and variance t — s. Summarizing, the process {Bt}t>o is a Brownian motion. 

The idea to simultaneously diagonalise the fields in the family {B{x[o,t))j i > 0} is implicit 
is some of the earliest work in quantum field theory. However, Segal [78] in the 1950s was 
the first to emphasise the connection with probability theory. See also [80] for a nice 
review of the importance of these ideas for Euclidean quantum field theory. 

Fix a in the interval [0,7r). Then, as in the previous paragraph, the family of operators 
5(e*"x[o,t)) is commutative. This means there is a probability space {Uaj'^a^^a) on 
which the operators -B(e*"x[o,t)) are realised as random variables i?(e*"x[o,t)), =• -^f • 
a similar fashion as in the previous paragraph it can be shown that the process }t>o 
is a Brownian motion. However, for two different alphas, i.e. ai and 0:2, it follows from 
equation (1.5.3) that the Brownian motions B^^ and B^^ do not commute. For instance 
we have 

[BlB^^^]=2imm{s,t}. 

Expression (1.5.3) is known as the canonical commutation relation. It is clear that the 
Weyl relation (1.4.2) is an exponentiated form of the canonical commutation relation. Weyl 
operators have the advantage that they are bounded, so that there are no domain problems 
involved. That is why we abstractly define the C* -algebra of canonical commutation 



k 

.e^^fcB(/fc)$^ = ^^>,exp(i^XiB(/i))$^ 

i=l 

1 ^ 



(1.6) 
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relations CCR{H, a) with respect to some symplectic space {H,a) as the C*-algebra 
generated by abstract elements W{f) (/ G H) satisfying the relations in (1.4) where H 
and Im(-, •) should be replaced by H and a, respectively. It follows from [81] that this 
C*-algebra exists and moreover that it is unique up to isomorphism. 

Equation (1.3) defines a representation of CCR{M, Im(-, •)) on the Fock space J^{M). This 
representation is the GNS-representation (cf. [75]) with respect to the vacuum state, i.e. 
the state that maps W{f) to exp(— l/2||/p). In Chapter 4 we will encounter GNS- 
representations of CCR(M,Im{-, •)) with respect to other states than the vacuum. Given 
a representation of CCR{Il,lm.{-, •)), we can define the von Neumann algebra W gener- 
ated by the W(f) (/ G H) represented as operators on the representation Hilbert space. 
For the vacuum GNS-representation it turns out that W is the whole algebra of bounded 
operators on J^(M.) cf. [75]. The algebra W = B[J-(M.)^ equipped with the vacuum state is 
the quantum probability space with which we will describe an electromagnetic field in the 
vacuum state, see also Section 1.3. More information about CCi?-algebras can be found 
in [25], [26] and [75]. 

Apart from the Brownian motions the algebra W = 6(.7-"(L^(]R))) contains Poisson 
processes [41], a point that we will investigate below. Let Pt : L^(M) L^(M) denote the 
projection given by / i— > X[0,t)f- The second quantisation of an operator A G i3(L^(M)) is 
the operator in B(jr(L^(M))) given by 

oo 

r{A) :=ie0^^". 

n=l 

Then for A,B e B{L^{R)) we have r{AB) = r{A)r{B). This means that the one- 
parameter group of unitaries generated by Pt, i.e. exp(isPt), leads to a one-parameter 
group of unitaries r( exp(zsPt)) . Stone's theorem 1.1.2 asserts the existence of a selfadjoint 
operator A{t) on J^[L^{R)) such that for all s G R 

T{exp{isPt)) = exp (isA(i)). 

The operator A{t) is interpreted as the random variable that counts how many particles, 
i.e. photons, are present in the interval [0,t). 

The family of selfadjoint operators {A(t); t > 0} is commutative, i.e. the bounded oper- 
ators := (A(t) + i)^^ and Tl := (A(t) — i)^^ generate a commutative von Neumann 
algebra C to which the operators A(t) are affiliated. For / G L^(M) a state p on C is 
defined by 

p{C) := m),C^{f)), 

where V'(/) G •^(-^^(^^)) is the coherent vector defined as the exponential vector e(/) 
normalised to unity 

m :=exp(-i||/||2)e(/). 



1.2. STOCHASTIC CALCULUS ON FOCK SPACES 



9 



The state given by the coherent vector ip{f) is the state of the field in a laser beam. 
Spectral theory provides a probability space (r2,E,P) such that {C,p) is isomorphic to 
L°°(r2,S,P). As before we can realise the selfadjoint operators A{t) as measurable func- 
tions A(t), on n taking values in M. 

For si < ti < S2 < t-2 the characteristic function of the increments A(ii), — A(si), and 
A(t2)« — A(s2), for x, y G M is given by 

]E g»a;(A(ti).-A(si).)gi2/(A(t2).-A(s2).) ^ f g«a:(A(ti)a,-A(si)a,)giy(A(t2)u.-A(s2)u.)p^^^ ^ 

e-i^i' (e(/), r (e-(^*i "^^i )) V (e^2/(^'*2 -^.2 )) e(/)) = 

e-ll/lP(^e(/),e(e-(^*i-^n)e^^(^*2-^.2)/)^ = e<^'(^'^^"*^""^^ = 
exp ( r(e- - l)|/|2dA) exp ( f\e'y - iM^'dX), 

i.e. {A(t),}j>o is a process with independent increments. Moreover, the characteristic 
function shows that it is a Poisson process with intensity measure |/p(iA, where A stands 
for the Lcbcsgue measure. Summarizing, when photons are counted in a laser beam 
they arrive Poisson distributed. Since V-'(/) = ^(/)c(0) = W{f)^, wc can just as well 
study the commutative von Neumann algebra W{f)*CW{f), generated by the operators 
W{f)* A{t)W{f), equipped with the vacuum state to arrive again at a Poisson process 
with intensity measure \f\^d\. 

The above exposition shows that the quantum probability space (W, cf)) contains many 
processes that do not commute with each other. Hudson and Parthasarathy showed [54] 
that it is possible to extend the definition of Ito's stochastic integrals to deal with all 
these processes simultaneously. To this end we introduce annihilation A{t) and creation 
operators A*{t) by 

Ait) := 1 {B{ixio,t)) - iB{x[o,t))) and A*{t) := ^ {B{ixio,t)) + ^5(X[o,t))) • 

Let Mf stand for one of the three processes A{t),A*{t) or A(t), all restricted to the domain 
v. The stochastic integrals will be defined with respect to Mf. 



Let us write L^(]R) as the direct sum ((—00, t)) ©L^([t, 00)) , then .?^(L^(R)) is unitarily 
equivalent with ((—00, t))) ([t, cx)))) through the identification e(/) = e(/i))(gi 

e{f[t) with := /x(-oo,t) and /[^ := /X[t,oo)- We will also use the notation /j^^^) for 
fX[s,t) and omit the tensor product signs between exponential vectors. Furthermore, the 
algebra W = B {j^ [L"^ (R))) splits as W = Wj) ® W[t where Wt) = B {J^{L'^ {{-00, t)))) 
and W^t = ;B(jr(L^([t, 00)))). The following factorisation property [54], [74] makes the 
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definition of stochastic integration with respect to Mt possible 

{Mt - MMf) = <fs)){iMt - MMf[s,t))}<f[t), 

with {Mt — Ms)e(/[s t)) G jr(L^([s, t))) . We first define the stochastic integral for the so 
called simple operator processes with values in B (g) W, where B is an n-dimensional von 
Neumann algebra called the initial system. 

Definition 1.2.1: Let {Ls}o<s<t be an adapted (i.e. Lg ^ B ® Wg) for all < s < i) 
simple process with respect to the partition {sq = 0,si, . . . ,Sp = t} in the sense that 
Ls = Lgj whenever sj < s < Sj+i- Then the stochastic integral of L with respect to M on 
C^^V is given by [54], [74]: 




p-i 

LgdMs fe{u) := ^ (L, ./eK.))) ((M,.^, - M, .)e(n[,.,,.^^)))e(«[, .^J. 



By the usual approximation by simple processes the definition of the stochastic integral 
can be extended to a large class of stochastically integrable processes [54], [74]. The 
notation is simplified by writing dXt = LtdMt for Xt = Xq + Jq LgdMg. 

The following theorem of Hudson and Parthasarathy [54] extends the Ito rule of classical 
probability theory. 

Theorem 1.2.2: (Quantum ltd rule [54], [74]) Let Mi and M2 each be one of the 
processes A{t),A*{t) or A{t). Then M1M2 is an adapted process satisfying the relation: 

d{MiM2) = MidM2 + M2dMi + dMidM2, 

where dMidM2 is given by the quantum ltd table: 



dMi\dM2 


dA*{t) 


dA{t) 


dA{t) 


dA*{t) 











dA{t) 


dA*{t) 


dA{t) 





dA{t) 


dt 


dA{t) 






This theorem will prove to be much more useful in calculations than the actual defini- 
tion of the stochastic integral. It reduces hard questions regarding analysis to algebraic 
manipulations with increments dA*{t),dA{t) and dA{t). 

For a fixed a G [0,7r) we saw that the fields B{e'"^X[o,t)) define a Brownian motion. Note 
that from the definition of the creation and annihilation operator (1.2) it follows that 
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B{e''^X[o,t)) = ie-'°'A{t) - ie^"^*(t), i.e. dB(e^"x[o,t)) = ie~'°'dA{t) - ie''^dA*{t). There- 
fore, using Theorem 1.2.2, ((ii?(e*"x[o,t)))^ = dt, which is exactly what would be expected 
for a Brownian motion. Note that in the non-commutative theory we can also calculate 
products of increments of two non-commuting processes. For instance for ai,a2 G [0,7r), 
dB(e^"iX[o,t))dS(e^"2^[o,t)) = e'i'^^-'^^Ut. 

Let / be an element of L^(M). It is not hard to see (cf. [74]) that the Weyl operator 
= exp (^A*{f^^) — satisfies the quantum stochastic differential equation 

dWift)) = {mdA*{t) - l{t)dA{t) - \\f{t)\^dt}W{ft)). (1.7) 

Denote Nt := W{fYK{t)W{f) = W {f\)Y K{t)W {f\)) then it follows from Theorem 1.2.2 
and equation (1.7) that the Poisson process Nt satisfies [41] 

dNt = dA{t) + f{t)dA*{t)+J{t)dA{t) + |/(t)|^di. 

It easily follows that dN^ = dNt, which is exactly what would be expected for a Poisson 
process. 



1.3 Quantum optics 

Quantum optics deals with the interaction between quantum systems and the quantized 
electromagnetic field. In this section we will point out that in some Markovian approxima- 
tion this interaction is governed by a stochastic differential equation in the sense of Hudson 
and Parthasarathy. Indeed, the Markovian approximation is justified if the time scale of 
the field evolution can be considered to be extremely much faster than the time scale of 
the quantum system with which it interacts. In this way the field can be considered as a 
(non-commutative) noise acting on the quantum system. 

A convenient starting point (cf. [88], [45]) for quantizing the classical free electromagnetic 
field is the vector potential A{r,t) in the Coulomb gauge, i.e. V • A = 0. The magnetic 
field B and the electric field E are determined by S = V x A and E = — , respectively. 
The Hamiltonian of the free field is given by 

where Siiid Eq are the magnetic permeability and electric permittivity of free space, 
respectively. Then a Fourier expansion of A{r,t) is made and E{r,t) and B{r,t) arc 
expressed in terms of this Fourier expansion. After filling these expressions for E{r,t) 
and B{r,t) into the Hamiltonian one easily recognizes the Hamiltonian for an assembly 
of independent harmonic oscillators, see [88] and [45]. An assembly of quantum harmonic 
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oscillators is described by its algebra of creation and annihilation operators obeying the 
usual commutation relations, which justifies our description of the quantized electromag- 
netic field by a CCi?-algebra. In this description the electric field of a mode g G L^(M) 
is given by a polarisation vector times the operator i[A{g) — A*{g)^, which is just the 
operator B{g) from the previous section, see also [88], [45]. 

Let B stand for the algebra of observables of an n-dimensional quantum system, i.e. the 
elements of B are operators on C". Let W be the von Neumann algebra of observables of 
the electromagnetic field, i.e. elements of W are operators on some representation Hilbert 
space for the CCi?-algebra. In the theory of open quantum systems the field W is often 
referred to as the reservoir. The time evolution of the system B and the electromagnetic 
field W together is governed by the Hamiltonian H from quantum electrodynamics, cf. 
[38]. Using shorthand notation Hg+Hr := Hs'S>lr+'^s'^Hr,we can write this Hamiltonian 
as 

H = Hs + Hf + Hsr, (1-8) 

where the subscripts s and r stand for the system B and the reservoir W, respectively. 
The interaction between system and field is given by the interaction Hamiltonian Hgr- 

Let p and 7 be states on B and W, respectively. We can define a time evolved state on 
i3® W by 

x\Z) :=/9®7(exp(^if)Zexp(- ^F)), ZeB^W, t>0. 

This leads for times i > to the following reduced state on system B 

p\X):=x\X^lr), XeB,t>0, 

and it defines an operation T^ : B* ^ B* ioT:t>0 by T^{p) := p*. The evolution Tf is still 
rather complicated, for instance it is not a semigroup, i.e. we do not have Tt+g = TfTg for 
all s, t > 0. We proceed by assuming that Hgr is small compared to Hg + Hr- We separate 
the slow evolution generated by Hgr from the rapid evolution generated by Hg + Hr by 
transforming to the interaction picture, i.e. we define for Z G H ^ W and X E B 

X\Z) := x*(exp ( - |(^. + i^r))^exp + Hr))), 

p\X) := x\X ® 1,) = p*(exp ( - |i/,)Xexp (|^.)), and 
Hi := exp + Hr))Hgr exp {-'{{Hs + Hr)) . 

It easily follows that = j^x^{[Hl,Z]) for all Z e B <^W. Expanding this up to 

second order we find for all X G B 



dt 



= ^p®7{[Hl,X^lr])-^ x''[[HlrAHi,X®lr\])dt'. (1.9) 
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We may take j^p 'j[[Hl^.,X 1,]) = in the above equation by including a term 
Id ^{Hsr) in the system Hamiltonian. 

In the coming two paragraphs a heuristic exposition, common in physics, on the so-cahed 
Markov approximation to equation (1.9), cf. [27], is given. The approximation consists out 
of two steps. In the first step we replace x* by 7, and in the second step we replace 
p* in its turn by p*. The first step of the approximation is called the Bom approximation. 
It is justified if the reservoir is a very large system, i.e. its state is virtually unaffected 
by the coupling to B, and the coupling between system and reservoir through Hgr is very 
weak, i.e. at all times x* shows only deviations of order Hgr from an uncorrelated state, 
see also [27]. The second step is reasonable if the reservoir correlation times are much 
shorter than the time scale of the evolution of the system. Then the past history of the 
system, imprinted in the reservoir through the interaction, can not influence the present 
state of the system since in the reservoir it gets lost very quickly, see also [27] . 

After identifying the states p* and 7 with their density matrices, equation (1.9) reduces 
in the Markov approximation to 

^ = -^jy,r{[Hlr,[Hi,P'^^]])dt', (1.10) 

i.e. it is of the form dp* = L*{p*)dt. In many concrete examples, see [27], it can be 

shown, again by using the fact that the reservoir correlation functions decay extremely 
much faster than the time scale of the evolution of the system, that the generator L* is 
actually time independent. This means we end up with a time evolution T^{p) ■= p* that 
is a one-parameter semigroup, i.e. T^Ts = Tf^s for all t,s > 0. The semigroup property 
reflects that there are no memory effects present. The equation dp^ = L{p*)dt is called 
the master equation. The following theorem of Lindblad characterizes the generator of a 
one-parameter semigroup of operations on ;S = M„(C). 

Theorem 1.3.1: (Lindblad [67]) Let {Tt}t>Q be a semigroup of completely positive 
identity preserving operators on M„(C) with generator L. Then there exist a self-adjoint 
element H G M„(C) and elements Vj G M„(C) for j = 1, 2, . . . , fc with k < n^, such that 

k 

L{X) = i[H,X]+Y,V*XVj - -{V;V,,X}, X G M,(C), (1.11) 

J=l 

where {X, y} stands for the anti-commutator XY + YX. Conversely, every operator L 
of this form generates a semigroup of completely positive identity preserving operators. 

Lindblad's result [67] is actually more general, it is valid for norm-continuous semigroups 
on the algebra of bounded operators B{M.) for some, possibly infinite dimensional, Hilbert 
space M. The commutator with H describes the evolution generated by some system 
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Hamiltonian and the part with the Vj's deseribes the decay into the field. The Vj^s are 
determined through equation (1.10) and the expression for Hgr in quantum electrodyna- 
mics. 

Another, more rigorous, approach to the master equation is via the weak coupling limit. 
In this limit the Hamiltonian H of equation (1.8) is replaced by 

= Hg + Hj. + XHgr- 

While the coupling constant A is sent to 0, the time variable t in the interaction picture 
has to be scaled to r := to compensate the slower decay of the system, i.e. for X e B 
and Z e B^yV we define 

pHX) := xliX ® Ir), where Xa(^) - P ® l{U^*ZU^) , with 
:= exp {j{Hs + Hr)) exp ( - '-I-H^). 

Davies [32] showed that under some technical assumption the time-evolution T^, given by 
T^'*{p) := p^, uniformly converges to a semigroup Tt when A goes to 0. 

Let us assume that the electromagnetic field is initially in the vacuum state, i.e. 7 = ^. 
Let Tt be the semigroup describing the irreversible evolution of the system in the weak 
coupling limit. Its generator L can be written in the form of equation (1.11). Accardi, 
Frigerio and Lu showed, see [2], that when A goes to 0, converges in distribution to 
the solution Ut of the quantum stochastic differential equation 

dUt = {VjdA*{t) - V*dAj{t) - {iH + ^V*Vj)dt}Ut, 

where repeated indices are summed. In this context the weak coupling limit is also called 
the stochastic limit. Just to support this result, let us check that Id (p{UtX Ir^^t) 
is the semigroup Tf. Using the continuous tensor product structure of the Fock space on 
which Ut is defined, it is not hard to see that the expression Id ® (p{UtX (gi IrUt) defines a 
semigroup, i.e. we only have to show that it is generated by L. Since the noises Aj and Aj 
are independent for different j and vacuum expectations of and Aj are zero, it follows 
from Ito's rules that 

did (t){U;X l,.Ut) = Id (>5 ^{d{U;X l,.Ut)) = 

Id (j){{dU;)X ^ IrUt) + Id «) 0(f/;X ^ IrdUt) + Id «) (j){{dU;)X ^ IrdUt) = 
(iHX - ]^V*VjX)dt + {-iXH - ^XV*Vj)dt + VjXVjdt = L{X)dt. 

This result was already obtained by Hudson and Parthasarathy in [54]. 

In this thesis we will always take the semigroup and its corresponding quantum stochastic 
differential equation, describing the interaction of the system B and the electromagnetic 
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field W in tlie vacuum, as tlic starting point. The Markov approximation or the weak 
coupling limit is a very good approximation for many phenomenon in quantum optics. 
Barchielli [7] was one of the first to see the relevance of quantum stochastic calculus for 
quantum optics. See for instance [8] for a description of the electron shelving effect with 
quantum stochastic calculus. In [76] the Mollow spectrum of fluorescence of a driven two 
level atom is derived using quantum stochastic calculus. 



1.4 Filtering and control, outline of results 



In this thesis a quantum system B in interaction with the electromagnetic field W is 
studied in the weak coupling limit. For simplicity we always assume B to be M„(C) for 
some n G N. Next, an observable Yf of the field W is measured continuously in time. 
For instance in Chapter 2 the photons emitted into the field by a laser driven two-level 
atom are counted continuously in time. Since the system B interacts with the field W, we 
gain information about the system B when the observable Yt in the field is measured. For 
instance, when the field is in the vacuum state we can infer immediately after a photon 
appears in the field and is counted that the two-level atom is in the ground state. The 
central question in this thesis is how to condition the state of B continuously in time on 
information gained by measuring some observable Yt in the field. 

Chapter 2 introduces the subject by analyzing a photon counting experiment, i.e. the 
observable Yt is the photon number operator A{t) of the field. The system H is a two-level 
atom driven by a laser. We explicitly write down the solution of the quantum stochastic 
differential equation using Maassen's integral-sum kernels. The conditioning is done by 
sandwiching with the projection corresponding to an observed event. In this way a whole 
family of evolutions for the two-level atom is obtained, i.e. for every possible observed 
event E in the photon counter there is a map that represents the reduced evolution of 
the two-lcvcl system conditioned on this event E. We prove that this family of maps 
satisfies the axioms of the processes studied by Davies [31], [33] in the late sixties and 
early seventies. We use his theory for these processes and the explicit solution, in terms of 
Maassen's kernels, of the quantum stochastic differential equation to calculate the jump 
operators that describe the evolution when a photon appears and is detected in the field. 
This leads to a continuous evolution of the two-level atom interrupted by jumps at the 
moments at which photons are detected. In quantum optics such an evolution is known 
as a quantum trajectory [27]. We use the trajectory evolution to show that the photons 
are detected according to a renewal process. 

In Chapter 2 the explicit solution of the quantum stochastic differential equation is used, 
i.e. we leave the strength of Ito's calculus for increments unused. In Chapter 3 the condi- 
tioned evolution of the system B is described infinitesimally, fully exploiting the quantum 
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ltd calculus summarized in Theorem 1.2.2. Some of the ideas of Chapter 2 directly carry 
over to this description. For every time t > there is a measure space (r2j,St,Pj) with 
Ut the set of all possible paths of the observed process up to time t. In both chapters we 
arc interested in the consistency of the family of measures {Pt}t>o, i-e. for all s > t and 
E ^ Tit IPs(-E') = ^t{E). Kolmogorov's extension theorem then states that this family 
extends to a single probability measure P on the paths observed up to infinity. 

The conditioning by sandwiching with projections in Chapter 2 is not that easily carried 
over to the infinitesimal setting of Chapter 3. In the third chapter this is done using 
the decomposition of a von Neumann algebra A, equipped with a normal state x aiid 
represented on some Hilbert space H, over its center C := {C G A; CA = AC, "iA G A}. 
Theorem 1.1.1 states that there exists a probability space (0,S,P) such that (C,x) = 
L°°(J1,E,P). The decomposition of A over C described below (see [56] for proofs) is an 
extension of this result. 

The Hilbert space H has a direct integral representation H = j^M.^f'{duji) [56] in the 
sense that there exists a family of Hilbert spaces {Ht^jt^eQ and for all ■0 G H there exists 
a measurable map u ^ '^^ & M.^ such that 

(V', 4>)= ! {ipuj, ^J)n<^), V', </> e 

The von Neumann algebra A has a central decomposition A = Au)'^{duj) [56] in the 
sense that there exists a family {»4t^,}u;en of von Neumann algebras with trivial center, 
called factors, and for all A G v4 there exists a measurable map u i— > A^ G Ai^ such 
that {Ail))ij = A^ipi^ for all ^ G H and for P-almost all a; G O. The state % on ^ has a 
decomposition in states Xuj [56] such that for all A G .4 its expectation is given by 

X{A) = [ xUA^Mdu^)- 
Jn 

We think of the state x a-^d an arbitrary operator A E A as maps %, : UJ i-^ Xui and 
A, : u A^, respectively. 

Let t>0 and let p be a state on B. Define a time evolved state x^i^) '■= p(^4'{UtZUt), Z G 
B<^yV in the interaction picture, with Ut the solution of the quantum stochastic differential 
equation describing the interaction between B and W. Define a commutative algebra Ct 
as the algebra generated by the process {Ys}s>o up to time t, i.e. Ct is generated by the set 
{Ys', < s < t}. Because of the consistency and Kolmogorov's extension theorem there 
exists an increasing family of (X-algcbras {'Et}t>o such that (Ct, x*) — L°°{Q, St,Pp) where 
Q, is the set of all paths of the observed process in the field {ys}s>o up to infinity. 

We define a subalgebra At of B^W diS the commutant of Ct, i.e. At := C't, and we restrict 
X* to this algebra At- Now Ct is the center of At and we can decompose the state x* over 
this center. In this way we get a random state x. from U to the states on At which for 
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each t > is measurable with respect to Sj. Since B <E> lyv is a subalgebra of At for all 
t > and since (X® lyv). is the constant map a; i— >■ X, we obtain a random reduced state 
from CI to the states on B given by 

The state is the state of the system B after t seconds of time evolution conditioned on 
having observed the first t seconds of the path uj. The state p*^ on B after t seconds of 
time evolution conditioned on some event E £T,t is then given by 

pUx) = I pl{x)¥p{dw), xeB. 

JE 

The goal of Chapter 3 is to derive a stochastic differential equation for the state in 
which the stochastic part of the equation is determined by the observed process {Fs}s>o. 
As an example, the stochastic differential equation for the two-level atom of Chapter 2 
when photons are counted in the field is given by 

dpiiX) = pi{L{X))dt + {jr^^ - pH^)) {dNt - Pi{m)dt), XeB, (1.12) 

where L is the Lindblad generator of the semigroup evolution when we would not be 
counting photons in the field, Nt is the random variable from J7 — > N that represents 
the number of photons counted up to time t, and J' is the jump operation of Chapter 2. 
The stochastic differential equations derived in Chapter 3 are called Belavkin equations 
[15], [17]. They are known in quantum optics under the somewhat misleading name of 
stochastic Schrddinger equations [27] since they are actually a stochastic variant of a 
master equation. 

Apart from the conditioning for which we use the central decomposition, our derivation 
of the Belavkin equation is rather close to the original derivation of Belavkin [17]. How- 
ever, instead of using the explicit construction of the quantum conditional expectation 
encountered in that derivation, our proofs try to exploit its characterizing properties. 
Furthermore, we work in the interaction picture, which makes the derivation more acces- 
sible. 

Let us return to the Belavkin equation, see the example in equation (1.12). It is an equation 
for the expectation pi{X) of some system operator X conditional on the observed process 
Yt in terms of dt and dYf. We only have acces to the process Yt in the field and since 
the field is coupled to the system B we indirectly gain information on X. This is similar 
to the situation in classical filter theory, developed mainly by R.L. Stratonovich in the 
late 1950s, cf. [83], [57]. There one is interested in some system process Xt but only has 
access to a process Yt which is the system process polluted with some noise. In classical 
filtering theory one derives a non-linear stochastic differential equation for the conditional 
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expectation of Xt on the observed process Yt in terms of dt and dYt, called the filtering 
equation. The Belavkin equation is a non-commutative analogue of this equation and is 
therefore some times called the quantum filtering equation. 

In Chapter 4 a similar setup with a system B in interaction with the field W is being 
studied. The problem is to control the state of the system B, i.e. we start with an unknown 
initial state p of the system B and then we try to control it in a way that depends on 
observation of some process Yt in the field, aiming to keep the state of the system B as 
close as possible to p. First the case is studied where the stochastic differential equation 
governing the interaction between B and W contains only commutative noise terms. This 
special situation is called essentially commutative [63]. 

To control the state it is first evolved over a time interval r according to the Belavkin 
equation for the measurement we are performing in the field, i.e. 



where the tilde is meant to indicate that pi is the state of B before control. In the time 
interval r some measurement result has been obtained for the process observed in the field. 
We then correct with a unitary that depends on this result i.e. denoting the state of 
B with its density matrix, we have 



We derive a stochastic differential equation for the controlled state evolution depending 
on the stochastic measurement process. It turns out that in the essentially commutative 
case it is possible to freeze the state evolution, i.e. dp* = 0. Such a control scheme is said 
to restore quantum information [51]. 

However, the interaction between the system B and the field can in general not be treated 
in an essentially commutative way. For instance, the interaction causing spontaneous 
decay of a two-level atom into a vacuum field is given by a quantum stochastic differential 
equation in which the two non-commuting noises Bf = S(x[o,t)) and B^^'^ = B{ix[o,t)) of 
Section 1.2 appear. The strategy when we are not in the essentially commutative case is to 
manipulate the field state such that the quantum stochastic differential equation resembles 
the essentially commutative situation more and more. This can be done by replacing the 
vacuum state of the field by a squeezed state. 

In the vacuum state the variances of Bf and B^^'^ are both equal to t. In a squeezed state 
the variance of one of these noises is decreased while the other one is increased, since we 
still have to satisfy Heisenberg's uncertainty relation. The noise with the large variance 
has the biggest disturbing influence on the evolution of the system B, therefore this noise 
is being observed in the field and the control scheme is based on the measurement results 
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of this observation. It turns out that in the hmit for squeezing to infinity the situation 
where quantum information can be restored is refound. 



Chapter 2 



The Davies process of resonance 
fluorescence 
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Abstract^ 

Starting point is a given semigroup of completely positive maps on the 2x2 matrices. This 
semigroup describes the irreversible evolution of a decaying two-level atom. Using the integral-sum 
kernel approach to quantum stochastic calculus we couple the two- level atom to an environment, 
the electromagnetic field. The irreversible evolution of the two-level atom stems from the reversible 
time evolution of atom and field together. Mathematically speaking, we have constructed a Markov 
dilation of the semigroup. 

Next step is to drive the atom by a laser and to count the photons emitted into the field by the 
decaying two-level atom. For every possible sequence of photon counts a map is constructed that 
gives the time evolution of the two-level atom inferred by that sequence. This family of maps 

^This chapter is an adapted version of [24]. 
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forms a so-called Davics process. In his book Davies describes the structure of these processes, 
which brings us into the field of quantum trajectories. Within our model we calculate the jump 
operators and we briefly describe the resulting counting process. 



2.1 Introduction 



In this paper we want to illustrate that quantum stochastic calculus together with the 
processes studied by Davies in his book [33] , and explained in his paper with Srinivas [82] , 
form a suitable mathematically rigorous framework for doing quantum trajectory theory 
[27]. As an example we consider here the case of resonance fluorescence. 

Our starting point is a semigroup of transition operators {Tt}t>o on the algebra M2 of 
all 2 X 2-matrices. This semigroup describes the irreversible evolution of a spontaneously 
decaying two-level atom in the Heisenberg picture. By coupling the atom to a quantum 
noise, we construct a stationary quantum Markov process having precisely these transition 
operators. If we impose the requirements that the external noise be a Bose field, and the 
quantum Markov process be minimal, then the latter is uniquely determined. It is called 
the minimal Bose dilation of {M2,Tt,g) [60], where g is the ground state of the two-level 
atom. 

Since this dilation is uniquely determined, any other reversible dynamical model which 
couples (M2 ,Tt,g) to some Bose field necessarily contains this Bose dilation as a subsys- 
tem. Therefore, without deriving our model from an explicit Schrodinger equation (by 
performing a Markovian limit) we may safely assume it to be a physically correct way to 
describe the interaction of the two-level atom with the electromagnetic field. 

We will couple the two-level atom to the electromagnetic field by using quantum stochastic 
calculus [74] , [72] . We use a version of quantum stochastic calculus based on integral-sum 
kernels [69], [68], [72], which has the advantage that we have an explicit construction for 
the solution of the quantum stochastic differential equation with which we will describe 
the coupling of atom and field. Having this explicit construction in our hands is important 
for doing the actual calculations we encounter later on. 

To be able to discuss resonance fluorescence a dilation with two channels in the electro- 
magnetic field is used. On one of them wc will put a laser state to drive the two-level 
atom. We will call this field the forward channel and the other one the side channel. Then 
photons are counted in both channels. We need the side channel, because we know that 
there all detected photons are fluorescence photons. In the forward channel a detected 
photon could just as well be coming directly from the laser. 
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For every event that can occur in the photon counters a map is constructed that gives 
the evolution of the two- level atom inferred by that event. We will see that the family of 
maps we obtain, fulfills the axioms for the processes discussed by Davies [33]. We have 
constructed the Davies process of resonance fluorescence. 

Using the structure theory for Davies processes [33] we can decompose the process into 
its trajectories [27]. Within our model we calculate the expression for the jump operators 
and for the time evolution in between jumps. Note that a jump in the system occurs the 
moment we detect a photon, since our knowledge concerning the system changes. Using 
the above apparatus we show that the resulting counting process in the side channel is a 
so-called renewal process. 



2.2 The dilation 



Let M2, the algebra of 2 x 2-matrices, stand for the algebra of observables of a two- 
level atom. On this algebra we are given a (continuous) semigroup {Tt}t^o of completely 
positive maps. This semigroup describes the, generally irreversible, evolution of the two- 
level atom. Lindblad's Theorem [67] then says that T^ = exptL where L : M2 — > M2 can 
be written as: for X G M2: 

k 

L{X)=i[H,X]+J2v*XVj - -{V;Vj,X}, (2.1) 

j=l ^ 

where the Vj and K are fixed 2 x 2-matrices, K being Hermitian. In this paper we will 
restrict to the simpler case where H = and there are just two VjS. This means there 
is dissipation only into two channels, the forward channel described by Vf, and the side 
channel described by V^. We choose Vf and Vg such that: 

^=(10)' ^/ = ''/^' ^s = >^sV, |/^/|2 + |/^,|2 = l. 

This exactly gives the time evolution for spontaneous decay to the ground state of the 
two-level atom into two decay channels, where the decay rates are given by and |/«sp. 

We want to see this irreversible evolution of the two-level atom as stemming from a 
reversible evolution of the atom coupled to, in this case, two decay channels in the field. 
So let us first construct the algebra of observables for these channels in the field. Let 
T be the symmetric Fock space over the Hilbcrt space L^(]R) of square intcgrable wave 
functions on the real line, i.e. := C (B 0^j^ -L^(M)®''"'. The electromagnetic field is 
given by creation and annihilation operators on J^, generating the algebra of all bounded 
operators. We need two copies of this algebra, which we denote by Wf, which will be the 
forward channel, and Wg, which will be the side channel in the field. 
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The evolution over a time t of a free field is given by the second quantization of the 
left shift, i.e. the second quantization of the operator on L^(R) which maps /(•) into 
/(• + We denote the second quantization of this operator by St- This means that in the 
Heisenbcrg picture we have an evolution on Wj® Ws mapping X into {S^ <Si Sf)X{St<Si St) 
(= (5_t S-t)X{St «) St)), also denoted by Ad[St «) St]{X). 

The presence of the atom in the field introduces a perturbation on the evolution of the free 
field. We let this perturbation be given by a certain family of unitary operators {Ut}t^m 
on 'C? T (^T, to be specified later, that forms a cocycle with respect to the shift S'j S'j , 
i.e. for all s G M : Ut+s = {S-s ^ S-s)Ut{Ss 'Si Ss)Us- Given this cocycle, we let the 
time evolution of the atom and the field together be given by the following one-parameter 
group {Tt}teR (i-^- the evolution is now reversible) of *-automorphisms on M2'S>Wf'S>Ws- 
for all X eM2SWfS Ws'. 



Tt{X) 



Ut\S-t ® S-t)X{St S St)Ut if t > 
■1 
t 



® S-t)U-tXUZ^ {StSSt) if t < 



The solution of the following quantum stochastic differential equation [54], [74] provides 
us with a cocycle of unitaries with respect to the shift: 



dUt = {VfdA}^ - VjdAf^t + VsdAl^ - V*dAs,t - -V*Vdt}Ut, Uo = I. 



(2.2) 



In the next section we will give an explicit construction for the solution Ut of this equation. 
It can be shown ([54], [40], [69], [74]) that if the cocycle satisfies equation (2.2) we have 
constructed a so-called quantum Markov dilation (M2 ® W/ Wg, {T}teM., id(8) 0®^) of the 
quantum dynamical system {M2,{Tt}t>o,g) [59], [60], where (/) is the vector state on W/,s 
given by the vacuum vector. This means that the following dilation diagram commutes 
for alH > (and that the resulting quantum process is Markov): 



Mo 



Tt 



Mo 



Id®7®2 



Tt 



M2 (g) W/ (g) 

i.e. for ah X e M2 : Tt{X) = (Id ® (f^(^x ® I^^)) . 



M2 (g) W/ (g) Ws 



(2.3) 



The diagram can also be read in the Schrodinger picture if we reverse the arrows: start 
with a state p of the atom M2 in the upper right hand corner, then this state undergoes 
the following sequence of maps: 

p ^ p (g <^®2 ^ (p (g <^®2) o ft = ft4p (g> ^^^) ^ Tr^,^^^ {ft^p (g cl>^^)) . 

This means that at t = 0, the atom in state p is coupled to the k channels in the vacuum 
state, and after t seconds of unitary evolution we take the partial trace over the 2 channels. 



2.3. GUICHARDET SPACE AND INTEGRAL-SUM KERNELS 

2.3 Guichardet space and integral-sum kernels 
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Let us now turn to giving the explicit construction for the solution of equation (2.2). For 
this we need the Guichardet space [52] of M, which is the space of all finite subsets of 
M, i.e. O := UneN^"' where := W <^ ^5 I'^'l = '^}- Let us denote by A„ the Lebesgue 
measure on R". If, for n G N, we let jVi '■ 1^"' denote the map that maps an ra-tuple 

(ti,t2, . . . ,tn) into the set {ti,t2, ■ ■ ■ ,tn}; then we can define a measure Hn on $7^ by: 
IJLn{E) := ^A„(j~^(£^)) for all E in the sigma field of Qn induced by jn and the Borel 
sigma field of M". Now we define a measure /i on $7 such that /x({0}) = 1 and = Hn 
on Qn- This means we have now turned the Guichardet space into the measure space 

The key to constructing the solution of equation (2.2) is to identify the symmetric Fock 

space with the space of all quadratically integrable functions on the Guichardet space 
L^(r2,/i). To see this identification note that L^(r2„,/i„) is, in the canonical way, unitarily 
equivalent with the space of all quadratically integrable functions on R" invariant under 
permutations of coordinates, denoted Lgy-^{W^). It is now obvious how to identify JF = 

c e e~ 1 Lly^iR^) with L^{n, ^i) = c® e^=i L^n^, fin). 

For every / G L^(R) we define the exponential vector e(/) £ in the following way: 

e(/) := 1 © / © ^/^^ © Te-^*^^ ® ^^^^ linear span of all exponential vectors 

forms a dense subspace of JT. For every / G L^(R) we define the coherent vector tp{f) to 
be the exponential vector of / normalised to unity, i.e. V'(/) = 6xp(— ^||/p)e(/). Under 
the above identification of J- with L'^(Q, fi), the exponential vector (of an / G L^(R)) e(/) 
is mapped into an element of L'^{Q,fi) which we denote by vr(/) and which is given by: 
7r(/) : Q ^ C : a; i-^ Ylseu /('^)' where the empty product Hse© f^^) defined to be 1. 
We will often choose for / the indicator function of a certain interval / C R, which we 
denote by xi- This is the function which is 1 on 7 and elsewhere. 

We are now ready to start the construction of the solution Ut of equation (2.2). Define the 
integral-sum kernel of Ut (name will become apparent in a minute) to be the map ut that 
maps four disjoint finite subsets of R, o"/, as,Tf,Ts (where / and s stand for "forward" and 
"side") to the following 2 x 2-matrix, where we write cr/UfTsUr/Urs also as {ti,t2, ■ ■ ■ , tk} 
such that ii < ^2 < • • • < ^fc and A; G N: 



utiaf,as,Tf,Ts) :=7r(x[o,t])(c^/ UcJs Ur/ Urs)exp(- 



^ - tfc 



V*V)VkX 



exp(-^^^y*F)Ffe_i ...V^ exp(-|y*F), 
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where 



Vf 


if 


tj 


£ af 


-Vf 


if 


tj 


G Tf 


Vs 


if 


tj 


G as 


-v: 


if 


h 


G Ts 



Then we have the following theorem of Maassen, see [69] , [68] : 

Theorem 2.3.1: After identifying C'^®J^®J^ with L^a (^^ x fi, jj), the space of all square 
integrable functions on CI X CI with values in C^, the solution Ut : L'^2{Cl X CI, /j, X /j,) ^ 
L'^^i^ X CI, ji X ji) of equation (2.2) is given hy: 

{Utf){ujf,UJs) = y2 / Uti'^f,Crs,Tf,Ts)f{iuJf\af)UTf,{uJs\as)UTs)dTfdTs. 



Ofdujf 

asCu). 



Now we have an explicit expression for the time evolution Tf = Ad[?7f], where Ut is given 
by {St St)Ut if t > and Uzl{St <8) St) if i < 0. The family {Ut}tes. forms a group 
of unitary operators on <C'^ ® ® T describing the time evolution of the two-level atom 
and the two channels in the field together. Stone's Theorem says that there must be 
a Hamiltonian associated to this time evolution. This Hamiltonian has been calculated 
recently [50], [49]. 



2.4 The Davies process 



We now return to the situation in figure 2.3. We wish to make some changes in this diagram 
and for this we need to introduce some more notation regarding Guichardet spaces. Let 
/ C M be an interval. Then the Guichardet space of I is the set Cl{I) = Un^o ^n(^)) where 
Cln{I) = {a d I; |cr| = n}. In a similar way as for Cl, which is r2(M), we can provide these 
sets with a measure structure: {Cl{I), S(/), //). Given a subset E of Ct{I) in the sigma field 
S(/), we can construct the projection M^^ : L'^{Cl, fi) — L^{Cl, fi) : f Xe/- 

Let / be [—t,0), then the events in S([— 1,0)), which we abbreviate to Sf, are events in 
the output field of the atom up to time t. Remember that the evolution of the free field 
was given by the left shift and that the atom is sitting in the origin. Since the Guichardet 
space representation corresponds to the photon number picture, we can give concrete 
interpretations to the subsets in Ef. For instance, the subsets r2„([— 1,0)), correspond to 
the events " there are n photons in the output of the atom into this channel of the field up 
to time t" . 

Now back to the situation in figure 2.3. Suppose we have been observing the output in 
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the forward and side channel of the atom up to time t with two photon counters. Then we 
are given two events Ef and Es in Ej. Since we know the outcome of the measurements 
we have to change the time evolution of the two-level atom, i.e. we have to project onto 
the observed events (see also [9]). This is summarized in the following figure: 



M2 



M2 (g) W/ (g) Ws 



Tt 



M2 



M2 (g) W/ (g) Ws 



where we have suppressed the capital letters M in the projections. The map £Q[Ef,Es] : 
M2 — ^ M2 : X I— ^ Id (g 4)®'^{Tt{X (g XEf <g XEs)) is the unnormalized time evolution of 
the two-level atom in the Heisenberg picture given that we see event Ef in the output 
of the forward channel and event Eg in the output of the side channel. If we are given 
a state on M2, i.e. a 2 x 2 density matrix p, then the probability of seeing event Ef in 
the forward channel and Eg in the side channel after t seconds of observation is given by: 
K\.iEf,Es)]=T,{p£l[Ef,Eg\{I)). 

The setting is still not complete for describing resonance fluorescence. Since we are not 
driving the atom, both the forward and the side channel fields are in the vacuum state, 
at most one photon can appear in the output. We change this by putting on the forward 



channel a coherent state with amplitude z G C, defined by: 7^ 



W 



C : X 



exTp{—t\z\^){TT{zx[Q^t\)->^''^{^X[o,i\))- Note that 70 is the vacuum state. Putting a coherent 
state on the forward channel mimics a laser driving the atom. We have suppressed its 
oscillations for the sake of simplicity. Now we are ready to do resonance fiuorescence, i.e. 
the diagram has changed into: 



M2 



M2 (g W/ (g Ws 



Tt 



M2 



M2 (g W/ (g Ws 



where the map Sl[Ef,Es] : M2 ^ M2 is now defined by Sl[Ef,Es]{X) := Id (g 7^^ (g 
7o(lt(X (g) XEf <g XEs))- It describes the unnormalized time evolution of the laser-driven 
atom given that wc see event Ef in the output of the forward channel and event Eg in the 
output of the side channel. Given a state p of the atom, the probability of seeing event 
Ef in the forward channel and Eg in the side channel after t seconds of observation is now 
given by: P* [(-E/, Eg)] = Tr[p6l[Ef,Es]{I)). To make the notation lighter we suppres the 
z in f * in the following. 



Since L^(r2,S,/x) (g L^(ri,i;,//) is canonically isomorphic to L^{Q, x i^, S (g) S, x /n) we 
can simplify our notation even a bit further. By identifying these spaces we can write 
f *[£■/, Eg] = S*[Ef X Eg], where the righthandside is defined by: for aU E G St (g Sj, X G 
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M2,t > : £^[E]{X) := Id^ ^fz^fl{ft{X xe)) , where j^tfi is an abbreviation for 7zj(8)7o- 
We will now study the properties of the family of maps we defined. 

Theorem 2.4.1: The family of maps {£*[E]}t>o,Ee'St®^t satisfies the axioms of a Davies 
process, [33]: 



1. For all t >0 and £^ G Sj, £^[E] is completely positive. 

2. For all t > and all countable collections of disjoint sets {En} in (8) 
and for all X e M2 : S'[\Jn En] {X) = En [En] {X) . 

3. For allt>0 we have f * ' 



n{[-t,o))xn{[-t,o)) (/) = /. 



4. For all X e M2 : limt_^o £^ 



n{[-t,o)) X n{[-t,o)) {x) = x. 



5. For all t,s>0 and G (g) S^, -F G St (g) and all X e M2 we have: 

£^ [F] o £^ [E] (X) = £'+* [F -sQE]{X), 

where F — s G t — s, —s) (gi t — s, — s)) is given by: 
F- s = {(// -sjs- s); iff, fs) G F)} and U is defined by: 
AOB = {{uf U af,LOs U (Ts); {ojf,LJs) € A, {erf, as) G B}. 



Proof. The only point where there is really something to prove is point 5. Let us first 
introduce some short notation which we shall only use in this proof. Let 7r{zt,0) denote 
7r(zx[o,t]) 'X' vr(0) and denote St St just by St- Further we use the notation (Jt{Us) for 
S-tUs 'St. Then for all X G M2, s, t > 0, £^ G and F G St we have: 

£^F]o£'[E]{X) 
exp(-(s + t)|z|2) 

'^7r{zt, 0),ft ((7r(z„ 0),tiX xeMz,, 0)> Xf)^^^ 0)) 
'n{zt, 0), U;{7r{zs,0),fs{X ® xeMzs,0)) ® XF+tUt7r{zt, 0)^ 

[TT{zt, 0), U*{S^tTT{Zs,0), S^ttiX ® XE)StS^tTT{Zs,0)) ® XF+tUtTT{zt, 0) ) = 



= f*[F] (7r(z„O),r,(X0X£W^.,O)) exp(t|zn = 



TT 



{zt, 0), t/; (5_t7r(z„ 0), cJt(?7,)*X XE+t+sCTt{Us)S-t7T{zs, 0)> ® XF+tC/t7r(zt, 0) 



Now we use the cocycle identity and the continuous tensor product structure of the sym- 
metric Fock space to obtain: 

fxp^L°(f+f)|iy) = {A^t+s,0),{at{Us)UtrX(^XF+tuE+t+s<^t{Us)Ut7r{zt+s,0) 



Tr{zt+s,0), Ut+sX ® XF+tOE+t+sUt+sT^izt+s,^) ^ 

£^+\F - s(jE]{X) 



,T^{zt+s,^),Tt+s{X ® XF-s\JE)'^{zt+s,^) 



exp(-(s + t)|z|2) 
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□ 

Define maps Yt : M2 ^ M2 : X ^ £:*[{(0,0)}] {X). They represent the evolution of the 
atom when it is observed that no photons entered the decay channels. Then for all t,s >0 

YtY,=S'[{M}] oS'iWM] =f*+1{(0,0)}-sU{(0,0)}] =f*+^[{(0,0)}] =Yt+,, 

i.e. the family {lt}j>o forms a semigroup. Now observe that for X G M2 and t>0 

Yt{X) = [({0}, {0})] (X) = Id 7,, 70 {ft{X xm ® xm)) = 
{TT{zt) 7r(0), U^X xm ^ XmUtTT{zt) 7r(0)> cxp(-t|z|2) = 
{Ut7r{zt) ® 7r(0)) *(0, 0)X(C/t7r(^t) ® 7r(0)) (0, 0) exp{-t\z\^), 

i.e. we can write Yt{X) = B^XBt where Bt is given by 

i?t= (exp(-^i|z|2)(C/i7r(zt)®7r(O)))(0,0). (2.4) 
Using Theorem 2.3.1 it follows that Bt is the following semigroup of contractions: 

Bt = exp ( - i(|z|2/2 + V*V + 2zVf*)t), 

If Yt can be written as B^ ■ Bt for some semigroup of contractions Bt, the Davies process 
£* is called ideal, see [33]. 

Lemma 2.4.2: The Davies process has bounded interaction rate i.e. there exists a 
constant K such that for allt>0 

s^[n[-t,o) X j^H,o)\{(0,0)}] (/) < tKi. 

Proof. From Theorem 2.4.1 point 2 it follows that 

[n[-t, 0) X n[-t, O)\{(0, 0)}] (/)=/- B^Bt, 
with Bt as in equation (2.4). Using Theorem 2.3.1 we find 

Bt = exp(-ltl.P) (^"pH) 2.7c,(exp(-i) - l)^| _ 

Therefore 

£'[n[-t,o)xn[-t,o)\{{$,$)}] (/) = 

/ i_e-*(N'+i) 2z^/e-*(l^l'+5)(l-e-|) \ 

\^2zK/e-*(N'+^)(l-e-5) l-e-*N'(l + 4|z|2|/^;|2(l-e-f)')J ■ 
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Since (5* 0) x Q,[—t, O)\{(0, 0)}] is completely positive this matrix has to be positive, 

i.e. its eigenvalues are positive. Let us denote them by and E- such that E^ > > 0. 
Define a := exp(— and f3 := exp(— |). Then for alH > 

a>0, l-/3>0, < 1 - a < \z\'^t, < 1 - a/?^ < + 

Therefore, for alH > 

f*[r2[-t,0) X f)[-t,O)\{(0,0)}](/) < E+I < 
{E+ + E^)I = Tv(^*[O[-t,0) X f][-t,O)\{(0,0)}] (/)) = 
((1 - a/32) ^ _ ^) _ 4^a\z\^\^^f\^{|3 - lf)l < 
{{l-a0^) + {l-a))l < (2|z|2 + l)il. 

□ 



2.5 Quantum trajectories 

In the seventies Davics studied the structure of what we now call Davies processes [33]. Let 
us first state his results, as far as relevant, in the context of the process we are studying. 

Lemma 2.5.1: (Davies [33]) Given any ideal Davies process with bounded interaction 
rate, as defined in the previous section, we have existence of the following limits: 

J/ :=lim^f*[Oi[-i,O),{0}] and J, := lim [{0}, J^i [-i, 0)] . 

These completely positive maps represent the action we have to apply on the two-level 
atom the moment we see one photon appear in the forward, respectively side channel. They 
are the jump operations for these channels. We will explicitly calculate these limits later 
on, but first we turn our attention to decomposing the Davies process into its trajectories 
[27]. For this we need the following definition. 

Definition 2.5.2: Let Yt : M2 M2 be the maps from the previous section, i.e. Yt = 
£^*[{0},{0}] and let J'f and J7s be the maps from lemma 2.5.1. Let Uf and lVs be disjoint 
elements of n[—t, 0) and denote cof U lOs also as {ti,. . . ,tk} where —t<ti<t2<---< 
ifc < for a A; e N. Then we define: 

WY,jf,jM,u;s) := Yt,+tJ''Yt,.t,J''^ . . . j'^Y_t^, 

where JT"*' denotes Jg if U € and J'f if ti e ujf. 
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Since Yt is the time evolution of the system when, both in the forward and the side 
channels, no photons are detected and J'f and J'g are the jump operations that we have 
to apply when a photon in the corresponding channels appears, it is clear that the string 
of maps Yt-^^tJ'^^Yf^^ti'J^^ ■ ■ ■ J^'^Y-t^ represents the trajectory of an observable X in M2 
when we find the outcomes 1^/ in the forward and lOs in the side channel during our 
counting experiment. The following theorem of Davies [33] shows how to decompose the 
Davies process into its trajectories. 

Theorem 2.5.3: (Davies [33]) Given any ideal Davies process £^ with hounded inter- 
action rate, as defined in the previous section, we have for all t > 0, Ef, Es & T,t and 
X € M2: 



£'[Ef,Es]{X) 



EfXEs 



CO. 



In the previous section we already found the expression for the time evolution in between 
jumps: Yt. We now turn to the calculation of J'f and Jg- For all X in M2 we have: 

^(Y^ 1- lp*rro r .^w^^^Ulr^ r Jo Ad[[/,7r(z) 7r(O)({g},0)] 

Now look at UtTr{z) (g) 7r(0)({s}, 0), we find by using Theorem 2.3.1: 

C/t7r(z)®7r(O)({s},0) = V / «t (a, 0, r, 0)^^-1^1+1^1 = ztit(0, 0, 0, 0) + 

f ut(0,0,{r},0)dr + 0, 0, 0) + z / 0, {r}, 0)dr+ 

Jo Jo 

Jo Jo V^/SXpt^— 2j Z J 

Therefore we get, for all X G M2: 



■^^ ' m texp(-t|z|2) 
Kd\zl2^Vf\{X). 



z 



{X) = 



Let us now turn to the calculation of Jg- We find for all X G M2: 

Ti^c^ V ^PmiQ.\o\ ,^^S\(^c^ v Jo Ad[^,7r(^) 7r(Q)(0,{5})] {X)ds 
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Taking a closer look at Ut'K{z) (g) 7r(O)(0, {s}), applying Theorem 2.3.1: 

C/i7r(z)®7r(O)(0,{s}) = / ut(0, {s}, r, 0)^1^1 = Ui(0, {s}, 0, 0) + 

JO. 

z[ ut{$,{s},{r},^)dr + z^ f P nt(0, {s}, {n, rs}, 0)dridr2 ^ 
JO Jo Jo 

Therefore we get, for all X £ M2: 

^ ^ 40 iexp(-t|z|2) [Vk, Oyj^ ^ 

= Ad[y,](X). 



Ksexp(-|) 



Since we are driving the atom with a laser now, the time evolution when we do not observe 
the side channel nor the forward channel is now given by 

T,' = s'[n[-t,Q),n[-t,o)] 

and no longer by Tj. We will now derive the Master equation for this new time evolution. 
For this we need the Dyson series: let Lq and J be linear maps from M2 — M2, then for 
all t > 0: 

exp (f(Lo + J)) = / exp (^{ui + t)LojJ e^p (^{u!2 — u!i)LojJ . . . J exp{—uJkLo)du, 
Jni-t,o) 

where we have written lo as {ui, . . . ,LOk} with — t < < . . . < < 0. 

Now remember that {lt}t>o is a semigroup, i.e. we can write = exp(tLo). Then, using 
Theorem 2.5.3 and twice the Dyson series 

= [i^[-t, 0), n[-t, 0)] = / WY,jf,j^iu;f, Us)dujfdujs = 

Jn[-t,o)xn[-t,o) 

exp{t{Lo + Jf + Js)). 
This means we get the following Master equation: 

±T%.) = {Lo + Jf + Js)T,^i-) = 

^* (2.5) 

--{v*v, r/( . )} + - zVf, T/( • )] + y*r/( • )v, 

which is exactly the Master equation for resonance fluorescence (see [27]) if we take z = 
n 



—I 



2=y with Q, the Rabi frequency, real. 

In the quantum optics literature (see for instance [27]), usually there is no photon counting 
measurement done in the forward channel, i.e. Ef = Cl[—t, 0). Prom here on we will do the 
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same, we define for alH > and £ T,t : M^[Es] := £*[n[-t,0), Es]. In the following 
we will also suppress the index s on Es- Using the Dyson series and Theorem 2.5.3 we 
find, for alH > and E G S*: 

M'[E] = f Wz,Mdti{uj), (2.6) 
Je 

where the time evolution in between sidc-channcl-jumps Zf is given by Zt = exp (t(Lo + 
J'f)) and Wz,js is defined in the obvious way analogous to Definition 2.5.2. Note that we 
have found exactly the same jump operator and time evolution in between jumps as in 
the usual quantum optics literature, see for instance [27], [29], i.e. we have succeeded in 
constructing the Davies process of resonance fluorescence with quantum stochastic calcu- 
lus. 



2.6 A renewal process 



We will now look briefly at some features of the process we obtained. It is easily seen 
from the fact that {Jsf' = (i.e. 52(0) = 0) that the photons in the side channel arrive 
anti-hunched: the probability to see two photons immediately after each other is 0. We 
will now show that the photon counting process in the side channel is a so-called renewal 
process. 

We denote S* := E[0, t) and, via a shift, we let events E in E* correspond to events E — t 
in the output sigma field Sj. This means that an element u = {coi, . . . ,LOk} in £^ G S* 
with < wi < . . . < ojfc < t should be interpreted as seeing the first photon appear in the 
side channel at time coi, the second at time 102 up to the fc'th photon at time Uk- 



Given that wc start the photon counting measurement in the initial state p, we define 
on the sigma fields {t > 0) probability measures in the natural way: for E G : 
P*[E] := Tr(pA^*[E — t]{I)). The family of sigma fields {S*}t>o generates a sigma-field 
1;°° of n[0, 00). Using that T^{I) = I, see equation (2.5), we find for all E G S*: 

¥'+'[E] = Tt(^pM'+' [{EOn[t, t + 8))-{t + s)] (/)) = 

Tr(^pM*+' [E-{t + s)On[-s, 0)] (/)) = Tr(^pM'[E - t]M' [n[-s, 0)] (/)) = 
Tr{pM^[E - t]T^{I)) = Ty{pM*[E - t]{I)) = P* [E], 

i.e. P*+* [i?] does not depend on s. Therefore the family {Pp}t>o on the sigma-fields {Il*}t>o 
is consistent, hence by Kolmogorov's extension theorem it extends to a single probability 



measure 



on S= 
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In the following, when we write uj £ Q[0, oo) as {wi, 0^2, • • •}, we imply that < wi < a;2 < 
. . .. For j = 1, 2, . . . we define random variables: 



Xj : J7[0, 00) — ^ R"*" : a; = {a;i, 0^2, ■ ■ ■} 1— > | 



w+ . , . _ r _ 1 . , / - ^j-i if 1^1 ^ j 

00 otherwise 



where we take loq to be 0. These random variables give the time elapsed between the 
{j — l)th and jth detection of a photon. To prove that the counting process is a (modified) 
renewal process we have to show that for i = 1,2, .. . the random variables Xi are inde- 
pendent and for i = 2, 3, . . . they arc identically distributed. This means we have to show 
that for i = 2,3, ... the distribution functions Fx^ix) := Fp[Xi < x] are equal and for 
i, J = 1, 2, . . . the joint probability distribution function Fxi,Xj{x, y) ■= Fp[^i < x A Xj < 
y] factorizes: Fxi,Xjix,y) = Fxi{x)Fxj{y). 

Let us first introduce some convenient notation. Note that, using equation (2.6), we have 
for all G S*: 



¥,[E] = P* [E] =Tr(p Wz,jMM^m) = 

J E — t 

T^^(^P J^Z^iJsZu,2-ojiJs ■ ■ ■ JsZt-ojki^)diJ,{u}fj . 
We will denote: xi := loi,X2 := L02 — loi, ■ ■ ■ , Xk+i ■= t — uj^, then we can write: 

Let P denote the matrix , then we have: 

JsZ^^iP) = (l'^^l'(^-(^))22 0) = |./(Z.,(P))^^P. 

Therefore, if we define z{x) := \Ks\'^{ZxiP))^2' ^last{x) := \ks\^ {Z^il))^^ and Zfirst{x) := 
Tv(^pZx{P)) , we can write (see also [29]): 

Z first ^last 

{xk+i)dp{uj). (2.7) 

•'^ 1=2 

We would like to stress that this formula is only valid for events E e T,^ and not for all 
events in S°°. 
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For t > we introduce the following random variables: 

Nt: J^[0,oo)^N: u \u H [0,t]\, 

counting the number of photons appearing in the side channel up to time t. Since, for 

strictly positive driving field strengths, i.e. \z\ > 0, the eigenvalues of the generator Lq+Jj 
of the semigroup Zt all have strictly negative real parts, wc have lim^^oo Zt = Using 
this, formula (2.7) and the fact that the event [Nt = 0] is an clement of S*, we obtain: 

lim Fp[Nt = 0] = lim Zfirst{t) = 0. 

Now suppose we have that lim^^oo Fp[Nt < n] = for a certain n G N. For s < t we use: 
Fp[Nt <n+l]= Fp[Nt <n + l\Ns < n]¥p[Ns < n] + Fp[Nt <n + l\Ns > ra]Pp[iV, > n]. 
Therefore we have: 

lim Fp[Nt < n + 1] = lim lim Fp[Nt < n + 1] = 

t— *oo s— ♦oot— »oo 

lim lim fPp[iVt <n+ llN^ < n]Fp[Ns < n] + FplNt <n + llNg > n]Fp[Ns > n]] = 
lim hm P* [iVt <n + l\Ns > n] = lim hm ziast{t - s) = 0. 

s— >oot— »oo s— +oot— >oo 

Now using induction, wc get for n G N: 

lim FplNt < n] = 0. 

Wc arc now ready to calculate the distribution functions Fx^ and Fx^^Xj- The problem is 
that for instance the event [Xi < x] £ is not an element of S* for a t € M. We solve 
this by conditioning on the event [A^^ > i] and taking the limit for t to infinity: 

Fx,{x)=Fp[X,<x] = 

lim (Fp[Xi < x\Nt > i]Fp[Nt > i] + Fp[Xi < x\Nt < i]Fp[Nt < i]) = 



hm F^lXi <xANt>i]. 



t- 
] 

Now we use again formula (2.7) to obtain for i > 2: 

k 



Fxi{x)=limY] Zfirstixi)(Ylz{xi))ziastixk+i)dxi...dxk+i 

/*X f n k \ 

lim / z{xi) \ y~] / Zfir.st{xi)dxi(Y\z{xi)dxi)ziast{xk+i)dxk+i ]dxi 

lim [ z{xi)Fi-^' [Nt-^^ >i-l\dxi. 



1=2 



t- 

Then we use dominated convergence to interchange the limit and the integral to obtain: 



Fx,{x)= r z{x')dx'. 

JQ 
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When i = 1 we can repeat the whole calculation to find the same result when for z we 
substitute z first- It is now obvious that for i = 2,3,... the random variables Xi are 
identically distributed. 

In a similar fashion, only extracting two integrals now, we find that for i,j = 2,3,... : 
Fxi,Xj{x,y) = Jq z{x')z{y')dx'dy' . If i or j is 1 we again have to substitute Zfi^st for 
z. It is now obvious that the random variables Xj and Xj are independent. We conclude 
that the family of random variables {Xi}i=i^2,... is a (modified) renewal process. 



Chapter 3 



Stochastic Schrodinger equations 
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Toemooiveld 1, 6526 ED Nijmegen, The Netherlands 

^^EURANDOM, PO Box 513, 5600 MB Eindhoven, The Netherlands 

Abstract^ 

A derivation of Bclavkin's stochastic Schrodinger equations is given using quantum filtering theory. 
We study an open system in contact with its environment, the electromagnetic field. Continuous 
observation of the field yields information on the system: it is possible to keep track in real time 
of the best estimate of the system's quantum state given the observations made. This estimate 
satisfies a stochastic Schrodinger equation, which can be derived from the quantum stochastic 
differential equation for the interaction picture evolution of system and field together. Throughout 
the paper we focus on the basic example of resonance fluorescence. 



3.1 Introduction 

It has long been recognized that continuous time measurements can not be described by 
the standard projection postulate of quantum mechanics. In the late 60's, beginning 70's, 

^This chapter is an adapted version of [23]. 
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Davies developed a theory for continuous time measurement [31] culminating in his book 
[33]. His mathematical work became known to the quantum optics community through 
the paper with Srinivas on photon counting [82]. 

The late 80's brought renewed interest to the theory of continuous time measurement. For 
instance the waiting time distribution of fluorescence photons of a two-level atom driven by 
a laser was obtained by associating a continuous evolution to the atom in between photon 
detections and jumps at the moments a photon is detected [29]. In this way every record of 
photon detection times determines a trajectory in the state space of the atom. Averaging 
over all possible detection records leads to the well-known description of the dissipative 
evolution of the atom by a master equation. Advantage of the trajectory approach is the 
fact that an initially pure state will remain pure along the whole trajectory. This allows 
for the use of state vectors instead of density matrices, significantly speeding up computer 
simulations [73], [30], [44], [27]. 

Infinitesimally, the quantum trajectories are solutions of a stochastic differential equation 
with the measurement process as the noise term. The change in the state is given by 
the sum of two terms: a deterministic one proportional with dt and a stochastic one 
proportional to the number of detected photons dN^ in the interval dt. For other schemes 
such as homodyne detection the corresponding stochastic differential equation is obtained 
as the diffusive limit of photon counting where the jumps in the state space decrease in 
size but become increasingly frequent [9], [27], [92]. In this limit the stochastic term in 
the differential equation is replaced by a process with continuous paths. 

The stochastic Schrodinger equations obtained in this way had been postulated before by 
Gisin [47], [48], [35], [46], in an attempt to generalize the customary unitary evolution 
in quantum mechanics. The stochastic terms are seen as randomness originating from 
the measurement process. However, in this approach the correspondence between the 
different quantum state diffusion equations and the measurements that can be performed 
is not emphasized. 

Another approach originated from the development of quantum stochastic calculus [54], 
[74], generalizing the classical Ito table to quantum noises represented by creation and 
annihilation operators (see Section 3.6). Barchielli saw the relevance of this new calculus 
for quantum optics [7]. Indeed, in the Markovian approximation the interaction between 
a quantum system and the electromagnetic field is governed by a unitary solution of a 
quantum stochastic differential equation in the sense of [54] . 

Bclavkin [11] was the first to sec the connection between quantum measurement theory 
and classical filtering theory [57], in which one estimates a signal or system process when 
observing a function of the signal in the presence of noise. This is done by deriving the 
filtering equation which is a stochastic differential equation for the expectation value of the 
system process conditioned on outcomes of the observation process. Belavkin extended the 
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filtering theory [17], [16] to allow for the quantum noises of [54]. Stochastic Schrodinger 
equations turn out to be examples of the quantum filtering or Belavkin equation [15], [20]. 

Aim of this paper is to give an elementary presentation of quantum filtering theory. 
We construct the expectation of an observable conditioned on outcomes of a given mea- 
surement process. The differential form of this conditional expectation is the stochastic 
Schrodinger equation associated with the given measurement. At the heart of the deriva- 
tion lies the Ito table of quantum stochastic calculus enabling a fast computation of the 
equation. The procedure is summarized in a small recipe in Section 3.7. 

To illustrate the theory we consequently focus on the basic example of resonance fluores- 
cence of a two-level atom for which we consider photon counting and homodyne detection 
measurement schemes. The stochastic Schrodinger equations for these examples are de- 
rived in two ways, once via the usual approach using quantum trajectories and a diffusive 
limit, and once using quantum filtering theory. In this way we hope to emphasize how 
conceptually different both methods are. 

This paper is organised as follows. Sections 3.2 and 3.3 serve as an introduction to the 
guiding example of this paper: resonance fluorescence of a two-level atom driven by a 
laser. In Section 3.2 we put the photon counting description of resonance fluorescence 
by Davies [24], [29], [28] into the form of a stochastic differential equation driven by the 
counting process. In Section 3.3 we discuss the homodyne detection scheme as a diffusive 
limit of the photon counting measurement, arriving at a stochastic differential equation 
driven by a diffusion process. The equations of Sections 3.2 and 3.3 will be rederived later 
in a more general way using quantum filtering theory. 

In Section 3.4 we introduce the concept of conditional expectation in quantum mechanics 
by first illustrating it in some simple, motivating examples. Section 3.5 describes the 
dissipative evolution of the open system within the Markov approximation. The joint 
evolution of the system and its environment, the quantized electromagnetic field, is given 
by unitaries satisfying a quantum stochastic differential equation. Given a measurement 
of some field observables it is shown how to condition the state of the system on outcomes 
of the measurement using the construction of Section 3.4. Section 3.6 is a short review 
of quantum stochastic calculus and its applications to open systems. Sections 3.5 and 3.6 
describe dilation theory and quantum stochastic calculus in a nutshell. 

Section 3.7 contains the derivation of the quantum filtering equation, the stochastic differ- 
ential equation for the conditional expectation. This equation is the stochastic Schrodinger 
equation for the given measurement. This part ends with a recipe for computing stochas- 
tic Schrodinger equations for a large class of quantum systems and measurements. The 
end of the article connects to Sections 3.2 and 3.3 by showing how the recipe works in our 
main example. 



40 



CHAPTER 3. STOCHASTIC SCHRODINGER EQUATIONS 



3.2 The Davies process 



We consider a two-level atom in interaction with the quantized electromagnetic field. The 
state of the atom is described by a 2 x 2-density matrix p, i.e. p > 0, and Trp = 1. Atom 
and field together perform a unitary, thus reversible evolution, but by taking a partial 
trace over the electromagnetic field we are left with an irreversible, dissipative evolution 
of the atom alone. In the so called Markov limit it is given by a norm continuous semigroup 
{Tt}t>o of completely positive maps. A central example discussed in this paper is resonance 
fluorescence. Here the atom is driven by a laser on the forward channel, while in the side 
channel a photon counting measurement is performed. For the time being we will suppress 
the oscillations of the laser for reasons of simplicity. In this case the Lindblad generator 
of Tt, or Liouvillian L is given by (cf. [27]): 



Ttip) = L{p) = -i[H,p]+i^[V + V*,p] - l{V*V,p} + VpV\ (3.1) 



1 



d 

di 

where 

V 



H := ^(Tz is the Hamiltonian of the atom, and J7 is the Rabi frequency. 
The master equation (3.1) can be unravelled in many ways depending on what photon 
detection measurement is performed. By unravelling the master equation we mean writing 
L as the sum C + J', where J represents the instantaneous state change taking place when 
detecting a photon, and C describes the smooth state variation in between these instants. 
The unravelling for photon counting in the side channel is given by [27] 

C{p) = -i[H,p]+i^[V + V*,p] - ^{V*V,p} + (1 - \Ks\^)VpV* and J{p) = |«.|VpF*, 
with the decay rate into the side channel. 

An outcome of the measurement over an arbitrary finite time interval [0, t) is the set of 
times {ti,t2, ■ ■ ■ ,tk} at which photons are detected in the side channel of the field. The 
number of detected photons can be arbitrary, thus the space of outcomes is 

oo oo 

^ ([0, t)) := U nn ([0, t)) =\J{aC [0, t); \a\ = n} 

n=0 n=0 

also called the Guichardet space [52]. In order to describe the probability distribution of 
the outcomes we need to make Q,{[0,t)) into a measure space. Let us consider the space 
of n-tuples [0, t)"' with its Borel cr-algebra and the measure where A„ is the Lebesgue 
measure. Then the map 

jn:[0,tT3{ti,...,tn)^{tl,...,tn}enn{[0,t)) 

induces the cr-algebra S„([0,t)) and the measure /i„ on f2„([0,t)). We define now the 
measure /i on O {[0,t)) such that /i({0}) = 1 and fi = fin on r2„ {[0,t)). We will abbreviate 
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O.{[0,t)) and S([0,t)) to and S*, respectively. 

Davies was the first to show [33] (see also [27] , [24] ) that the unnormalized state of the 2- 
level atom at time t with initial state p, and conditioned on the outcome of the experiment 
being in a set E e E* is given by: 

M'[E]{p) = [ Wt{u){p)dfi{u), 
Je 

where for a; = {ii , . . . , tfe } G $7* with < ti < . . . < tk < t we have 

Wt{u){p) := exp {{t - tk)C)j ... J exp ((ts - h)C)jeicp {tiC){p). 

Furthermore, P* [£;] := Tt{M*[E\{p)) is the probability that the event E occurs if the initial 
state is p. The family of prabability measures {P*}t>o is consistent, i.e. P^+^[£;] = P* [£■] 
for all -E G S*, s > 0, see [24], hence by Kolmogorov's extension theorem it extends to a 
single probability measure on the cj-algebra of the set 
On the measure space (0°°, Pp) we define the following random variables: 

Nt: 0°°^N: w ^ n [0, 

counting the number of photons detected in the side channel up to time t. The counting 
process {A^t}t>o has differential dNt := A^t+^t — Nt satisfying dNt{io) = 1 if t G w and 
dNt{uj) = otherwise. Therefore we have the following Ito rules: dNfdNt = dNt and 
dNtdt = 0, (cf. [9]). 

To emphasise the fact that the evolution of the 2-levcl atom is stochastic, we will regard 
the normalized density matrix as a random variable {pi}t>o with values in the 2x 2-density 
matrices defined as follows: 

f f Wt(u; n\0, t))(p) , , 

p* : n~ ^ M2 : a; ^ p^ := —j^ V ^ ^3 2) 



Tv(m('^n[o,t))(p))' 



The processes Nt and p^ are related through the stochastic differential equation dp^ = 

atdt + f3tdNt. Following [9] wc will now determine the processes at and Pt by differentiating 
(3.2). If t G a; then dNt{uj) = 1, i.e. the differential dt is negligible compared to dNt = 1, 
therefore: 

(3.3) 

On the other hand, li t ^ u then dNt(uj) = 0, i.e. dNt is negligible compared to dt. 
Therefore it is only the dt term that contributes: 

d exp{{s-t)ji:){pl) 
at(w) = — 



ds -tTv(exp((s-t)£)(p*,)) 



(3.4) 
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where we used that Tr(£(p^)) = —Tr(i7(p^)) , as a consequence of the fact that Tr(L(cj)) = 
for all density matrices a. Substituting (3.3) and (3.4) into dp* = atdt + (5tdNt we get 
the following stochastic Schrddinger equation for the state evolution of the 2-level atom if 
we are counting photons in the side channel (cf. [9], [28], [19]): 

dpi = L{pi)dt + {t^^^ - pi) {dNt - Tr{j{pi))dt) . (3.5) 

The differential dMf := dNt — Tr(j'(p*))dt and the initial condition Mq = define an 
important process Mf called the innovating martingale, discussed in more detail in Section 
3.7. 



3.3 Homodyne detection 



We change the experimental setup described in the previous section by introducing a local 
oscillator, i.e. a one mode oscillator in a coherent state given by the normalised vector in 

/2(N) 

i;{at) := exp (^^) (1, . . .), (3.6) 

for a certain at G C. Wc take at = where wt is a complex number with modulus 
\wt\ = 1, and e > 0. The number e is inversely proportional to the intensity of the 
oscillator. Later on we will let the intensity go to infinity, i.e. e ^ 0. The phase (pt of the 
oscillator is represented by wt = exp(z^(), with (pt = (f)o + uiot, where ujio is the frequency 
of the oscillator. 

The local oscillator is coupled to a channel in the electromagnetic field, the local oscillator 
beam. The field is initially in the vacuum state. The local oscillator and the field are 
coupled in such a way that every time a photon is detected in the beam, a jump on the 
local oscillator occurs, given by the operation 

Jioip) = AiopAl, (3.7) 

where Ai^ is the annihilation operator corresponding to the mode of the local oscillator. 
The coherent state V'(cKt) is an eigenstate of the jump operator Ai^ at eigenvalue at- 
Now we are ready to discuss the homodyne detection scheme. Instead of directly counting 
photons in the side channel we first mix the side channel with the local oscillator beam 
with the help of a fifty-fifty beam splitter. In one of the emerging beams a photon counting 
measurement is performed. A detected photon can come from the atom through the side 
channel or from the local oscillator via the local oscillator beam. Therefore the jump 
operator on states a of the atom and the oscillator together, is the sum of the respective 
jump operators: 



Ja&o{(T) = {ksV ® / + / (g) Aio)a{KsV* (8) / + / ® Al). 
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An initial product state pfS> \'ip{at)){ip{ctt)\ of the 2-levcl atom and the local oscillator will 
remain a product after the jump since tp{at) is an eigenvector of the annihilation operator. 
Tracing out the local oscillator yields the following jump operation for the atom in the 
homodyne setup: 

Ja{p) = Ticio[ja^io{p<^ \^{at)){i;{at)\)) = {ksV + ^)p{ksV* + ^). 

In the same way as in Section 3.2, we can derive the following stochastic Schrodinger 
equation for the state evolution of the two-level atom when counting photons after mixing 
the side channel and the local oscillator beam [9] [28] : 

dpi = L{pi)dt + - pi)e(dNt - Tr{Mpi))dt), (3.8) 

where the extra e's are introduced for future convenience. We will again use the abbrevi- 
ation: dM^ = dNt — Tr [Jaipi)) dt for the innovating martingale (see Section 3.7). In the 
homodyne detection scheme the intensity of the local oscillator beam is taken extremely 
large, i.e. we are interested in the limit e —> [9], [27], [92]. Then the number of detected 
photons becomes very large and it makes sense to scale and center Nt, obtaining in this 
way the process with differential dW^ := edNt — dt/e and Wq = 0. We find the following 
Ito rules for dW^: 

dWfdWf = {edNt - ^dt) {edNt - ^dt) = e^dNt = edW^ + dt, 
dWidt = 0. 

In the limit e ^ this becomes dWtdWt = dt and dWtdt = 0, i.e. the process Wt := 
lim£_>o Wt is a diffusion. It is actually this scaled and centered process that is being 
observed and not the individual photon counts Nt, see [27]. We pass now to the evaluation 
of the limit of (3.8): 

hm - ( f;.^^'! - pi) = wtKspiV* + wtKsVpi - Tr{wtKspiV* + wtKsVpi)pi. 

This leads to the following stochastic Schrodinger equation for the homodyne detection 
scheme [9], [28], [92], [19] 

Api = L{pi)dt + {wtKspiV* + WtKsVpi - TT{wt-KspiV* + wtKsVpi)pi ) dMt^'^, (3.9) 

for all states p G M2, where 

dM^"^ := dWt - Ti{wt-KspiV* + wtKsVpi)dt. (3.10) 



Let as{t) and ab{t) denote the annihilation operators for the side channel and the local 
oscillator beam, respectively. They satisfy the canonical commutation relations 



[ai{t),a*{r)] = SijS{t-r), i,j G {s,b}. 
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Smearing with a quadratically integrable function / gives 

Mf) = J f{t)ai{t)dt, ie{s,b}. 

By definition, the stochastic process {Nt}t>o counting the number of detected photons 
has the same law as the the number operator A(t) up to time t for the beam on which the 
measurement is performed. Formally we can write 



Ht)= / {a*sir)(g)I + I(g>al{r)){asir)(g>I + I(g)abir))dr. 
Jo 

The oscillator beam is at time t in the coherent state ('^)' where ft G L^(]R) is the 
function r i— WrXio,t]{f)- Since the state of the local oscillator beam is an eigenvector of 
the annihilation operator ab{r) 

<r)^ (f ) = (f ) , 

we find 

sAit) -i = sAs{t) ®I + e f C^alir) + ^a,(r)) ®I+ - - 

E Jq EE EE 

= eAs{t)®I+ {Al{ft)+As{ft))®I. 

The operator := A* [ft) + As{ft) is called a field quadrature. We conclude that in 

the limit e — > the homodyne detection is a setup for continuous time measurement of 
the field quadratures X^{t) of the side channel, (cf. [27]). 



3.4 Conditional expectations 



In the remainder of this article we will derive the equations (3.5) and (3.9) in a different 
way. We will develop a general way to derive Bclavkin equations (or stochastic Schrodinger 
equations). The counting experiment and the homodyne detection experiment, described 
in the previous sections, serve as examples in this general framework. The method we 
describe here closely follows Belavkin's original paper on quantum filtering theory [17]. 
The construction below, however, uses explicitly the decomposition of operators over the 
measurement results. In the next section it will turn out that this is done most naturally 
in the interaction picture. 

Let us remind the concept of conditional expectation from probability theory. Let (fi, S, P) 
be a probability space describing the "world" and S' C S a cr-algebra of events to which 
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"we have access". A random variable / on (17, S,P) with IE(|/|) < oo can be projected to 
its conditional expectation £{f) which is measurable with respect to S' and satisfies 

[ fdP= [ £if)dF 
Je Je 

for all events E in S'. Our information about the state of that part of the world to 
which we have access, can be summarized in a probability distribution Q on S'. Then the 
predicted expectation of / given this information is £{f)dQ. We will extend this now 
to quantum systems and measurements. 

The guiding example is that of an n level atom described by the algebra B := Mn un- 
dergoing a transformation given by a completely positive unit preserving map T : B ^ B 

with the following Kraus decomposition T{X) = ^jg^ V*XVi. The elements of Q can be 
seen as the possible measurement outcomes. For any initial state p of B and measurement 
result i e Q, the state after the measurement is given by 

= V,pV;/Ti:{VipV*), 

and the probability distribution of the outcomes is p = "^^^q Pi^i where 6i is the atomic 
measure at i, and Pi = TYiVipV^ ), which without loss of generality can be assumed to be 
strictly positive. We represent the measurement by an instrument, that is the completely 
positive map with the following action on states 

Let X G ;B be an observable of the system. Its expectation after the measurement, given 
that the result i G has been obtained is Tr(piX). The function 

£{X) : n^C: Tr(piX) 

is the conditional expectation of X onto £°°{^). If q = ^qidi is a probability distribution 
on Vl then ^ qi£{X){i) represents the expectation of X on a statistical ensemble for which 
the distribution of the measurement outcomes is q. We extend the conditional expectation 
to the linear map 

£: B^e°°{n) ^ e^^in) cB^e'^in) 

such that for any element A : i Ai in B i^i {^) = i°° ^ B) we have 

£{A) : i ^ Ti{piAi). 

This map has the following obvious properties: it is idempotent and has norm one. More- 
over, it is the unique linear map with these properties preserving the state M.{p) on 
B ® For this reason we will call £, the conditional expectation with respect to 

M{p). Its dual can be seen as an extension of probability distribitions q G to states 

onB(^i'^{n) 

£* ■ q^'^Pi® qA- 
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Thus while the measurement (3.11) provides a state M{p) on B <^i°^{^}), the conditional 
expectation with respect to M.{p) extends probability distributions q G £^{Q.) of outcomes, 
to states on B (g) £°°{n), and in particular on B which represents the state after the 
measurement given the outcomes distribution q. 

With this example in mind we pass to a more general setup which will be needed in 
deriving the stochastic Schrodinger equations. Let ^ be a unital *-algebra of bounded 
operators on a Hilbert space H whose selfadjoint elements represent the observables of a 
quantum system. It is natural from the physical point of view to assume that A is strongly 
closed, i.e. if {^n}n>o is a sequence of operators in A such that H^nV'll ~^ II^V'll ^oi any 
vector in H and a fixed bounded operator A, then A e A. Prom the mathematical 
point of view this leads to the rich theory of von Neumann algebras inspired initially by 
quantum mechanics, but can as well be seen as the generalization of probability theory 
to the non-commutative world of quantum mechanics. Indeed, the building blocks of 
quantum systems are matrix algebras, while probability spaces can be encoded into their 
commutative algebra of bounded random variables L°°(J7, S, P) which appeared already in 
the example above. A state is described by a density matrix in the first case or a probability 
distribution in the second, in general it is a positive normalized linear functional ip : A ^ C 
which is continuous with respect to the weak*-topology, the natural topology on a von 
Neumann algebra seen as the dual of a Banach space [56] . 

Definition 3.4.1: Let ;8 be a von Neumann subalgcbra of a von Neumann algebra A of 
operators on a (separable) Hilbert space H. A conditional expectation of A onto ;S is a 
linear surjective map 6 : A^ B, such that: 

1. £^ = £ {£ is idempotent), 

2- "^AeA '■ 11^(^)11 ^ ll^ll (^ is normcontr active) . 

In [85] it has been shown that the conditions 1 and 2 are equivalent to £ being an identity 
preserving, completely positive map, and satisfying the module property 

£{BiAB2) = Bi£{A)B2, for all Bi,B2 G B, and AeA, (3.12) 

generalizing a similar property of conditional expectations in classical probability theory 
(cf. [89]). 

In analogy to the classical case we are particularly interested in the conditional expecta- 
tion which leaves a given state p on ^ invariant, i.e. po £ = p. However such a map does 
not always exist, but if it exists then it is unique [84] and will be denoted £p. Using £p we 
can extend states cr on H to states ao£p of A which should be interpreted as the updated 
state of A after receiving the information (for instance through a measurement) that the 
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subsystem B is in the state a (of. [62]). 

In the remainder of this section we will construct the conditional expectation £p from a 
von Neumann algebra A onto its center C := {C e A; AC = CA for all A e A} leaving 
a given state p on ^ invariant. The center C is a commutative von Neumann algebra 
and is therefore isomorphic to some L°°{Q,Ti,¥). In our guiding example the center of 
B (8) is £°°(r2). Later on (see section 3.6) this role will be played by the commu- 

tative algebra of the observed process with CI the space of all paths of measurement records. 



Theorem 3.4.2: There exists a unique conditional expectation £p : A^ C which leaves 
the state p on A invariant. 

Proof. The proof is based on the central decomposition of A [56]. In our guiding example, 
B (g) is isomorphic to (Bi^nBi where the BiS are copies of B. In general we can 

identify the center C with some L°°(ri, S,P) where P corresponds to the restriction of p 
to C. We will ignore for simplicity all issues related with measurability in the following 
constructions. The Hilbert space H has a direct integral representation H = M.i_jF{dLo) 
in the sense that there exists a family of Hilbert spaces {Mj^jt^gn and for any ijj eM there 
exists a map u ip^ e such that 

(•0,0) = / {^|J^,cP^)F{du;). 
Jn 

The von Neumann algebra A has a central decomposition A = AuiF{duj) in the sense 
that there exists a family {Auj}iueft of von Neumann algebras with trivial center, or factors, 
and for any A E A there is a map co i— > Ai_^ E Au} such that {Aip)i_j = A^^^^^ for all ^ G EI 
and P-almost all a; G The state p on A has a decomposition in states p^ on Au, 
such that for any ^ G .4 its expectation is obtained by integrating with respect to P the 
expectations of its components A^,: 

p{A) = I p^(A^)P(da;). (3.13) 

The map Ep-. A —>■ C defined by 

£p{A) : u ^ Puj{Auj) 

for all A^ A is the desired conditional expectation. One can easily verify that this map is 
linear, identity preserving, completely positive (as a positive map onto a commutative von 
Neumann algebra), and has the module property. Thus, £p is a conditional expectation 
and leaves the state p invariant by 3.13. Uniqueness follows from [84]. □ 

It is helpful to think of the state p and an arbitrary operator A as maps p, : u ^ 
Pu,, and respectively A, : u ^ A^,. The conditional expectation £p{A) is the function 
p,{A,) : u ^ Puj{Au,). 
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3.5 The dilation 

Let B be the observable algebra of a given quantum system on the Hilbert space H. In the 
case of resonance fluorescence B will be all 2 x 2 matrices M2, the algebra of observables 
for the 2-level atom. The irreversible evolution of the system in the Heisenberg picture is 
given by the norm continuous semigroup {Tt}t>o of completely positive maps Tt : B ^ B. 
By Lindblad's theorem [67] we have Tt = exp(fL) where the generator L : B ^ B has the 
following action 

k 

L{X)=i[H,X]+J2v*XV^ - -{V;Vj,X}, (3.14) 

i=i 

where H and the Vj 's are fixed elements of B, H being selfadjoint. 

We can see the irreversible evolution as stemming from a reversible evolution of the system 
B coupled to an environment, which will be the electromagnetic field. We model a channel 
in the field by the bosonic or symmetric Fock space over the Hilbert space L?{M) of square 
integrable wave functions on the real line, i.e. 

00 

J^:=C©0L2(R)®»". 

n=l 

The algebra generated by the field observables on F contains all bounded operators and 
we denote it by W. For the dilation we will need k independent copies of this algebra 

The free evolution of the field is given by the unitary group St^ the second quantization 
of the left shift s{t) on L^(M) , i.e. s{t) '■ f ^ f{- + t). In the Heisenberg picture the 
evolution on W is 

W ^ S*tWSt := Ad[St]{W). 

The atom and field together form a closed quantum system, thus their joint evolution is 
given by a one-parameter group {TijtgK of *-automorphisms on i3 W®'': 

X ^ UtXUt := Ad[Ut]{X). 

The group Ut is a perturbation of the free evolution without interaction. We describe this 
perturbation by the family of unitaries Ut := S®tUt for all t G R satisfying the cocycle 
identity 

Ut+s = S^^UtSf'Us, for all t,seR. 

The direct connection between the reduced evolution of the atom given by (3.14) and the 
cocycle Ut is one of the important results of quantum stochastic calculus [54] which makes 
the object of Section 3.6. For the moment we only mention that in the Markov limit, Ut 
is the solution of the stochastic differential equation [54] , [74] , [72] 



dUt = {VjdA* (t) - V*dAj (t) - {iH + ^V*Vj)dt}Ut, Uo = 1, (3.15) 
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where the repeated index j is meant to be sTimmed over. The quantum Markov dilation 
can be summarized by the following diagram (see [59], [60]): 

^ B 

'id®^®'^ (3.16) 

i.e. for all X G ^ : Tt{X) = (id (g) c/)'^^) {ft{X ® 1®^)), where ^ is the vacuum state on 
W, and 1 is the identity operator in W. Any dilation of the semigroup Tt with Bosc fields 
is unitarily equivalent with the above one under certain minimality requirements. The 
diagram can also be read in the Schrodinger picture if we reverse the arrows: start with 
a state p of the system B in the upper right hand corner, then this state undergoes the 
following sequence of maps: 

p®(t>®^^{p® (j>®^) oft = ft,{p ® <^®*^) ^ T¥^«fe {ft,{p ® (1)^'')) . 

This means that at t = 0, the atom in state p is coupled to the k channels in the vacuum 
state, and after t seconds of unitary evolution we take the partial trace taken over the k 
channels. 

We would now like to introduce the measurement process. It turns out that this can be 
best described in the interaction picture, where we let the shift part of Ut = Sf^Ut act on 
the observables while the cocycle part acts on the states: 

p\X):=p®ct>'^\u;XUt) (3.17) 

for all X G H (g) W®^. It is well known that for the Bose field for arbitrary time t we can 
split the noise algebra as a tensor product 

with each term being the algebra generated by those fields over test functions with support 
in the corresponding subspace of L^(M): 

l2(R) = l2((-oo,0)) ®L^{[Q,t))®L^{[t,oo)). 

Such a continuous tensor product structure is called a filtration and it is essential in 
the development of quantum stochastic calculus reviewed in Section 3.6. The observables 
which wc measure in an arbitrary time interval [0, t) form a commuting family of selfadjoint 
operators {Ys}o<s<t whose spectral projections belong to the middle part of the tensor 
product VV[o,t). In the Davies process Yg = A(s), i.e. the number operator up to time 
s, while in the homodyne case Yg = X^{s). Notice that the part Wq) will not play any 
significant role as it corresponds to "what happened before we started our experiment" . 



B 



Id(8ii® 
25 0>V®fc 
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Let Ct be the commutative von Neumann generated by the observed process up to time 
t-i {^s}o<s<t (t > 0), seen as a subalgebra of B ® W®'^. By a theorem on von Neumann 
algebras, Ct is equal to the double commutant of the observed process up to time t: 
Ct = {Y^; < s < t}", with the commutant S' of a subset 5 of i3 (g) W^'^ being defined by 
S' := {X eB(S) W'^''; XS = SXyS E S}. The algebras {Ct}t>o form a growing family, 
that is Cg C Ct for all s < t. Thus we can define the inductive limit Coo '■= limt^ooCt, 
which is the smallest von Neumann algebra containing all Ct- On the other hand for each 
t > wc have a state on Ct given by the restriction of the state p* of the whole system 
defined by (3.17). We will show now that the states p* for different times "agree with each 
other" . 

Theorem 3.5.1: On the commutative algebra C^o there exists a unique state p°° which 

coincides with when restricted to Ct C Coo, for all t > 0. In particular there exists 
a measure space (r2,S,Pp) such that {Coo,p°°) is isomorphic with L°°(r2,S,Pp) and a 
growing family {St}t>o of a-subalgebras of S such that {Ct, p*) = L°°(0, St, Pp). 



Proof. In the following we will drop the extensive notation of tensoring identity operators 
when representing operators in W[s,t) for s,t eM. Let X G Cs, in particular X G 

By (3.15), Ut & B® VVj®^^, because the coefficients of the stochastic differential equation 

lie in H ® >Vg^). This implies that S^^UtSf" G H Wg^^,). Using the tensor product 

structure of W^'^, we see that Wjq^^ and B (g) W[®^_|.^) commute, and in particular X 

commutes with S^^UtSf". Then 

p*+^(x) = p'{Ut%xut+s) = p'\u:isf%s"frxs"^%sfUs) 

= p\u:XUs)=p''{X). (3.18) 

This implies that the limit state p°° on Coo with the desired properties exists, in analogy 
to the Kolmogorov extension theorem for probability measures. As seen in the previous 
section, (Coo, P°°) is isomorphic to L°°(p,, S, Pp) for some probability space (O, S, Pp). The 
subalgebras {Ct, p^) are isomorphic to L°°(J7, Ej, Pp) for some growing family {T,t}t>Q of 
(T-subalgebras of E. □ 



Remark. From spectral theory it follows that the measure space (fi*, E^) coincides with 
the joint spectrum of {Ys}s<t, i-e. f^* is the set of all paths of the process up to time t. 
For the example of the counting process this means that is the Guichardet space of the 
interval [0,t), which is the set of all sets of instants representing a "click" of the photon 
counter, i.e. it is the set of all paths of the counting process. 

We define now At '■= C't for alH > , i.e. At is the commutant of Ct, then Ct is the center 
of the von Neumann algebra At- Notice that the observable algebra of the atom B is 
contained in At- By Theorem 3.4.2 we can construct a family of conditional expectations 
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{£^^t ■ -^t Ct}t>o- For each t, Sj^t depends on the state of the "world" at that moment 
p*, keeping this in mind we wih simply denote it by An important property of £^ is 
that p°° o f * = p* o f * = p*, since the range of £* is Ct and £* leaves p* invariant. 

For an element X G At, is an clement in Ct, i.e. a function on i}f Its value in a 

point Lu £ Qt, i-e. an outcomes record up to time t, is the expectation value of X given the 
observed path to after t time units. We will use the notation £'^{X) := p* (X,) defined in 
the end of Section 3.4 to emphasise the fact that this is a function on Clf. When restricted 
to H (g) Ct the conditional expectation is precisely of the type discussed in our guiding 
example in Section 3.4. 

There exists no conditional expectation from B W onto Ct since performing the mea- 
surement has demolished the information about observables that do not commute with 
the observed process [17]. We call At the algebra of observables that are not demolished 
[17] by observing the process {Ys}o<s<t- This means that performing the experiment and 
ignoring the outcomes gives the same time evolution on At as when no measurement was 
done. 

From classical probability it follows that for all t > there exists a unique conditional 
expectation : Coo ~^ Ct that leaves the state p°° invariant, i.e. p°° o = p°°. These 
conditional expectations have the tower property, i.e. E* o E* = E^ for alH > s > 0, which 
is often very useful in calculations. E° is the expectation with respect to P^, and will 
simply be denoted Ep. Note that the tower property for s = is exactly the invariance of 
the state p°°(= Ep). 



3.6 Quantum stochastic calculus 



In this section we briefly discuss the quantum stochastic calculus developed by Hudson 
and Parthasarathy [54]. For a detailed treatment of the subject we refer to [74] and 
[72]. Let T{Ti.) denote the symmetric (or bosonic) Fock space over the one particle space 
n := C'' (g> L2(R+) = L2({1, 2, ... A;} x M+). The space describes the k channels we 
identified in the electromagnetic field. As in the previous section we denote the algebra 
of bounded operators on the one channel Fock space by W, and on the k channels 

T{n) by W®^ 

For every f eH we define the exponential vector e(/) € J^iji.) in the following way: 

oo ^ 

which differs from the coherent vector by a normalization factor. The inner products of 
two exponential vectors e(/) and e{g) is {e{f),e{g)) = exp((/, 5)). Note that the span 
of all exponential vectors, denoted V, forms a dense subspace of J^{H). Let fj be the 
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j'th component of f £ 7i for j = 1,2, ... ,k. The annihilation operator Aj{t), creation 
operator A*{t) and number operator Ajj(i) are defined on the domain V by 

Mt)e{f) = iXm, fMf) = ffj{s)ds e{f) 

Jo 

{e{g),A*{t)e{f)) = {gj,X[o,t]){e{9), e(/)) = / gj{s)ds eMif, d)) 

Jo 

{e{g), Aij{t)e{f)) = {gi, X[o,t]fj){e{g), e(/)) = / gi{s)fj{s)ds ex-p{{f,g)). 

Jo 

The operator Ajj(t) is the usual counting operator for the i'th channel. Let us write 
L^(R+) as direct sum L^([0, t]) © oo]), then J^(L^(R+)) is unitarily equivalent with 

J^{L^{[0,t])^J^{L^[t,oo)) through the identification e(/) ^ e(/t])®e(/[j), with /^j = fx[o,t] 
and /[t = fx[t,oo)- We will also use the notation /[^ for fx[s,t] ^^'^ omit the tensor product 
signs between exponential vectors. The same procedure can be carried out for all the k 
channels. 

Let Mf be one of the processes Aj(t),Aj(t) or Ajj(t). The following factorisation property 
[54], [74] makes the definition of stochastic integration against Mt possible 

{Mt - Ms)e{f) = e(/,]){(Mt - M,)e(/[,,t])}e(/[0, 

with (Mt - Ms)e{fis,t]) G •^(C' L'^{[s, t])) . We first define the stochastic integral for the 
so called simple operator processes with values in the atom and noise algebra B ® H"^*^ 
where B := M„. 

Definition 3.6.1: Let {Ls}o<s<t be an adapted (i.e. Ls € B yVg] for < s < t) 
simple process with respect to the partition {sq = 0, si, . . . , Sp = t} in the sense that 
Lg = Lg. whenever sj < s < Sj+i. Then the stochastic integral of L with respect to M on 
C" (g) P is given by [54], [74]: 

/ LsdMs fe{u) := ^ (L,,/e(u,.])) ((M,,^, - M,,)e(n[,.,,.^,]))e(^X[, J. 

^0 



By the usual approximation by simple processes wc can extend the definition of the 
stochastic integral to a large class of stochastically intcgrablc processes [54], [74]. We 
simplify our notation by writing dXt = LtdMt for Xt = Xq + Jq LgdMg. Note that the 
definition of the stochastic integral implies that the increments dMg lie in the future, i.e. 
dMg G y^[s- Another consequence of the definition of the stochastic integral is that its 
expectation with respect to the vacuum state cp is always due to the fact that the incre- 
ments dAj , dA* , dAij have zero expectation values in the vacuum. This will often simplify 
calculations of expectations, our strategy being that of trying to bring these increments 
to act on the vacuum state thus eliminating a large number of differentials. 
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The following theorem of Hudson and Parthasarathy extends the Ito rule of classical prob- 
ability theory. 

Theorem 3.6.2: (Quantum ltd rule [54], [74]) Let Mi and M2 each be one of the 
processes Aj,A* or Ajj. Then M1M2 is an adapted process satisfying the relation: 

d{MiM2) = MidM2 + M2dMi + dMidM2, 

where dMidM2 is given by the quantum ltd table: 



dMi\dM2 


dA* 


dAij 


dAi 


dAl 











dAki 


5udAl 


^lidAkj 





dAk 


Skidt 


SkidAj 






Notation. The quantum Ito rule will be used for calculating differentials of products of 
Ito integrals. Let {^i}i=i,...,p be Ito integrals, then 

d{ZiZ2...Zp)= Yl M 

uC{l,...,p} 

where the sum runs over all non-empty subsets of {1, ... ,p} and for any ly = {ii, . . . i^}, the 
term [i'] is the contribution to d{ZiZ2 ■ ■ ■ Zp) coming from differentiating only the terms 
with indices in the set {ii, . . . ik} and preserving the order of the factors in the product. 
For example the differential d{ZiZ2Z^) contains terms of the type [2] = Zi{dZ2)Z^, [13] = 
{dZi)Z2{dZ^), and [123] = {dZi){dZ2){dZ^). 

Let Vj for j = 1,2, ... ,k, and H be operators in B with H is selfadjoint. Let S" be a 
unitary operator on C" /^({l, 2, . . . , k}) with Sij = {i, Sj) G B the "matrix elements" in 
the basis {\i >: i = l,...,k} of C'^. Then there exists a unique unitary solution for the 
following quantum stochastic differential equation [54], [74] 

dUt = {VjdA*{t) + {Sij - 5ij)dAij{t) - V*SijdAj{t) - {iH + ^V*Vj)dt}Ut, (3.19) 

with initial condition Uq = 1, where again repeated indices have been summed. Equation 
(3.15), providing the cocycle of unitaries perturbing the free evolution of the electromag- 
netic field is an example of such an equation. The terms dAjj in equation (3.19) describe 
direct scattering between the channels in the electromagnetic field [10]. We have omitted 
this effect for the sake of simplicity, i.e. we always take Sij = Sij. 

We can now check the claim made in Section 3.5 that the dilation diagram 3.16 commutes. 
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It is easy to see that following the lower part of the diagram defines a semigroup on B. 
We have to show it is generated by L. For all X G B we have 

d Id (g) /(Ti(X (g) l®'^)) = Id (g) U^X (g l^f^Ut). 

Using the Ito rules we obtain 

d U;x (g l®^Ut = {dU;)X (g l®^Ut + U;x ® l^'^dUt + idU;)X (g l^'^dUt. 

With the aid of the Ito table we can evaluate these terms. We are only interested in the 
dt-terms since the expectation with respect to the vacuum kills the other terms. Then we 
obtain: d Id ® (l)''{U^X (g l^'^C/j) = Id (g (j)^{U^L{X) ® l^''Ut)dt, proving the claim. 



Now we return to the example of resonance fluorescence. Suppose the laser is off, then we 
have spontaneous decay of the 2-level atom into the field which is in the vacuum state. 
For future convenience we already distinguish a forward and a side channel in the field, 
the Liouvillian is then given by 

L{x) = i[H,x] + ^ v;xv, - ^{v;v,,x}, 

where 

^=(10)' ^f = '^f^^ V, = K,V, \Kf\' + \Ks\^ = l, 

with and |ksP the decay rates into the forward and side channel respectively. 

The dilation of the quantum dynamical system {M2,{Tt = cxp{tL)}t>o), is now given by 
the closed system (M2 ® W/ (g) Ws, {Tt}teM.) with unitary cocycle given by 

dUi"^ = {VfdA}it) - V^dAfit) + VsdA*{t) - V;dA,{t) - (iH + ^V*V)dt}Ut"^, U^'^ = 1, 

where the superscript sd reminds us of the fact that the laser is off, i.e. we are considering 
spontaneous decay. We can summarize this in the following dilation diagram 



^ Tt=exp(tL)^ ^ 



Id(g)i(g)i 



H (g W/ (g — H (g W/ (g 



where U^'^ is given by St ® StUi'^ for t > 0. 



We change this setting by introducing a laser on the forward channel, i.e. the forward 
channel is now in a coherent state (see 3.6) 7/1 := {il^{h), -ipih)) for some h G L^(M+). 
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This leads to the following dilation diagram 



'T'/l 

B — ^ B 



Id(g)i(g)i 



Id(S)7/j(g)0 (3.20) 



B^Wf^Ws — B^Wf^Ws 

i.e. the evolution on B has changed and it is in general not a semigroup. Denote by W{h) 
the unitary Weyl or displacement operator defined on V by: W{h)ilj{f) = exp(— 2zlm(/i, /)) 
tpif + h). Note that W{h)(l) = W{h)tp{0) = tpih), so that we can write 

Tj'iX) = Id 0-fhCS> Huf*x 1 ® IC^t"^) = 
Id ® ® (l){Wf{h)*Uf'^*X 1® lUf'^Wfih)) = 
Id ® ^{Wfiht]TUt"^*X ® 1 ® lJ7/'^W/(/it])), 

where /i^] := hx[o,t] and := 1 (g) W{h) (g) 1. Defining [/"t := C//'^W/(/it]), together 

with the stochastic differential equation for Wf{h^) [74] 

dWfiht]) = {h{t)dA}{t) - h{t)dAf{t) - ^\h{t)\^dt}Wf{ht]), Wf{ho) = 1, 
and the Ito rules leads to the following quantum stochastic differential equation for Uf. 

dUt= {{Vf + h{t))dA*{t)-{V* + h{t))dAf{t) +VsdA*{t)-V;dA,{t) - 
-{iH + l{\h{t)\'^ + V*V + 2h{t)Vf*))dt }Ut, Uo = l. 

Define Vf := Vf + h{t), Vs := K and H := H + i\(h{t)Vf - h{t)VJ) then this reads 

dUt = Yl i^adAUt) - KdA^ - {iH + ^v;V^)dt}Ut, Uo = 1. (3.21) 

The time dependent generator of the dissipative evolution in the presence of the laser on 
the forward channel is 

L{x) = i[H,x] + v;x% - ^-{v;v„,x}. (3.22) 

Therefore the diagram for resonance fluorescence (3.20) is equivalent to 

T'h 

B B 

1(S)1 
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where Ut is given by St StUt for t > 0. For h{t) = —i-^, we find tlic master equation 
for resonance fluorescence (3.1). From now on we will no longer suppress the oscillations 
of the laser, i.e. we take h{t) = —iexp{iLut)-^. Then we find 

L{X) = i[H,X] - z^[e-^'^*y + e''^*V*,X] - ^{V*V,X} + V*XV, 
note that the laser is resonant when lo = coq. 



3.7 Belavkin's stochastic Schrodinger equations 

Now we are ready to derive a stochastic differential equation for the process £^ (X) . In the 
next section we will see that this equation leads to the stochastic Schrodinger equations 
(3.5) and (3.9), that we already encountered in Sections 3.2 and 3.3. 

Definition 3.7.1: Let X be an clement of B := Define the process {M^}t>o in the 
algebra Coo — -^°°(f^, S,Pp), generated by the observed process {Yt}t>o (see Section 3.5) 

by 

:=f*(X)-fO(X)- fs^{L{X))dr, 
Jo 

where L : B B is the Liouvillian. In the following we suppress the superscript X in 
M/*- to simplify our notation. 

Note that from the above definition it is clear that Mf is an element of Ct for all t > 0. 
The following theorem first appeared (in a more general form and with a different proof) 
in [17] and is at the heart of quantum filtering theory. We prove it using the properties of 
conditional expectations. For simplicity we have restricted to observing a process in the 
field W^''. The theory can be extended to processes that are in B0 W®^ , transforming it 
into a more interesting filtering theory. For the stochastic Schrodinger equations arising 
in quantum optics our approach is general enough. 

Theorem 3.7.2: The process {Mt}t>o of deEnition 3.7.1 is a martingale with respect to 
the filtration {St}t>o ofO, and the measure Fp, i.e. for allt > s > we have: E*(Mt) = Mg. 

Proof. From the module property of the conditional expectation it follows that E*(Mt) 
= Mg for i > s > is equivalent to E^^Mt — Mg) = for t > s > 0. This means we have 
to prove for alH > s > and E G Sg^ 

/ WAMt-M,){u;)Fp{dw) = Q, 

JE 
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which, by the tower property, is equivalent to 

/ (Mi - Ms)iiu)¥p{dco) = 0, (3.23) 
Je 

i.e. Ep(xf;(Mt — Ms)) = 0. Now using Definition 3.7.1 and again the module property of 
the conditional expectation we find, writing E also for the projection corresponding to xe 

Ep{xE{Mt - Ms)) = p'^(8\X ®E)- S'{X ® E) - j S'{L{X) E)dr'^ 

J S 

= p\X ®E)- p'{X ®E)- p'\L{X) E)dr. 

This means we have to prove: dp*{X ® E) — p^{L{X) ® E)dt = 0, for all t > s. Note 
that p\X ®E) = p^{U^X ® EUt) = p® (t>®^{U^X ® EUt). Therefore dp\X E) = 
p(^(t)'^^{d{UtX®EUt)). We will use the notation below Theorem 3.6.2 with Zi = Ut and 
Z2 = X EUt- Using the quantum Ito table and the fact that only the dt terms survive 
after taking a vacuum expectation, we find: 

cip°(f/;X®^[/i) =pO([l]) +p0([2]) +p0([12]), where 
pO([l]) + p0([2]) = p\u;{i[H,X] ®E- l{V*Vj,X} ® E)Ut)dt 
p°([12]) = p%U:iV*XVj) ® EUt)dt. 

This means dp^{X ® E) = p*(L(X) (g) E)dt, for all t>s, proving the theorem. □ 

Note that in the proof of the above theorem we have used that the projection E & Cs 
commutes with the increments dAj (s) , dA* (s) , ds and with the processes in front of the 
increments in equation (3.15), i.e. Vj, Vj, V^Vj and H. If the theory is extended to a 

more general filtering theory [17], then these requirements become real restrictions on the 
process {Yt}t>o- If they are satisfied the observed process {Yt}t>o is said to be self non 
demolition [17]. 

Definition 3.7.1 implies the following stochastic differential equation for the process S*{X) 

d£\X) = S\L{X))dt + dMt, (3.24) 

called the Belavkin equation. The only thing that remains to be done is linking the 
increment dMt to the increment of the observed process dYj. 

Let us assume that the observed process {Yt}t>o satisfies a quantum stochastic differential 
equation 

dYt = aj{t)dA*{t) + Pij{t)dAij{t) + a*it)dAj{t) + S{t)dt, 

for some adapted stochastically integrable processes aj,Pij, and S, such that aj{t), Pij{t), 
6{t) G W?*^ for aU t>0, and P*j = (3ji, S = S* since Yf is selfadjoint. Furthermore, since 
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the observed process {lt}t>o is commutative, we have [c/YjjFs] = for all s < t, which 
leads to 

[a, (t) , Y^dA* (t) + [A,- (t) , y,] dA,,- (t ) + [a* {t) , Ys]dAj (t) + [6{t),Ys]dt = => 
[aj{t),Ys]=0, [Pij{t),Ys]=0, [a*{t),Ys]=0, m,Ys]=0, 

i.e. aj{t), Pij{t),aj{t),S{t) G At- This enables us to define a process Yt by 

dYt = (aj{t)dA*{t) - £'{V*aj{t))dt\ + (/3ij{t)dAij{t) - £'{V*Pij{t)Vj)dt\ + 
[a* {t)dAj it) - (a* {t)Vj)dtj , Yo = 0, 

i.e. we have the following splitting of Yj: 

Yt = Yo + Yt + 1^ (^£'{V*aj{s)) + £'{V*f3ij{s)Vj) + S'{a*{s)Vj) + S{s))ds, (3.26) 

which in view of the following theorem is the scmi-martingale splitting of Yj. The process 
Yf is called the innovating martingale of the observed process Yj. 

Theorem 3.7.3: The process {Yt}t>o is a martingale with respect to the filtration {'Bt}t>o 
of Q and the measure Pp, i.e. for allt> s >0 we have: ^p{Yt) = Yg- 



Proof. We need to prove that for alH > s > : ^p{Yt — Yg) = 0. This means we have to 
prove for alH > s > and E G E^: 

J^Ki^t - Ys){uj)¥p{duj) = ^ JjYt - Ys){u;)Fp{du;) = ^ Ep(YtE - Y^E 

- f^[s^{v;a^{r))E + S^{v:p,^{r)V^)E + 8\a]{r)V^)E^ =0 ^ 

p\YtE) - p'{YsE) = 

^ p'(£'\V*aj{r))E + 8\V* l3ij{r)Vj) E + {a*{r)Vj)E + <5(r)E)dr. 

For t = s this is okay, so it remains to be shown that for alH > s > and £^ G S^: 

dp\YtE) = {e* {V*aj (t)) E + £' {V*Pij it)Vj)E + {a* it)Vj)E + 6it)E^ dt ^ 
dp\Ut*YtEUt) = p*(£\V*aj{t))E + £^ {V* Pij{t)Vj)E + £\a*{t)Vj)E + 

We define: Zi{t) := , Z2{t) := YfE and Z^{t) := Ut then we find, using the notation 
below Theorem 3.6.2: dp^{U^YtEUt) = + [2] + [3] + [12] + [13] + [23] + [123]). Re- 

member p^ = p® (j)®^ , i.e. we are only interested in the dt terms, since the vacuum kills 
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all other terms. The terms [1], [3] and [13] together make up the usual Lindblad term and 
since L(l) = we do not have to consider them. 

Furthermore, term [2] contributes U^d{t)EUtdt, term [12] contributes U^V*aj{t)E Utdt, 
term [23] contributes U^a*{t)VjEUtdt and term [123] contributes U^V*pij{t)Vj Utdt, 
therefore we get 

dp\u:YtEUt) = 

p'^{U*a*{t)VjEUt + U*V* f3ij{t)VjUt + U*V*aj{t)EUt + U;5{t)EUt)dt = 
p^{a*{t)VjE + V*pij{t)Vj + V*aj{t)E + S{t)E)dt = 

* (y/Q, {t))E + {V*p.,, {t)Vj)E + 8' {a* {t)Vj)E + 5{t)E) dt, 
proving the theorem. □ 



Remark. In the probability literature an adapted process which can be written as the 
sum of a martingale and a finite variation process is called a semimartingale [77]. The 
Theorems 3.7.2 and 3.7.3 show that Mt and Yt are semimartingales. 

We now represent the martingale Mt from Definition 3.7.1 as an integral over the inno- 
vating martingale (cf. [57]) by 

dMt = VtdYt (3.27) 

for some stochastically integrable process rjt, which together with equation (3.26) provides 
the link between dMt and dYt. We are left with the problem of determining r]t, which we 
will carry out in the next section for the examples of Section 3.2 and 3.3. Here we just 
give the recipe for finding rjt- 

Recipe. Define for all integrable adapted processes bt and cj a process Bt in by 

dBt = btdYt + ctdt. (3.28) 

These processes form a dense subalgebra of Coo- Now determine rjt from the fact that 6^ 
leaves p* invariant [17], i.e. for all Bt 

p\S\BtX)) = p\BtX). 

Prom this it follows that for all Bt 

dp''{UtBt{£\X) - X)Ut) = 0. (3.29) 

We evaluate the differential d(U^Bt{£^{X) — X)Ut) using the quantum Ito rules. Since 
p^ = p® (f)^'' we can restrict to the dt terms, since the others die on the vacuum. We will 
use the notation below Theorem 3.6.2 with Zi{t) = U^, Z2{t) = Bt, Z3{t) = £\X) - X 
and Ziit) = Ut- The following lemma simplifies the calculation considerably. 
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Lemma 3.7.4: The sum of all terms in which Z2 is not differentiated has zero expectation: 
p%[l] + [3] + [4] + [13] + [14] + [34] + [134]) = 0. 

Proof. The dt terms of [3] are Bt8\L{X))Utdt and -U^ Btr]t{£\V*aj) + £\V*Pij 
Vj) + £^{a*Vj)^Utdt. Using the fact that £"* leaves invariant we see that the term 
UlBt£^(^L{X)^Utdt cancels against the dt terms of [1],[4] and [14], which make up the 
Lindblad generator L with a minus sign. The other term of [3] is cancelled in expectation 
against the dt terms of [13], [34] and [134], since 

p°([13]) = p\Bt7]tV*a,)dt = p\£\Btr]tV*aj))dt = {BtVt£\V*aj))dt 

p0([34]) = p\BtVta*Vj)dt = p'{£\Bt7]ta*Vj))dt = p' {Bt7]t£\a*Vj)) dt 

p0([134]) = p\BtmV*l3ijVj)dt = p\£\BtritV*l3ijVj))dt = p'{Btm£\V* pijVj))dt. 

□ 



Using equation (3.25), the fact that £* leaves invariant and the module property, we 
find that the term [2] has expectation zero as well 

p\[2])=p' (btdYti£\X)-X)) = 

- p' {bt£\V*aj{t) + a*{t)Vj + V* PijVj){£\X) - X)) dt = 

- p* {bt£\V*aj{t) + a*{t)Vj + V* PijVj)£\£\X) - X)) dt = 0. 

Thus, only the terms containing no Bt nor ct can contribute non-trivially. This leads to 
an equation allowing us to obtain an expression for rjt by solving 

p°([12] + [23] + [24] + [123] + [124] + [234] + [1234]) = 0. (3.30) 

Although this can be carried out in full generality, we will provide the solution only for our 
main examples, the photon counting and homodyne detection experiments for a resonance 
fluorescence setup, in the next section. 



3.8 Examples 

We now return to the example considered in Section 3.2. We were considering a 2-level 
atom in interaction with the electromagnetic field. The interaction was given by a cocycle 
Ut satisfying equation (3.21). The observed process is the number operator in the side 
channel, i.e. Yt = Ass{t). Therefore dYt = dAss{t) - £\V*Vs)dt. Recall now the notation 
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Zi{t) = U^, Z2{t) = Bt,Z2,{t) = £\X) - X and Zi{t) = Ut, their differentials are given 

by 

dU: = U* {KdAait) - V,dA*,{t) - i-iH + ^V:V,)dt} 

dBt = btdAssit) + (ct - bt£\v;Vs))dt 
d{£\X) -X) = ritdKss{t) + {£\L{X)) - rit£\V*Vs))dt 

dUt = [y'^dAlit) - V;dA„{t) - {iH + ^-v;v„)dt]Ut. 

Following the recipe of the previous section we now only have to determine the dt terms 
of [12], [23], [24], [124], [123], [124] and [1234]. Ah of these terms are zero in expectation 
with respect to except for [124] and [1234] 

p\[m) = p^(u;btv:{£\x) - x)VsUt)dt 
p°([1234]) = p''{u;btvtv:VsUt)dt. 

For all bt the sum of these terms has to be in expectation, i.e. 

V6t : (^bt (v: {£' {X) -X)Vs + 7?*^ Vs)jdt = ^ 
V6t : (^£' (^bt(v:{£\X) -X)Vs + ^V/F,) j j dt = ^ 
ybt : p* (^bt (£\X)£\v:Vs) - £\v:XVs) + rit£\v:Vs))j dt = 

Substituting the expressions for r]t and Yj into equation (3.24) we obtain the Belavkin 
equation for photon counting in the side channel 

d£\X) = £\L{X))dt + {^^^^ - {dKs{t) - £\v:Vs)dt). (3.31) 

Now recall that £^{X) = p*(X,), i.e. it is the function — C : ui ^ p^(Xt^). For ah 
X ^ B = M2, the M2 valued function X, is the constant function uj X. Therefore for 
all X in B, the Belavkin equation (3.31) is equivalent to 

dpiiX) = p\{L{X))dt + {^^^ - {dkssit) - pi{V:Vs)dt), 
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which is equivalent to the Belavkin equation of Section 3.2, equation (3.5). In simulating 
the above equation we can take for Yt = Ass{t) the unique jump process with independent 
jumps and rate pi{V*Vs), since Agsit) — Jq pl{V*Vs)dr has to be a martingale. 

Let us now turn to the homodyne detection scheme which we already discussed in Section 
3.3. The observed process is now Yt = X^{t) = AKft) + As{ft) (see Section 3.3 for the 
definition of ft). This means the innovating martingale Yt satisfies dYt = e*'^*c?A*(i) + 
e-''l'*dAs{t) - £\e''t'^V* + e-''f'*Vs)dt, where 0t = 0o + ujiot with uio the frequency of the 
local oscillator. Therefore we find different differentials for Bt and £^{X) — X than we 
had in the photon counting case 

dBt = bt{e''^*dAl{t) + e-''l'*dA,{t)) + (q - bt£\e''l'*V; + e-'^*Vs))dt 
d{£\X) -X) = r,t{e''f'*dA*{t) + e-^'^*d^(i)) + 

(£\L{X)) - rit£\e'^*V* + 6-'"^*^,))^* 

Following the recipe of the previous section we now only have to determine the dt terms 
of [12], [23], [24], [124], [123], [124] and [1234]. Ah of these terms are zero in expectation 
with respect to p^, except for [12], [23] and [24] 

p0([12]) = p'{u;e''t'^V:bt{£\X) - X)Ut)dt 

p%[23]) = p\Ut*btVtUt)dt 

p0([24]) = pO(c/;6t(f*(X) - X)e-"l'*VMt)dt. 

For all bt the sum of these terms has to be in expectation, i.e. 



V6t : p* (e'^'y;(^:*(X) -X) + {£\X) - X)e-'^'Vs + r?*) j = 

ybt : p' (^£' (^bt (e'^' V: {£' {X)-X) + {£' {X) - X) e"^'^* + 77*) j j = ^ 

V6t : p^(bt{ - £\e'^W*X + e-''>>'XVs) + 

£\e}'t''V* + e-'^'Vs)£^{X) + r]t))dt = ^ 
= £\e'^'V;X + (T''^' XVs) - £\e'^'V; + e-'^Ws)£\X). 

Substituting the expressions for -qt and Yt into equation (3.24) we obtain the Belavkin 
equation for the homodyne detection scheme 

d£\X) = £\L{X))dt + {£\e''^*V*X + e-''^*XVs) - £\e''t'*V* + e-''t'Ws)£\X)) x 
X {e''>'*dAl{t) + e-''^*dAs{t) - £\e''^*V* + e-''^*Vs)dt) . 

(3.32) 
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Now recall that £^{X) = p^X,), i.e. it is the function $7^ ^ C : lo /9^(X^). For all 
X £ B = M2, the M2 valued function X, is the constant function lu X. Therefore for 
all X in B, the Belavkin equation (3.32) is equivalent to 

dpiiX) = pi{L{X))dt + {p\{e'^W:X + e-'^'XVs) - pUe'^'F/ + e-'^^Vs)pi{X)) x 
X {e'^'dAl{t) + e-''t''dAs{t) - p\{e'^'V; + e-'^'Vs)dt), 

which is equivalent to the Belavkin equation of Section 3.3, equation (3.9). Since A* {ft) + 
^sift) — /o + e~'^'^''Vs)dr is a martingale with variance t on the space of the 

Wiener process, it must be the Wiener process itself. 
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Abstract^ 



We study an open system in contact with its environment, the electromagnetic field. The infor- 
mation gained by measuring a quadrature of the field is used to send control pulses to the system. 
Goal is to fix the unknown state of the system in time. We show that in the special case of an 
essentially commutative interaction this goal can be achieved. In dealing with spontaneous decay 
we approximate the essentially commutative situation by bringing the field in a squeezed state. 
We show that when squeezing goes to infinity, the state can again be kept fixed. 



4.1 Introduction 

The last decade there have been rapid developments in quantum information theory, ini- 
tiated mainly by some fundamental papers [79], [39] showing the increased possibilities 
when quantum features are exploited in computations. However, implementation of the 

^This chapter is based on [22]. 
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proposed algorithms on real physical qubits still poses a great challenge. One of the 
problems is the interaction with the environment, i.e. the electromagnetic field, and the 
decoherence that goes along with it. Dealing with this problem motivates the development 
of theory and methods for coherently manipulating, or controlling, quantum systems. 

Decoherence is a result of ignoring information lost from an open quantum system to 
its environment via their interaction. However, the lost information can be retrieved, 
at least partially, by observing the environment, i.e. by performing measurements on it. 
The decoherence can be combatted by using the retrieved information in a scheme for 
controlling the quantum system, see also [51]. 

Since the electromagnetic field and the open system are in interaction, information on the 
system itself is gained when measuring some obscrvables of the field. Hence conditioning 
on the obtained measurement results provides a back-action of the measurement in the 
field on the open system. One of the pioneers in this area is Belavkin who extended many 
ideas in classical filtering theory, cf. [83], to the quantum regime [11], [17]. Quantum 
filtering theory [14], [17] explains how the state, conditioned on the result of a continuous 
time measurement in the environment, evolves in time. Note that since the results of 
the continuous time measurement are random, the conditioned state is also a random 
state. Quantum filtering theory provides a stochastic differential equation, the Belavkin 
equation, for the state evolution in which the measurement process is one of its driving 
terms [14], [17]. 

Another approach to the back-action due to conditioning, is via quantum trajectory the- 
ory as developed in quantum optics in the late 1980's and early 1990's [27], but already 
envisioned by Davies [82], [33] in the 1970's. In this approach photon counting measure- 
ments are analysed to obtain a continuous time evolution of the open system interrupted 
by jumps the moments at which photons are detected. Differentiation of the trajectory 
evolution leads to a stochastic Schrddinger equation [28], which is a stochastic differential 
equation for the evolution of the state conditioned on the outcomes of the counting exper- 
iment. A diffusive limit of photon counting in which the jumps in the state space decrease 
in size but become increasingly frequent, makes it possible to incorporate homodyne and 
heterodyne detection schemes into quantum trajectory theory [9], [27], [92]. The stochas- 
tic Schrodinger equations encountered in quantum optics are equivalent to the Belavkin 
equations from quantum filtering theory [17], [23]. 

The result of the continuous time measurement in the field can be used to exert control 
over the system. The solution to the quantum filtering problem [11], [14] makes it possible 
to directly carry over many ideas in classical control theory [83], [65], [66] to the quantum 
regime, [12], [13], [14], [18], [36], [86]. Coming from the quantum trajectory approach, 
other pionering work in quantum control was done by Wiseman and Milburn in the first 
half of the 1990's, see [91], [90], [93]. Two different objectives in control problems can be 
distinguished, one where the state is controlled in order to let it follow a certain path in 
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time [18], [36], and one where the semigroup describing the dissipative evolution of the 
open system, i.e. the channel itself, is being controhed [71], [51], [3], [4]. 

In this paper a problem of the second type is considered. The question addressed here is 
how to keep an unknown state of an open system fixed in time, i.e. how to keep its dy- 
namical semigroup as close to identity as possible. In this article we will not be concerned 
with optimality results. The main issue is to find or engineer situations where the control 
is perfect, in which case the control scheme is said to restore quantum information [51]. 
Furthermore, we will not be concerned here with encoding our system into the code space 
of a larger system and then protecting just this code space [5], [4]. 

The control scheme consists of two parts. The first part is an evolution over a period of 
r time units in which a quadrature of the field is observed. This evolution is governed by 
the Belavkin equation corresponding to this measurement. In the second part the result 
of the measurement is used to construct a laser pulse designed, if at all possible, to take 
the system through a Rabi cycle that corrects the evolution of the past r time units. This 
scheme is studied in the limit for very small r, i.e. the control pulses are sent at very high 
frequency. 

In general the above control scheme will not be able to restore quantum information. 
Since the interaction of the field and the system is studied in the weak coupling limit, the 
field acts as two classical noises. However these two noises, represented by two different 
quadratures of the field, do not commute with each other. Therefore only one of these 
noises can be observed and its disturbing effect on the system corrected. An idealised 
interaction of system and field in which there is only one instead of two classical noises 
present is called essentially commutative [63]. In the essentially commutative case it will 
turn out that the above control scheme restores quantum information. 

For the more realistic situation where both noises are present our strategy will be to 
manipulate the state of the field in order to approximate the essentially commutative 
case. This is done by putting the field in a squeezed state, i.e. one quadrature's variance 
increases while the other quadrature's variance decreases, [42], [43], [58]. The idea is to 
measure the noise with the large variance and correct its disturbing effect on the system. 
It will turn out that when squeezing goes to infinity the control scheme described above 
will restore quantum information. 

This paper is organised as follows. Section 4.2 describes the dissipative evolution of the 
open system within the Markov approximation. The joint evolution of system and field 
is given by unitaries satisfying a quantum stochastic differential equation in the sense of 
[54]. In the next section a brief exposition of quantum stochastic calculus [54] is given. 
This enables us to make sense of the quantum stochastic differential equation providing 
the unitaries of section 4.2. Sections 4.2 and 4.3 describe dilation theory and quantum 
stochastic calculus in a nutshell. Section 4.4 is a brief exposition of quantum filtering 
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theory. It contains a derivation of the Belavkin equation for field quadrature measurement. 

Sections 4.5 and 4.6 deal with controlling the state of an open system in the essentially 
commutative case and the more realistic situation of spontaneous decay of a two-level sys- 
tem, respectively. Here we show that for the essentially commutative case it is possible to 
restore quantum information. For spontaneous decay, however, problems are encountered 
motivating the investigation in the remainder of the paper. 

Section 4.7 shows how to describe the interaction of system and field when the field is 
in a squeezed state. To do this we have to do quantum stochastic calculus in the GNS- 
representation space of the squeezed state. In the last section the Belavkin equation for 
measuring a quadrature of a squeezed field is given and a control scheme based on this 
measurement is presented. It turns out that when squeezing goes to infinity, the scheme 
restores quantum information. 



4.2 The dilation 



Let B := M„ stand for the algebra of obscrvablcs of an n-dimensional quantum system. 
On this algebra {Tt}t>Q is a semigroup of completely positive identity preserving opera- 
tors. It represents the irreversible time evolution of the system in the Heisenberg picture. 
Lindblad's theorem [67] asserts that = exp(tL) where the generator L is given by 

k 

L{x) = i[H,x] + Y^v;xvj - :^{V;Vj,X}, XeB, 

i=i 

with H and the VjS fixed elements of B, H being selfadjoint. The notation {X, Y} stands 
for the anticommutator XY + YX. For simplicity, we take k = 1 and H = 0, i.e. 

L{X) = V*XV -^{V*V,X}. (4.1) 

This paper deals mainly with two special cases of the above situation. In the first special 
case we have cither V = V* ot V = —V* . This case is called essentially commutative [63], 
see section 4.5. In the second special case we have 

.^(;»). 

Then the semigroup Tt describes spontaneous decay to the ground state of a two-level 
system, see sections 4.6 and 4.8. 

The system B and its environment, the electromagnetic field, evolve reversibly in time. 
The irreversible evolution Tf of B is the result after tracing out the field. Up to section 4.7 



4.2. THE DILATION 



69 



the electromagnetic field to which the system B is coupled, will be taken in the vacuum 
state or a coherent state. Then, see section 4.7 for more details, a decay channel in the 
field can be modelled by the bosonic or symmetric Fock space over the Hilbert space (M) 
of square integrable wave functions on the real line, i.e. 

oo 

J^:=C©0l2(M)®»". 

n=l 

The algebra generated by the field observables on contains all bounded operators and 
it is denoted by W. 

For future convenience we already distinguish two decay channels in the field, i.e. we 
rewrite L as 

L{X) = VjXVf - \{VlVf, X} + V^XVs - ^{V:Vs,X}, XeB, (4.3) 

where Vf = KfV, Vs = KsV and Kf,Ks G M such that + = 1. The subscripts / 
and s stand for forward and side channel, respectively. On the forward channel in the field 
we will put a laser with which we want to control the system, while in the side channel of 
the field we are going to perform a measurement. The decay rates into the forward and 
side channel are given by and respectively. Since the field is modelled by these 
two decay channels, we need two copies of the algebra W, denoted (8) W*. 

The free evolution of a channel in the field is given by the unitary group St, the second 
quantization of the left shift s{t) on L^(M), i.e. s{t) : f ^ f{- + t). In the Heisenberg 
picture the evolution on ® is 

W ^{St® StYW{St ® St) ■■= Ad[St (g) St]{W), W eW^ (E) W. 

The system B and field together form a closed system, thus their joint evolution is given 
by a one-parameter group {Tt}teM. of *-automorphisms on B (E) W'^ (E) 

X ^U;XUt:= Ad[Ut]{X), X eB(^W^ (^W. 

The group Ut is a perturbation of the free evolution without interaction. We describe this 
perturbation by the family of unitaries Ut := {S-t ® S-t)Ut for all i G M satisfying the 
cocycle identity 

Ut+s = iS-s<^S-s)Ut{Ss®Ss)Us, for ah t,seR. 

The direct connection between the reduced evolution of B given by (4.3) and the cocycle 
Ut is one of the important results of quantum stochastic calculus [54] which is the object 
of the next section. For the moment we only mention that in the weak coupling limit [2], 
Ut is the solution of the stochastic differential equation [54] , [74] , [72] 



dUt = {VfdA}{t) - Vf*dAf{t) + VsdAlit) - V;dAs{t) - ^V*Vdt}Ut, Uq = 1. (4.4) 
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We will see in the next section that if Ut satisfies (4.4) the following dilation diagram [59], 
[60] commutes: 

Tt 



B 



B 



Idigiiigii 



(4.5) 



i.e. for all X e B : Tt{X) = (id (g) ^ (g) ^) {ft{X (g) 1 (g) 1)) , where is the vacuum state 
on W, and 1 is the identity operator in W. Any dilation of the semigroup Tf with Bosc 
fields is unitarily equivalent with the above one under certain minimality requirements. 

The dilation diagram can also be read in the Schrodinger picture if we reverse the arrows: 
start with a state p of the system B in the upper right hand corner, then this state 
undergoes the following sequence of maps 



This means that at i = 0, the atom in the state p is coupled to the electromagnetic field 
in the vacuum state, after t seconds of unitary evolution the partial trace over the field is 
taken. 



4.3 Quantum stochastic calculus 



Here, we briefly discuss the quantum stochastic calculus developed by Hudson and Par- 
thasarathy [54]. For a detailed treatment of the subject we refer to [74] and [72]. The 
exposition here is a bit broader than strictly necessary for the construction of the cocycle 
of the previous section. However, the general description [74] presented here is needed in 
section 4.7. 



Let 7i be a Hilbert space. We define the bosonic or symmetric Fock space over Ti by 

oo 

fc=i 

In the previous section we had 71 = L^(M). For every f & Ti. we define the exponential 
vector e(/) G ^(W) in the following way 
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The inner product of two exponential vectors e(/) and e{g) is {e{f), e{g)) = cxp{{f, g)). 
We denote the vacuum vector e(0) = 1 © 0© 0© . . . also by The span of all exponential 
vectors, denoted V, forms a dense subspace of TiH). 

Let ^ be a projection (on the Hilbert space H) valued measure on M with no jump 
points, i.e. = for all t G M. Denote by 'Ht\^,'Hys,i\ and Tiyi the ranges of the 

projections ^((— oo, t]), ^([s, t]) and ^([t, oo)), respectively. For a vector f eH we denote 
ft] ■= ^{{-oo, t])f,f[s,t] ■= ^{[s,t])f and f[t := ^{[t,oo))f. Let us write H as the direct 
sum Ti.^ © Wff, then J-{T-i) is unitarily equivalent with J^{Ti.^) ^TiTi^t) through the identi- 
fication e(/) = e(/j] ) (8) e For notational convenience the tensor product signs between 
exponential vectors are often omitted. The algebra W := B(^J^{H)) also splits as a tensor 
product Wt] Wit where Wt] := and := 

A map m : M+ 7i : t ^ rrit is called a (^-martingale if rrit € TC^ for all t and 
^([0, s])mj = nis for all s < t. For m and m' ^-martingales, there exists a complex 
valued measure (of finite variation on every bounded interval), denoted {{m,m')) on M+, 
satisfying 

{{m,m')){[0,t]) ={mt,m't), (4.6) 

for all t > 0. Let m be a ^-martingale. The annihilation operator A{mt) and creation 
operator A*{mt) are defined on the domain V by 

A{mt)e{g) = {mt, g)e{g), g eH, 

{e{h),A*{mt)e{g))^,^.= {h,mt){e{h),e{g))^,^y h,gen. 



Let Mt be one of the processes A{mt) or A* (nit) for some ^-martingale m. The following 
factorisation property [54], [74] makes the definition of stochastic integration against Mt 
possible 

{Mt - Ms)e{f) = e{f,^){{Mt - MMfls,t])}e{fit 



with {Mt — Ms)e{fis,t]) ^ -^(^[s.t])- We first define the stochastic integral for the so-called 
simple operator processes with values in the atom and noise algebra W where B := Mn 
and W is the algebra of all bounded operators on the Fock space T{H). 

Definition 4.3.1: Let {Ls}o<s<t be an adapted (i.e. Lg G B®B{T{Hg^)) for all < s < i) 
simple process with respect to the partition {sq = 0, = t} in the sense that 
Lg = Lgj whenever Sj < s < Sj+i. Then the stochastic integral of L with respect to M on 
I V is given by [54], [74]: 



-in 



/ LsdMs fe{u) := ^ (L,./e(u,.])) ((M,.^, - M,.)e(n[,.,,.^^]))e(ti[, J. 
•^0 j=o 



72 



CHAPTER 4. SQUEEZING ENHANCED CONTROL 



By the usual approximation by simple processes we can extend the definition of the 
stochastic integral to a large class of stochastically integrable processes [54], [74]. We 
simplify our notation by writing dXf = LfdMt for Xt = Xq + LgdHg. Note that the 
definition of the stochastic integral implies that the increments dMg lie in the future, i.e. 
dMs G y^[s- Another consequence of the definition of the stochastic integral is that its 
expectation with respect to the vacuum state • <&) is always due to the fact that the 
increments dA{mt) and dA*{mt) have zero expectation values in the vacuum. This will 
often simplify calculations of expectations, our strategy being that of trying to bring these 
increments to act on the vacuum state thus eliminating a large number of differentials. 

The following theorem of Hudson and Parthasarathy extends the Ito rule of classical 
probability theory. 

Theorem 4.3.2: (Quantum ltd rule [54], [74]) Let Mi and M2 each he one of the 
processes A(mt) or A*(m[). Then M1M2 is an adapted process satisfying the relation: 

d{MiM2) = MidM2 + M2dMi + dMidM2, 

where dMidM2 is given by the quantum ltd table: 



dMi\dM2 


dA*{m'^) 


dA{m'^) 


dA*{mt) 








dA{mt) 


d{{m, m')) 






Notation. The quantum Ito rule will be used for calculating difi^erentials of products of 
Ito integrals. Let {^j}j=i,...,p be Ito integrals, then 

diZiZ2...Zp)= Yl M 

uC{l,...,p} 

where the sum runs over all non-empty subsets of {1,. . . ,p} and for any 1/ = {ii,. . . i^}, the 
term [u] is the contribution to d{ZiZ2 ■ ■ ■ Zp) coming from differentiating only the terms 

with indices in the set {ii, . . . i^} and preserving the order of the factors in the product. 
For example the diff'erential d{ZiZ2Z^) contains terms of the type [2] = Zi{dZ2)Z^, [13] = 
(dZi)Z2(dZ3), and [123] = {dZi){dZ2){dZ^). 

Let us return to the setup of section 4.2. We now make sense of equation (4.4). Note that 
the Hilbert space TL is L^(IR) © L^(]R). The forward and side channel both have their own 
copy of L^(]R). The projection valued measure ^ is given by 



miff ® fs) = {XI ff) e ixifs), ffjs e L\R), 
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for all Borcl subsets / of R. Here xi denotes the indicator function of /, i.e. the function 
that takes the value 1 on I and is elsewhere. 

The maps mf : ^ H : t X[o,t] ® and : R+ ^ H : t X[o,t] 

are ^-martingales. Wc denote the annihilation A{m(), A{mf) and creation operators 
A*{m{),A*{ml), defined on V by (4.7), more compactly by Af{t), As{t), A*j:{t) and 
respectively. The calculus for stochastic integrals with respect to Ao-(t) and Al(t), a, G 
{/, s} is then given by the Hudson-Parthasarathy ltd table [54], [74]: 



dMi\dM2 


dAlit) 


dA^{t) 


dA*{t) 








dA^{t) 


5 audi 






Let us introduce the selfadjoint quantum noise f3t describing the interaction between the 
quantum system B = M„(C) and the electromagnetic field 

dpt := -i {VfdA} (t) - Vf*dA f (t) + VsdA* (t) - V^dA^ (t)), Po = 0. (4.8) 

It is clear in our example of spontaneous decay of a two-level system that this noise 
represents an interaction consisting of creations of excitations of the two-level system 
accompanied by annihilations of photons in the decay channels and vice versa. It describes 
the interaction of the electromagnetic field, in which we distinguished two decay channels, 
and the two-level system in the weak coupling limit [2]. We let the cocycle Ut of section 
4.2, providing the evolution in the weak coupling limit of the two-level system and field 
together, be given by the quantum stochastic differential equation 

dUt = {idpt - \{d(3t?}Ut = 

{VfdA}{t) - Vf*dAf{t) + V.dAlit) - V^dA.it) - ^V*Vdt}Ut, 
Uo = l. 

We can now check that the dilation diagram (4.5) commutes. Using the continuous tensor 
product structure of the Fock space J^{H), it is easy to see that following the lower part 
of diagram (4.5) defines a semigroup on B, i.e. we only have to show that it is generated 
by the Lindblad operator L of equation (4.3). For all X & B 

did (g) (/) (g) (f){ft{X (8) 1 ® 1)) = Id ® (g) (f){d{U^{X 1® l)Ut)) . 

Using the notation below Theorem 4.3.2 with Zi = and Z2 = {X (g) 1 (g) l)Ut, we find 



did (g (g (f){ft{X (g 1 (g 1)) = Id (g (g (^([1] + [2] + [12]) . 
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With the aid of the Hudson-Parthasarathy Ito table we can evaluate these terms. We 
are only interested in the dt-tcrms since the expectation with respect to the vacuum kills 
the other terms. The terms [1] and [2] provide the anticommutators —^{VjVf,X}dt and 
— ^{V*Vs,X}dt and [12] provides the terms V^XVfdt and V*XVsdt, proving our claim. 



We now change the situation in diagram (4.5) by introducing a laser on the forward 
channel, i.e. the forward channel is now in a coherent state 7/1 := (-(/'(/i), • V'(^)) where 
il^{h) := exp(— i||/tp)e(/t), the exponential vector e{h) for some h G L^(R+) normalised 
to unity. The laser will be used to send control-pulses to the system B. This leads to the 
following dilation diagram 



B 



B 



Id(g)i(g)i 



Id(g)7h<5 



(4.9) 



i.e. the evolution on B has changed and it is in general not a semigroup. Denote 
by W{h) the unitary Weyl or displacement operator defined on V by: W{h)il){g) = 
exp{-2ilm.{h, g))tp{g + h). Note that W{h)^ = VF(/i)V'(0) = i^ih), so that we can write 
for all X G H 

Tl'iX) = Id ® 7ft ® 0(C/;X ® 1 ® lUt) = 
Id (g) ® <p{Wf{hyU^X ® 1 (g) lUtWfih)) = 
Id ^(t)^(l){Wf{ht])*U;x 1® lUtWfiht])), 

where /i^j = hx{o,t] aiid ^/(^) := 1 ® W{h) (g) 1. Defining := UtWf{ht]), together with 
the quantum stochastic differential equation for Wf{ht]) [74] 

dWfiht^) = {h{t)dA}{t) - h{t)dAfit) - l\h{t)\'dt}Wfiht^), Wfiho]) = 1, 
and the Ito rules leads to the following quantum stochastic differential equation for 

dU^ = {{Vf + h{t))dA}{t) - [VJ + h{t))dAf{t)) + VsdA*{t) - V*dAs{t) - 

1 , (4.10) 

-{\h{t)\^ + V*V + 2hit)Vf*)dt\Ut\ = 1. 

Therefore, the dilation diagram (4.9) is equivalent to 



B — ^ B 



Idigiiigii 



Id(g)<^(g)</> (4-11) 



B^Wf^W — B^Wf^W 
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In the following, we will often omit the superscript h to simplify the notation. Define a 
Hamiltonian by H := i(/i(i)V/ — h{t)V^), then following the lower part of diagram (4.11) 
and using Ito's rules, see Theorem 4.3.2, shows that the time dependent generator of the 
dissipative evolution T/* in the presence of the laser on the control channel is given by 

L{X) = i[H, X] + V*XV - ^{V*V, X}. (4.12) 
Later on we will choose /i in a suitable way in order to exert control on the system B. 



4.4 The Belavkin equation 

Let us now turn our attention to the side channel. In this channel an observable is mea- 
sured continuously in time. Goal is to briefly show how to derive a stochastic differential 
equation for the stochastic state evolution of the system B conditioned on the outcome 
of the measurement process. The method described below is known as quantum filtering, 
see [17] and [23] for a more detailed treatment. 

In this paper the observable Y-^ of the field that is measured continuously in time will 
always be a field quadrature, i.e. 

F/ := 1b ® Iwf ® {{e-''^Mt) + e''^A*sit)) ® IwfJ e B^W^ (>V| Wf^) , (4.13) 

for some phase (p G [0, 27r). Such a field quadrature measurement can be performed by a 
homodyne detection experiment. See [17], [23] for measurement of other observablcs. Let 
p be the initial state of the quantum system B. We describe the measurement process in 
the interaction picture, i.e. the shift part of Ut := {St'S>St)Ut acts on the observables while 
the cocycle part Ut, given by equation (4.10) with the superscript h suppressed, acts on 
the states 

p\X) := p^^{U;XUt), XeB^W^^W. 

Let Ct be the von Neumann algebra generated by the family of observables {y/; <r <t}. 
Since and Yf commute for all r,t > the algebra Ct is commutative. The algebras 
{Ct}t>o form a growing family, that is Cs C Ct for all s < t. Thus we can define the inductive 
limit Coo '■= fimt-»ooCt, which is the smallest von Neumann algebra containing all Ct. It 
follows via Kolmogorov's extension theorem, see [23] Theorem 5.1, that there exists a 
unique state p°° on C^o which coincides with p* when restricted to Ct C C^o for all t > 0. 
From spectral theory it follows that there exists a measure space (Q, S,Pp) and a growing 
family {T,t}t>o of cr-subalgcbras of S, such that (Coo,/0°°) and {Ct,p^) are isomorphic to 
L°°(ri,S,Pp) and L°°(ri, St,Pp), respectively. The space Jl should be interpreted as the 
paths of the observed process Y^ when the measurement is continued infinitely long. The 
(7-algebras contain the events up to time t. 
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In the Hciscnbcrg picture, when a measurement of an observable Y with discrete spectrum 
Sp{Y) has been performed, all observables in B ^ W-^ ^ have to be sandwiched with 
the projection corresponding to the observed measurement result. If the result of the 
measurement is unknown, but the measurement has taken place, an observable takes the 
form of a direct sum over all possible outcomes of the original observable sandwiched with 
the projections corresponding to the outcomes, i.e. 

^aftermeas.= © ^v^Py XeB^W^^W^. 

yeSp(Y) 

Note that this procedure destroys all coherences between different measurement results. 
Moreover, it maps all observables in B (g) (8) to the commutant of the algebra 
generated by the measured observable. 

Therefore, in analogy with the above, when a process {l^/}o<r<t has been measured con- 
tinuosly in time, we can restrict to the algebra At C B<^W^ <SiW^ which is the commutant 
of the observed process 

At := C't:={XeB^W^<^ XC = CX, VC G Ct}. 

We call At the algebra of observables that are not demolished [17] by observing the process 
{YJ,*}o<r<t- Note that from the double commutant theorem it follows that Ct is the center 
of At 'i.e. Ct = {Ce At; AC = CA, e At}. 

We investigate the situation of the previous paragraph more abstractly for a moment, i.e. 
let ^ be a von Neumann algebra of operators on some Hilbert space H and let C be its 
center. Let p denote a state on the algebra A. We will now explain the decomposition of 
A over its center C, see [56] for all details and proofs. We can identify the center C with 
some L°°(r2,E,P) where P corresponds to the restriction of p to C. The Hilbert space 
H has a direct integral representation H = J® Mi^F{duj) in the sense that there exists a 
family of Hilbert spaces {Hj^jj^gn and for any ijj e M there exists a map u iJjoj & 
such that 

{ij,<P) = I {^P^,<P^)F{dlo) ^,cj)eM. 
Jn 

The von Neumann algebra A has a central decomposition A = J® AujF{duj) in the sense 
that there exists a family {Au)}uje^i of von Neumann algebras with trivial center, or factors, 
and for any A e A there is a map lo <—>■ A^j ^ A^j such that (Atl;)^ = A^jip^^ for all eM 
and P-almost all u> G Q,. The state p on A has a decomposition in states on Au, 
such that for any A G ^ its expectation is obtained by integrating with respect to P the 
expectations of its components A^^: 

piA) = [ p^{A^)F{dj). 
Jn 

Loosely speaking the component A^ G Au, is the operator A e A sandwiched with the 
projection corresponding to a measurement result lo. Moreover, the state p^ is the state 
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p conditioned on the measurement result uj. For all X G ^ we denote by p,(X,) the 
function to ^ pijj{X^^). The complex number Puj{Xi_j) is the expectation of the observable 
X in the state p conditioned on measurement result to. 

Define a map : A C = L°°(J7, S,P) by £p{X) := p,{X,) for all X ^ A. It is easily 
verified, see also [23], that this map is linear, surjective, identity preserving, completely 
positive, it satisfies the module property 

Sp{CiXC2) = CiSp{X)C2, Ci, C2eC, XeA, 

and it leaves the state p invariant, i.e. p[£p{X)) = p{X) for all X ^ A. These properties 
uniquely determine the map £p, see [84]. It is called the conditional expectation of A onto 
C with respect to p. Returning to the original problem, i.e. a whole family of algebras At 
with center Ct, we get a family of conditional expectations £pt : At —>■ Ct- We denote £pt 
more compactly by 

Apart from the family of quantum mechanical conditional expectations there is also 
a family of conditional expectations in the classical sense that plays an important role 
in the following. Denote by E* the unique classical conditional expectation from Coo — 
L°°(17,S,Pp) onto Ct = L°°(n,Et,Pp) that leaves the state or equivalently, the ex- 
pectation with respect to Pp invariant, i.e. p"^ o E* = p°^. These conditional expectations 
satisfy the tower property, that is E^(E*(C)) = E^(C) for ah C £ Coo and t > s > 0. E^ 
is the expectation with respect to Pp and will simply be denoted Ep. Note that the tower 
property for s = is just the invariance of the state p°°{= Ep). 

For all i > and X G ^ the operator X (g, 1 1 e B (gi W'f W commutes with the 
observed process {Y^}o<r<t up to time t, i.e. B C At- Therefore we can define for all 
X eB a process {M^}t>o in the algebra Coo = L'^i^, S,P(t) by 

:= £\X) - £\X) - r r (L(X))dr, (4.14) 
Jo 

where L : B ^ B is the Liouvillian of equation (4.12). Prom the definition it is clear that 

is an element of Ct for all t > 0. The process {M^}t>o is a martingale, i.e. for all 
< s < t we have E^(M/^) = , see [17], [23] for details and a proof. In differential 
form equation (4.14) reads 

dpi{X) = pi {L{X))dt + diWf , 

where we have used that X, is the constant function u ^ X. This equation is the Belavkin 
equation [14], [17], [23]. 

Denote by F/ the process given by the following stochastic differential equation 



dYf = dYf - £\e''^V* + e-''t'Vs)dt, = 0. 
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The process is a martingale, i.e. for all < r < t we have W^{Yf) = Y^, see [17], [23] 
for details and a proof. We call Yf the innovating martingale of the observed process 
y^*. The link between the martingale and the observed process Y^ is provided by the 
martingale representation theorem which states that there exists a stochastically integrable 
process rjf such that 

dMf = nfdY^^ = r,f{dY,^ - £\e^'^V: + e-^^Vs)dt) . 

The process ry^ can be calculated by using that 6* leaves invariant [17]. We refer to 
[23] for the details, the result is 

rif = E\e''t'V*X + e~''t'XVs) - £\e''t'V* + e-''t'Vs)E\X) . 

This leads to the Belavkin equation [17], [23] 

dpi{X) = p\{L{X))dt + (p*.(e*V;X + e-"t'XVs) - p\{e''t>V: + e-*V,)p* (X)) x 

/ X (4-15) 

X Uy^' - pi{e"t'V; + e-'^'Vs)dt\ X e B. 

This equation tells us how the state of the system B evolves over an infinitesimal time dt 
depending on what we observe for the measurement process dYf. Since Y^ is a martingale 
with variance t on the space of the Wiener process, it must be the Wiener process itself. 



4.5 Control: the essentially commutative case 



In this section we focus on dilations that are essentially commutative [63] . We will use the 
results of the measurement of Yff to control the time evolution Tf of the system B in order 
to bring it as close to the identity map as possible. For essentially commutative dilations 
this can be done (nearly) perfectly. This section serves as a guiding example for the more 
realistic situations described in sections 4.6 and 4.8. 

Let V be selfadjoint, i.e. V = V*. The discussion below can easily be adapted to fit the 
situation where V = —V*. Define for a = f,s field observables Yf := i(A'^(t) — A„{t)) G 
W^y Using V = V*, equation (4.4), i.e. the laser on the forward channel is off, simplifies 
to 

dC/t = I - iVfdY/ - iVsdYt' - ^V'^dt^Ut, Uo = l. (4.16) 

This means that for t > the solution Ut is an element of B^Ct, with Ct the commutative 
von Neumann algebra generated by the process {y/ <S) ^*}o<r-<t- (We have dropped the 
extensive notation with the identities tensored to the Y^'s.) This means that we can 
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restrict the dilation of diagram (4.11) to B (8) Coo, i-e. 



B 

Id^i^i 



Tt 



B 



(4.17) 



A dilation for which the relative commutant of the embedding of the algebra B into the 
subalgebra of B(g>Wf generated by {U^X(g>l(g,lUt; X eB,t>0} is commutative, is 
called essentially commutative [63]. Although we restrict the discussion to the essentially 
commutative dilation determined by equation (4.16), the results of this section can be 
extended to all essentially commutative dilations [63] . 



If the dilation is essentially commutative the derivation of the Belavkin equation is ex- 
tremely simple. Since Ut is not demolished by observing {y/}o<r<t, i-e. it is an element of 
the commutant of Ct, wc can just calculate diU^X i^i I® IC^t) using the quantum Ito rules 
and decompose it over the paths of the measurement process. It is clear that this leads 
exactly to the Belavkin equation of the previous section (4.15) with (p = ^ and V = V* 

dpiiX) = pi{L{X))dt + ipi{[Vs,X])dYf, (4.18) 

where L is as in equation (4.12) with H is 0, i.e. there has been no control yet. In general, 

however, we do not have the decomposition of Ut over the center and we have to resort 
to the methods of the previous section. Note that for = ^ and V = V* wc have 
dYf = dYf, i.e. Y^ is the Wiener process. This means that the measurement process is 
non-informative [51], i.e. since here there is no state dependent drift term, we do not gain 
information about the state p* by observing Yf. 

Let be the density matrix of the initial state of the system B. We observe Yf from time 
to time r. Suppose that the laser is off in that time interval, i.e. h{t) = for < t < r. 
Then the stochastic density matrix at time r is given by (4.18) 

p: = p''+ l\piv*dt-^{v*v,pi}dt+i[pi,Vs]dY,', 

where the tilde has been introduced to remind us that this is the state before control has 
taken place. In the time interval from to r we have observed Yf and therefore we can 
determine the difference A(r) := Yf — Yq at time r. Then we want to control the state 
pI with a unitary 

C/J:=exp(zA(r)K), 

i.e. the density matrix after control is given by pi = U^plU^*. This can be done by 
supplying a very sharply peaked laser pulse to the system, i.e. take 

h{t) = -i'^^^^5r{t), < i < 2t, 
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where 5r is the delta function at time r. Then H = — A(t)5t-K' in equation (4.12), i.e. at 
time T all terms in equation (4.12) are negligable with respect to the commutator with H. 
At time r this commutator performs a Rabi oscillation exactly of size U^. After having 
applied the control unitaries the state pi is taken as the new initial state p° and the control 
scheme is repeated after every r time units. 



Note that the control unitary satisfies the following stochastic differential equation 



dC/J = {iVsdY^^ - \vfdT}Ul = UliWsdY^ - \vfdT}, 



= 1. 



Recall that wc have the ltd rules dYfdY/ = dt, dY/dt = dtdYf = 0, dY/dt = dtdY/ 



and dYfdY/ = dY/ dYf = 0. Using the notation below Theorem 4.3.2 with Zi = 
U^, Z2 = pI and Z3 = U^* we find infinitesimally at r = 0, i.e. r should be very small or 
equivalently we should correct with extremely high frequency 



dpi 



T=0 



([1] + [2] + [3] + [12] + [13] + [23] + [123]) 



T=0 



(4.19) 



Note that it immediately follows from Ito's rules that [123] = 0. 



For W & B wc denote by Lw the Lindblad operator corresponding to W acting on density 
matrices p, i.e. 

Lw{p) := WpW* -^{W*W,p}. 

Then we can write ([l] + [3] + [13])|^=o = LvAp^)dT-i[p^ ,Vs]dYo' and ([2] + [12] + [23])|^=o = 
Ly{p^)dT + i[p^, Vs]dYQ — 2Lv^{p^)dT. Therefore we get dp1\^=Q = Ly^[p^)dT and since 
we repeat the control every r time units with r very small, i.e. we take r infinitesimal, 
this leads to the following deterministic state evolution 

dp* = Lvf{p*)dt. 

This means we only have dissipation into the forward channel. We can take Kf arbitrarily 

small which means we have succeeded in freezing the state evolution nearly perfectly, i.e. 
the control scheme restores quantum information in the sense of [51]. 



4.6 Control without squeezing 



We now return to the more realistic situation of spontaneous decay of a two-level atom to 
its ground state. We are again interested in controlling the state of a system in order to 
get as close as possible to freezing its state evolution. However, in trying to do this, we 
encounter problems that motivate the investigation put forward in the sections to come. 
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Guided by the previous section we write V of equation (4.2) as the sum V = Vr + iVi 
with Vr and Vi selfadjoint, i.e. 

Denote for a = /,s : := k^Vr, Vf := k^Vi, Y^{t) := z(A*(t) - A^{t)) and Yf{t) := 
A*{t) + Acr{t). Then equation (4.4), i.e. the laser is off, can be written as 

■= { ( E ^yfdmt) - iyRdYKit)) - ^V*Vdt]Ut, Uo = 1. (4.21) 

Since the noises Y^{t) and Yf(t) in the side channel do not commute we can not observe 
them both simultaneously. 

In the following we choose to observe Y^{t) and to keep notation simple we denote it by 
Yf. The Belavkin equation for observation of Yj follows from equation (4.15) 

dpi = L{pi)dt +i{piv: - Vspi - Tr(p*y; - Vspi)pi) x 
X {dYt - iTv{piv: - Vspi)di) . 

where L is given by equation (4.12) with H = 0. Using the relation p^V* — Vspi = 
[p.) y^] ~ ^{pI) ^/}) this equation simplifies to 

dpi = L{pi)dt + {i[pl Vi] + {pi, Vf} - 2Tr{piVf)pi) {dYt - 2Tr{piVf)dt) . (4.22) 

Note that is a Wiener process plus a stochastic drift term that depends on the state 
of the two-level atom. By observing Yt we can estimate this drift term and in this way 
obtain information about the state 

We run a control scheme similar to the one in the previous section, i.e. wc choose h{t) := 
— i^^^j^(5T-(i) for < i < 2r. Then we get a control unitary U^ = exp (iA(T)y^) , 
satisfying the stochastic differential equation 

dU^ = {iV^dYr - IvfdrjU^ = U^iiV^dYr - ^F^'dr}, C/° = 1. 

The state after control is again given by pi := U^plU^* where pi is given by the Belavkin 
equation (4.22). We use the notation below Theorem 4.3.2 with Zi = U^ , Z2 = pi and 
Z3 = Ul*. For infinitesimal r evaluated at r = 0, this leads to equation (4.19), i.e. 



dpi ^ = ([1] + [2] + [3] + [12] + [13] + [23] + [123]) 



r=0 
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Again [123] = and further ([1] + [3] + [13])|^=o = Lv^{p^)dT - V^jdYo. Furthermore 

we have 

[2] 1^^^ = (pO)dT + {i[p', V^] + {pO, Vf} - 2Tr(pV/)pO) (dYo - 2TV(p0y/)dr) , 

([12] + [23]) 1^^^ = -2Lv^{p°)dT + i[Vi, {p^, Vf}]dT - 2Tr{p%y[V^, p^r. 

A calculation shows that 

Lvip') - Lvsip^) + z[V^, {p', Vf}] = Ly^ip') + Lvs{p') + '-[VfV^ + V^Vf, p% 

Since for spontaneous decay VfV^ + V^Vf = 0, we get 

dpl\^^^ = Lv,{p'')dT + Lvf{p'')dT + ({p\ Vf} - 2Tr(pV/)pO) (dY^ - 2TV(pV/)dr) . 

Since we repeat the control every r time units with r very small, i.e. we take r infinitesimal, 
this leads to the following stochastic time evolution for the density matrix of the two-level 
atom 

dpi = Lv,{p\)dt + Lvs{p\)dt + ({p* , Vf} - 2T,{p\V!)pi) (dYt - 2T,{p\V!)dt) . 

The first term on the right hand side is harmless, we already encountered it in the previous 

section, by taking Kf small enough we can make it as small as we want. The third term 
is also harmless. Since Yt — fl^2TT{p'',Vf)d,r is a martingale it vanishes when we average 
over all possible outcomes for Yt. However, the second term reflects the fact that we can 
not observe Yf and correct it simultaneously with Y^. The next sections are devoted to 
finding a way around this problem. 



4.7 Squeezed states and their calculus 



In this section we drop the assumption that the side channel of the field is initially in the 
vacuum state. We take a step back and rethink our model for (a channel in) the field. For 
the vacuum state we are going to end up with the description we have already used this 
far. Goal of the description below is to incorporate the situation where the initial state of 
the side channel is a so-called squeezed state. In a sqcczcd state the variance of one of the 
quadratures and Yf decreases while the other one increases as a result of Heisenberg's 
uncertainty relation. In the next section we will observe the increased quadrature and 
correct it. The disturbing effect of the other quadrature has decreased as a result of the 
squeezing. 

Let H be the real space of quadratically integrable M^-valued functions on M. On H we 
define a symplectic form a : H x H ^Why 
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where A denotes the Lebesgue measure on M. For notational convenience we define 




We will describe (a channel in) the electromagnetic field by the C*-algebra of canonical 
commutation relations CCR{H,a) over the symplectic space {H,a). 



The algebra CCR{H,a) is defined as the C*-algebra generated by abstract elements 
{W{f); f G H} satisfying relations 

1. wifr = Wi-f), feH, 

2. W{f)W{g) = e^p{-ia{f,g))W{f + g), f,geH. ^ - ' 

The second relation is called the Weyl relation. It follows from [81] that the C*-algebra 
CCR{H, a) exists and moreover that it is unique up to isomorphism. Furthermore it 
immediately follows from (4.23) that W{f) is unitary for all f E H. 



Let a: iifxiJ— >Mbea symmetric positive bilinear form satisfying 

<j{f,9? < aif,f)aig,g), f,geH. (4.24) 

It is well known (cf. [75]) that if a satisfies (4.24) then there exists a unique state 7 on 
the C*-algebra CCR{H, a) satisfying 

7(W^(/))=exp (-!«(/,/)), feH. (4.25) 
Such a state 7 on CCR{H, a) is called a quasifree state. 



In this paper we focus on a particular class of quasifree states 7„c indexed by a parameter 
n G M and a complex parameter c = a + ib,a,b E These states will turn out to be the 
squeezed white noise states of the field as they are encountered in quantum optics after a 
Markov approximation is made (cf. [45]). They are defined through equation (4.25) with 
a symmetric positive bilinear form anc given by [53] 

For notational convenience we define 

_ /2n + 1 + 2a 26 
-"1^ 26 2n + 1 - 2a 

Condition (4.24) leads to the restrictions n(n + 1) > |c|^ and n > 0. For n = c = we get 
the usual vacuum state and for c = we end up with a chaotic temperature state. More 
details on the interpretation of this class of states will follow below. 
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A real linear map J : H ^ H is called multiplication by i if it satisfies = —id. Then 
H is a complex vector space with the usual addition and the scalar multiplication given 
by {x + iy)f = xf + yJf for all x,y eM.. 

Lemma 4.7.1: Let n > and n(n + 1) > |cp. H can be considered as a complex vector 
space equipped with an inner product given by 

if, 9)nc = ancif, g) + ^(T(/, g), f,geH, (4.26) 
if and only if n{n + 1) = |cp. In this case multiplication by i is given by Jnc — JoQnc- 



Proof. Since the inner product (4.26) is linear in its second argument Jnc has to satisfy 
(^{f,Jnc9) = (Xnc{f,g) for all f,g e H. It easily follows from n > and n(n + 1) > |cp 
that Qnc is non degenerate. Therefore Jnc has to satisfy —JoJ,,^. = Qnc which is equivalent 
to Jnc = JoQnc- Jnc is multiplication by i if and only if J^^ = —id, which is equivalent to 

-l\/2n+l + 2a 2b \ /O -l\/2n + l + 2a 2b 

1 o)[ 2b 2n+l-2ayVl o)\ 2b 2n + 1 - 2a 

-1 
-1 

46^- (2n + l + 2a)(2n + l-2a) 

46^- (2n + l + 2a)(2n + l-2a) 



-1 
-1 



I |2 2 , 

c =n + n. 



□ 



In the following we will always be in the situation of Lemma 4.7.1, i.e. n > and n(n+l) = 
|cp. The states of Lemma 4.7.1, i.e. states that allow for the definition of an inner product 
on H through (4.26), arc called Fock states (name will become apparent in a minute). We 
denote the complex Hilbert space given by the pair (H, Jnc) equipped with the inner 
product of (4.26) by Hnc- Note that Hqq is just the space L^(M) of all quadratically 
integrable functions on the real line M. The representation of CCR{H, a) discussed below 
is actually the GNS-representation with respect to a Fock state jnc, see [75] for the details. 

Fix n > and c G C such that n(n + 1) = |cp. Recall that the bosonic Fock space over 
Hnc was defined as 

oo 

-r- — (p /TN fr(S>sk 
k=l 

and that for all / in Hnc the exponential vector is given by e(/) := 1 "^^^ 
span of all exponential vectors was denoted V and the vacuum vector e(0) = 1©0©0© . . . 
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was also written as On the dense domain V we define for all / G -ff„c operators Wndf) 
by ^ 

Wnc{f)e{g) := exp ( - (/,5r)„c - 2"nc(/,/))e(/ + 5), f,geHnc- 

They are isometric and therefore uniquely extend to unitary operators on JP"„c- The map- 
ping n„c : ^(/) Wncif) uniquely defines a linear map II^c from CCR{H,a) into the 
bounded operators on the bosonic Fock space. The map n„c preserves the relations 1. and 
2. of (4.23) defining CCR{H,a), i.e. it is a representation of the canonical commutation 
relations on J^nc- The state 7„c is now given by the vector $ € J-'nc, i-e. 

7„e(^) = n„e(X)$)„„ X e CCR{H, a). 

The triple (J^nc^nc,^) is the GNS-triple corresponding to the state ^nc, cf. [75]. The 
algebra of observables for the electromagnetic field in the Fock state jnc is modelled by 
the von Neumann algebra Wnc generated by {Wndf)', f ^ Hnc}, which is just all bounded 
operators on J^nc- 



Remark. We can reduce the case of a non Fock quasifree state to a Fock state by doubling 
the space H to H (B H. We can embed the algebra of canonical commutation relations 
over H into the algebra of canonical commutation relations over H®H and view the state 
on CCR{H,a) as the restriction of a Fock state on this bigger algebra (cf. [75]). In this 
way we get representations on a doubled up Fock space. Then the algebra of observables 
is not the whole algebra of bounded operators but a true subalgebra. 



The dilation of the semigroup Tt of diagram (4.5) serves as our starting point. We change it 
by replacing the vacuum state (p = 7oo on the side channel by the Fock state jnc described 
above. The dilation diagram then changes to 



B 



IdiS)iig)i 



B 



7nc 



(4.27) 



Coupling the quantum system to a field in another state than the vacuum has changed its 
reduced dynamics to T""^. Changing the representation space of the algebra of canonical 
commutation relations from = J^qq to J^nc also means that we have to describe the joint 
evolution of the system and (the two channels in) the field in this representation. Making 
sense of the group T"*^ will be our main concern for the remainder of this section. 



For all f E H the family of operators {Wnc{tf )}^^^ forms a one-parameter group, conti- 
nuous in the strong operator topology. Therefore it follows from Stone's theorem that for 
all / G if there exists a selfadjoint -Bnc(/) such that 

Wnc{tf)=e^p{itBncif)). 
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The operators i?nc(/) are called field operators. The domain of the operator Bndfk) ■ ■ ■ 
Bncih) contains V for every fi, . . . £ H and k £N (cf. [75]). For f,g e H and t eR 
it follows from the Weyl relation that on the domain P 

1. Bnc{tf) = tBncif), 

2. Bnc{f + g)=Bnc{f)+Bnci9), (4.28) 

3. [Bnc{f),Bnci9)] = 2ia{f,g). 



Let Ho be the real Hilbert space {f & H; f = (/i,0)}. From (4.28.3) it immediately 
follows that the family of operators {Bndf)', f £ Hq} is commutative. Using spectral 
theory, they can be realised as random variables on a single measure space. If the field 
described by the algebra CCR{H, a) is in the Fock state jnc: then the joint characteristic 
function of the random variables Bndfi), Bnc{f2), ■ ■ ■ , Bndfk) is for ti, . . . , G M given 

by 



$,exp [itiBndh)) exp {it2Bndh)) • • - exp [itkB^dlk])^ ^ 

^,exp {iY^tiBndfi))^) = lnc(w (Y^tifi)) =exp( - ^ ^ Utjandfi, fj)) , 

i=l i=l ij=l 

i.e. their joint distribution is Gaussian with covariance matrix andfi^fj)- ^ similar 
way it can be shown that the family {BndJof)i f S Hq} is commutative and the joint 
distribution of the random variables -Bnc (<^o/i)) • • • BndJofk) is Gaussian with covariance 
matrix OindJofii Jofj)- The Gaussianity of these fields, the covariance matrix and the 
condition |cp = + n are exactly the defining properties of a squeezed vacuum state in 
the quantum optics literature, cf. [45]. 

Definition 4.7.2: Fix n G M and c G C such that |c|^ = + On the domain V C Tnc 
we define creation and annihilation operators by 

Keif) ■■= I {Bndf) - iBndJncf)) , ^nc(/) := ^ {Bndf) + iBndJncf)) , f & H, 
Alif) := ^{Bndf) - iBndJof)), ^o(/) := ^ (i?nc(/) + iS„,( Jq/)) , feH. 



It immediately follows from equation (4.28.3) that these operators satisfy the following 
commutation relations [Adf),Ao{g)] = [A*(/), A*(<7)] = [Andf),And9)] = [KdD^KM] 
= 0, [Andf),A*^d9)] = {f^9)nc and [Ao{f),A*{g)] = {f,g)oo for ah f,gGH. Moreover, it 
is a standard result (cf. [75]) that for Fock states Andf)^ = 0, f & H. Furthermore, we 
can build up the symmetric Fock space by acting with creation operators on the vacuum. 
From all these properties it easily follows that for all h, f,g & H 

Andf)49) = {f,9)nce{g), and {e{h), A*^df)e{g)) = /)nc(e(M, 6(5))^ , 
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i.e. Ancif) and satisfy the relations of section 4.3. This means we can define 

stochastic integrals with respect to A^c and A"^^. 

Define the following (non-atomic) projection valued measure ^ on the direct sum Hilbert 
space H = L^(M) © Hnc consisting of a copy of L^(M) for the forward channel and a copy 
of Hnc for the side channel by 



^{I) : L\R) © Hnc ^ L'i^) ®Hnc: 9® (^^^J ^ 9Xl - yj^^^j 

for all Borel subsets / of M. Here xi denotes the indicator function of the set /. Define 
^-martingales by 



/ 
m-' 



The measure ((m'*, m*)) is then given by ((m*, m*)) ([0, t]) = {mf, m'Dnc = (2n + 1 + 2a)t. 
For A„,.(—^===^=) and A*,,.( ,^ '"' ^ ) wc introduce the shorthand notation A^(t) and 
A*g{t), respectively. Note that for stochastic integrals with respect to As(t) and we 
find the Hudson-Parthasarathy Ito table. We denote Ao(mf) and AQ(mf) more compactly 
by Ao{t) and ^o(*)- "^^^ following lemma enables the definition of stochastic integrals 
with respect to AQ{t) and ^o(0- 

Lemma 4.7.3: Let n G M and c G C such that n{n + 1) = |cp. Then for allt>0 we can 
write Ao(t) and Ao(*) ^ following linear combinations of As{t) and A*{t) 

Mt) = //J;' , ^*sit) + /^WV, Mt), 

■\/2n + l + 2a ■\/2n + l + 2a 

A*{t) = , A,{t) + ^±i±^A:(t), 

" ^ V2n + 1 2a ^ ^ V2n + 1 + 2o ' 
where a is the real part of c. 



Proof. Prom Defintion 4.7.2 and J„c = JoQnc it follows that for all / G 

Mf) = \{Bnc{f) + iBnciJof)) = liBncif) + iBnciJncQnl f)) 
= liKcif) + Ancif) - KciQncf) + Anc{Qnc f)) ■ 

Using Jnc = JoQnc and 

Que = 2'^+!+^" ^ 



\ -26 2n -Fl -I- 2a 
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we find for the ^-martingale mt 

n-^m, - ( 2«+i+2a X[o,f] ] _ ^ ( Xlo,t]\ _ j ( Xlo,t] 

'^r^c^'-y -26x[o,t] )~2n + l + 2a\ J 2n + 1 + 2a ' ^ 







( 1 2b X 

V2n + l + 2a 2n + 1 + 2a "V"^*" 



From Definition 4.7.2 and equation (4.28) we see that Anc{Jncf) = —iAndf) and A* j,( J„c/) 
for all f e H. Therefore it follows that 



(2n + i + 2a + 2n +T+ 2a)"^-^^"^*^ ) = 2n + ! + 2a^-^"^*) + 2n + J + 2a^"^^"^*) 
n + c n + l + c 



, --Al{t) + ^(t). 

^271 + 1 + 2a * V2n + 1 + 2a 

□ 



Clearly, we now define for all stochastically integrable processes Lt stochastic integrals 
LtdAoit) and LtdAlU) by Lt , , rf^tft +Lt ,"+\+^, dAsit) and Lf + 

Lf--^^^^==(iA*(t), respectively. Using the Hudson-Parthasarathy Ito table it easily fol- 
lows that the calculus of these stochastic integrals is given by the squeezed noise ltd table 
[45], [53]: 



dMi\dM2 


dAlit) 


dAo{t) 


dAlit) 


cdt 


ndt 


dAo{t) 


{n + l)dt 


cdt 



We are now in a position to explain the construction of Tp'^ in the dilation diagram (4.27). 
The free evolution of the side channel is again given by the unitary group Sf, the second 
quantization of the left shift s{t) on Hnc, i-e. 

In the Heisenberg picture the free evolution on Wnc is then given by X i-^ S*XSt- The 
system 13 and field together form a closed system, thus their joint evolution is given by 
a one-parameter group Ut of unitaries, leading to a Heisenberg picture evolution T^''^ := 
Ad[Ut] on B^Wf® W^^. The group tit is a perturbation of the free evolution. As in the 
vacuum case of section 4.2, we let this perturbation be given by the cocycle of unitaries 



4.8. CONTROL WITH SQUEEZING 



89 



Ut := {S-t S-t)Ut- The stochastic differential equation (4.4) that was satisfied by the 
cocycle Ut when the side channel was still in the vacuum state is now changed. The 
quantum noise of equation (4.8) takes the form 

d(3t = -i{VfdAf{t) - V^dAfit) + VsdAlit) - V;dAo{t)), po = 0, 

in the squeezed noise representation. If the field is in the vacuum state the operators 
AQ^t) and As(t) coincide. AQ^t) should be interpreted as the annihilation operator of a 
photon in the side channel and As{t) should be interpreted as the annihilation operator 
of a squeezed excitation of in the side channel, i.e. a quasiparticle consisting out of many 
photons. Using Lemma 4.7.3 we find 

(..r.?;„T; - .„.".^;.„.F.)dA.(» 



V2n + l + 2a ^ ^/2n + l + 2a ) 
Define 

y - 7^ + 1 + ^ . y _ ^ + g /4 29) 

■ V2n + 1 + 2a ' V2n + 1 + 2a ^ ' ^ ' ' 

then the quantum stochastic differential equation for the cocycle Ut is given by 

dUt = {VfdA}{t) - Vf*dAf{t) + VncdAlit) - V^JAsit) - ^{V^^Vnc + V^Vf)dt}Ut, 
Uo = l. 

(4.30) 

In a similar way as in section 4.3 this leads to the Lindblad operator for the semigroup 
LUX) = Vf*XVf - l{V;Vf, X} + V:,XVnc - ^{KTc^nc, X}, XeB. 



4.8 Control with squeezing 



Note that the operator Vnc of equation (4.29) for strongly squeezed fields, i.e. n and c 
are big, is very close to being skew-selfadjoint. Therefore for strongly squeezed fields the 
dilation is very close to being essentially commutative. In this section we exploit this idea 
and control the skew-selfadjoint part of Vnc- 

Write again V = Vr + iVj with and Vj the selfadjoint operators of equation (4.20). 
We will again use for X G {R,I} and a G {f,s} the notation := k^Vx- Furthermore 
we introduce: 

W, := ^=JS= and W, V^/ + % + c)F.- - i(„ + c)y 
V2n + 1 + 2a V'^n + 1 + 2a 
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i.e. Vnc = Wr + iWi with Wr and Wi selfadjoint. Defining Y^{t) := i{Al{t) - A^{t)) and 
Yf{t) := + Af^{t), a G {f,s} equation (4.30), i.e. the laser is off, becomes 

dUt = {iVfdY/ - iVldvl + iWidV/it) - iWRdY^it) - l{Vf*Vf + V:,Vnc)dt]Ut, 
Uo = l. 

Using a homodyne detection scheme we can observe the quadratures ^^(i) := e~"^Ao(t) + 
e''t'Al{t) for 4> G [0, 27r). With the help of Lemma 4.7.3 this can be written as 

e-'t>{n + c)+e"l>{n+l + c) e'^'f' {n + 1 + c) + e^'t' {n + c) 

V2n + 1 + 2a V2n + 1 + 2a 

For simplicity we assume that c is real, i.e. c = a. Note that the variance of Xq has 
increased due to the squeezing, while the variance of has decreased. Therefore we 
choose to observe Yt := Xo{t) = y/2n + l + 2aYf{t). 

The Belavkin equation for observing Y-t when the laser is still off, follows from equation 
(4.15) (cf. [94]) 

V2n + 1 + 2a V / 

Note that the observed process 1^ is a drift, represented by the term 2Tr(p*V^)(it plus 
an amplified Wiener process, i.e. amplified up to a variance of (2n + 1 + 2a)i. Through 
the drift term we gain information on the state of the two-level system. However, for 
strong squeezing, i.e. n and a big, this information gets lost in the noise of the amplified 
Wiener process. In the limit for squeezing to infinity, the measurement scheme is again 
non-informative, just as in the essentially commutative case. 

We run a control scheme as in section 4.6 only now based on the observation of Yj. We 
correct with the control unitary given by 

= exp — I , 

V V2n + l + 2ay 

where A(t) := 1^ — Yq. Note that for c real, i.e. c = a we have 



W,= '-!hX^ll^±^i^^^ ;j=V2n + l + 2ay/, 

i.e. we can realise this control unitary by applying a laser pulse determined by h{t) = 
'^"^^P Sr{t) for < t < 2r. The control unitary satisfies the following quantum stochastic 
differential equation 

dU- = {-iVfdYr - ^!^L±l±l^vfdT}U^ = Uli-iVfdY, - ^^ + ^ + ^V fdr}, 
U'^ = 1. 



;V2n + l + 2a f 1 
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The state after control is again given by := U^p^U^* where p^ is given by the Belavkin 
equation (4.31). We use the notation below Theorem 4.3.2 with Zi = UJ., = pi and 
Z3 = Ul*. For infinitesimal r evaluated at r = 0, this leads to equation (4.19), i.e. 



dpi 



r=0 



([1] + [2] + [3] + [12] + [13] + [23] + [123]) 



r=0 



Again [123] = and further ([1] + [3] + [13])|^=o = LwAp^)dT + i[p^ , Vf]dYQ. Furthermore 
we have 

[2]|^^^ = LUp')dr + {^[Vf,p^] + ^""^^'^17+2^'' ) ('^° " ^^^p"^^^'-) ' 
([12] + [23]) = -2LwAp'')dT - i [Wi, {Wr, p^}] dr + 2Tv{p''WR)i[Wu p°]dr. 

T=0 

A calculation shows that 

L„e(p°) - LwAp"") {Wr, p°}] = Ly^(p°) + Lw^ip"") - '-[WiWr + WrWi, p\ 

Since WiWr + WrWi = for real c and Ty{p^Wr) [Wi,p^] = Tr(pOyg)[l^/, p°] , we get 



dpi 



T=0 



Lvf{p^)dT+Lwn{p^)dT + 



Since we repeat the control every r time units with r very small, i.e. we take r infinitesimal, 
this leads to the following stochastic time evolution for the density matrix of the two-level 
atom 

The first term on the right hand side is again due to the fact that we did not measure and 
correct the forward channel. It is harmless, since wc can take Kf arbitrarily small. The 
other two terms converge to when squeezing goes to infinity. Therefore, in the limit, 
this control scheme restores quantum information. 
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Summary 



In this thesis the time evolution of open quantum systems is studied. We focus on the 
time evolution of the state of the system conditioned on the result of some measurement 
performed continuously in time in its environment. The goal is to derive a stochastic 
differential equation for the conditional state in which the observed measurement process 
is a stochastic driving term. 

The derivation is carried out within the framework of non- commutative or quantum prob- 
ability theory. In this generalised probability theory there exists a non-commutative gen- 
eralisation of Ito's stochastic calculus. The interaction between an open system and the 
quantized electromagnetic field in the weak coupling limit is governed by a stochastic 
differential equation in this non-commutative sense. The stochastic differential equation 
for the conditional state can be derived from it. It is interpreted as a non-commutative 
generalisation of the Kushner-Stratonovich filtering equation and is called the Belavkin 
filtering equation. 

In Chapter 2 we study a photon counting experiment. The open system studied here, 
is a two-level atom driven by a laser. It emits fluorescence photons in its environment, 
the electromagnetic field. The emitted photons are detected continuously in time. We 
explicitly condition on the observed event by sanwiching with the corresponding projection 
in the field. Using the processes introduced by Davies [33], this leads to a continuous 
time evolution of the reduced system, i.e. the two-level atom, interrupted by jumps at 
the moments at which photons are detected. In quantum optics such trajectories of the 
conditioned state are known as quantum trajectories [27]. 

The third chapter focusses on the derivation of the infinitesimal description of the time 
evolution of the conditional state by the Belavkin filtering equation. Two separate deriva- 
tions of the filtering equations are given. In the first approach we simply differentiate along 
the trajectories for the conditioned state that we obtained in Chapter 2. A diffusive limit 
of the photon counting description enables us to incorporate the situation where instead 
of counting photons in the field we are performing a homodyne detection experiment. 

The second approach uses the decomposition of a von Neumann algebra over its center. 
The observed process in the field determines a commutative algebra C which is the center 
of its commutant A := C . This commutant is the algebra of observables that are not 
demolished by the observation in the field. The central dccompostion enables us to de- 
compose the state restricted to A over the possible paths the observed process yields. In 
this way we are able to condition the state on the observed path. 

The central decompostion is an existence result. It doesn't show how to construct the dc- 
compostion. However, using martingale techniques familiar from classical filtering theory. 
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it is possible to derive a stochastic differential equation for the conditioned state which is 
called the Belavkin filtering equation. Chapter 3 concludes with a recipe for the derivation 
of the Belavkin equation for a wide class of possible continuous time observations in the 
field. 

In the fourth and last chapter we are not only interested in the time evolution of the state 
of an open system conditioned on the result of a measurement in its environment, but we 
also want to use the measurement result for controlling the state evolution. The objective 
is to keep the unknown state of a qubit fixed in time, i.e. to stop the decoherence. 

We show that in a special case where the interaction with the field is such that the 
qubit only couples to commutative noise the above objective can be met. This special 
case is called essentially commutative. However in reality the qubit couples to two non- 
commutative noises. By putting the electromagnetic field in a squeezed state the variance 
of one of the noises increases while the variance of the other decreases. In this way we 
approach the essentially commutative case by stronger and stronger squeezing. We show 
that in the limit for the squeezing strength to infinity the unknown state of the qubit can 
again be kept fixed. 
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Samenvatting 



In dit proefschrift wordt de tijd evolutie van open quantum systemen bestudeerd. We 
concentreren ons op de tijd evolutie van de toestand van het systeem geconditioneerd op 
de uitkomst van een meting continu in tijd in zijn omgeving. Het doel is de afleiding van 
een stochastische differentiaal vergelijking voor de geconditioneerde toestand waarin het 
geobserveerde meetproces als een drijvende stochastische term voor komt. 

De afleiding wordt uitgcvocrd binncn lict kadcr van niet-commutatieve of quantum kans- 
theorie. In deze veralgemeniseerde kanstheorie bestaat een niet-commutatieve genera- 
lisatie van Ito's stochastische calculus. De interactie tussen een open systeem en het 
gequantiseerde electromagnetische veld in de zwakke koppelings limiet wordt beschreven 
door een stochastische differentiaal vergelijking in deze niet-commutatieve zin. De sto- 
chastische differentiaal vergelijking voor de geconditioneerde toestand kan eruit worden 
afgeleid. Deze vergelijking wordt geinterpreteerd als een niet-commutatieve generalisatie 
van de Kushner-Stratonovich filter vergelijking en heet de Belavkin filter vergelijking. 

In Hoofdstuk 2 bcstuderen we een foton detectie experiment. Het bcstudccrde open sys- 
teem is hier een twee niveau systeem gedreven door een laser. Het systeem zendt fluores- 
centie fotonen uit in zijn omgeving, het electromagnetische veld. De uitgezonden fotonen 
worden continu in de tijd gedetecteerd. We conditioneren expliciet op een geobserveerde 
gebeurtenis door de corresponderende projecties om de toestand heen te zetten. Gebruik- 
makend van de processen gei'ntroduceerd door Davies [33] kan de tijd evolutie van het 
twee niveau systeem continu in de tijd afgeleid worden. Dit leidt tot een continue evolutie 
tussen de momenten waarop fotonen worden gedetecteerd en sprongen op de momenten 
waarop fotonen worden gedetecteerd. In quantum optica staan deze trajectorieen van de 
geconditioneerde toestand bekend als quantum trajectorieen [27]. 

Het derdc hoofdstuk richt zich op de infinitesimale beschrijving van de geconditioneerde 
toestand door dc Belavkin filter vergelijking. De filter vergelijkingen worden op twee 
verschillende manieren afgeleid. De eerste methode bestaat uit het eenvoudigweg differen- 
tieren langs een trajectorie van de geconditioneerde toestand. Een diffusieve limiet van de 
foton dctcctic beschrijving maakt het mogelijk om ook homodyne detectie experimenten 
in onze beschrijving op te nemen. 

De tweede methode maakt gebruik van de ontbinding van een von Neumann algebra over 
zijn centrum. Het geobserveerde proces genereert een commutatieve algebra C die het 
centrum is van zijn commutant A := C Deze commutant is dc algebra van observabelen 
die niet gedemoleerd zijn door de observatie in het veld. De centrale ontbinding maakt 
het mogelijk om de toestand beperkt tot A te onbinden over de mogelijke paden van 
het geobserveerde proces. Op deze manier kunnen we de toestand conditioneren op het 
geobserveerde pad. 
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De centralc ontbinding is een existentie rcsultaat. Hct laat nict zien hoc dc ontbind- 
ing geconstrueerd kan worden. Gebruikmakend van martingaal technieken, bekend uit 
klassieke filter theorie, is het mogelijk om een stochastische difFerentiaal vergelijking voor 
de geconditioneerde toestand af te leiden, dc Belavkin filter vergelijking. Hoofdstuk 3 sluit 
af met ccn rcccpt voor dc aflciding van dc Belavkin vergelijking voor een grote klasse van 
mogelijke continue tijd observaties in het veld. 

In het vierde en laatste hoofdstuk zijn we niet alleen gei'nteresseerd in de tijd evolutie van 
de geconditioneerde toestand van een open systeem, maar wc willcn het rcsultaat van de 
continue tijd observatie ook gebruiken om de toestand van het systeem te controleren. 
Het doel is om een onbekende toestand van een qubit vast te houden in de tijd, dat wil 
zeggen de decoherentie tegen gaan. 

We vinden dat dit mogelijk is in het geval dat de interactie met het veld zodanig is dat 
de qubit slechts koppelt met commutatieve ruis. Dit speciale geval heet essentieel com- 
mutatief. In werkelijkheid koppelt de qubit echter met twee onderling niet-commutatieve 
ruizen. Door het electromagnetische veld in een samengedrukte of gesqueezdc toestand te 
nemen, kunnen we de variantie van een van de twee ruizen verkleinen terwijl de variantie 
van de andere wordt vergroot. Op deze manier kunnen we door sterker en sterker squeezen 
het essentieel commutatieve geval beter en beter benaderen. We laten zien dat in de limiet 
voor de squeezing naar oneindig de onbekende toestand van de qubit weer vast kan worden 
gehouden. 



96 



Curriculum Vitae 



Op 11 oktober 1976 word ik, Luc Boutcn, als tweede in het gezin van Piet Bouten en Mia 
Bouten-van Lier geboren te Heldcn. Ik doorlicp de lagere en middelbare school te Pannin- 
gen. Op de middelbare school, het Bouwens van der Boije cohege, werd mijn interesse voor 
wiskunde en natuurkunde gewekt. In 1995 begon ik aan de studie natuurkunde aan de 
universiteit te Nijmegen. Na het eerste studiejaar werd hot zwaartcpunt verlegd naar de 
studie wiskunde die ik vervolgens begin 2000 voltooide. In 2002 volgde de voltooiing van de 
studie natuurkunde. Tot mijn leraren behoorden Hans Maassen, Gert Heckman, Arnoud 
van Rooij, Ted Janssen, Thomas Rijken en Hubert Knops. Tijdens de studie was ik actief 
in een aantal universitaire raden en voorzag ik gedeeltehjk in mijn levensonderhoud door 
middel van het geven van werkcolleges. 

Tegen het einde van mijn studie periode was ik voornamelijk geinteresseerd geraakt in 
quantum mechanica. In het bijzonder de mathematische zijde van de theorie, te weten 
spectraal thcoric, kanstheorie en groepentheorie. Dit leidde in april 2000 tot de keuze 
promotie onderzoek te doen in het gebied van de quantum kanstheorie onder supervisie van 
Hans Maassen. Het promotie onderzoek bestond uit het toepassen van niet-commutatieve 
stochastische analyse op problemen afkomstig uit de quantum optica. Bijzondere aandacht 
ging hierbij uit naar Belavkin's quantum filter theorie. Kennis over quantum optica werd 
opgedaan tijdens een bezoek van drie maanden aan Howard Carmichael in Oregon VS. 
Aansluitend aan de voltooiing van het promotie onderzoek in april 2004 ben ik een aantal 
maanden post-doc geweest bij Slava Belavkin in Nottingham UK. 



97 



98 



Bibliography 



[1] L. Accardi, A. Prigerio, and J. Lewis. Quantum stochastic processes. Publ. RIMS, 
18:97-133, 1982. 

[2] L. Accardi, A. Prigerio, and Y. Lu. The weak couphng hmit as a quantum functional 
central limit. Commun. Math. Phys., 131:537-570, 1990. 

[3] C. Ahn, A. Doherty, and A. Landahl. Continuous quantum error correction via 
quantum feedback control. Phys. Rev. A, 65:042301, 2002. 

[4] C. Ahn, H. Wiseman, and G. Milburn. Quantum error correction for continuously 
detected errors. Phys. Rev. A, 67:052310, 2003. 

[5] G. Alber, C. Charnes, A. Delgado, M. Grassl, and M. Mussinger. Stabilizing dis- 
tinguishable qubits against spontaneous decay by detected-jump correcting quantum 
codes. Phys. Rev. Lett, 86:4402-4405, 2001. 

[6] A. Aspect, J. Dalibard, and G. Rogers. Experimental test of Bell's inequalities using 
time-varying analyzers. Phys. Rev. Lett, 49:1804-1807, 1982. 

[7] A. Barchielli. Measurement theory and stochastic differential equations in quantum 
mechanics. Phys. Rev. A, 34:1642-1649, 1986. 

[8] A. Barchielli. Quantum stochastic differential equations: an application to the elec- 
tron shelving effect. J. Phys. A: Math. Gen., 20:6341-6355, 1987. 

[9] A. Barchielli and V. Belavkin. Measurements continuous in time and a posteriori 
states in quantum mechanics. J. Phys. A: Math. Gen., 24:1495-1514, 1991. 

[10] A. Barchielli and G. Lupieri. Quantum stochastic models of two-level atoms and 
electromagnetic cross sections. J. Math. Phys., 41:7181-7205, 2000. 

[11] V. Belavkin. Quantum filtering of Markov signals with white quantum noise. Ra- 
diotechnika i Electronika, 25:1445-1453, 1980. 

[12] V. Belavkin. Theory of the control of observable quantum systems. Autom. remote 
control, 44:178-188, 1983. 



99 



[13] V. Belavkin. Nondcmolition measurement and control in quantum dynamical systems. 
In A. Blaquiere, S. Diner, and G. Lochak, editors, Proceedings of CISM, Udine, 
Information complexity and control in quantum physics, pages 331-336. Springer- 
Verlag, Wien-New York, 1987. 

[14] V. Belavkin. Nondemolition stochastic calculus in Fock space and nonlinear filtering 
and control in quantum systems. In R. Guelerak and W. Karwowski, editors. Proceed- 
ings XXIV Karpacz winter school, Stochastic methods in mathematics and physics, 
pages 310-324. World Scientific, Singapore, 1988. 

[15] V. Belavkin. A new wave equation for a continuous nondemolition measurement. 
Physics Letters A, 140:355-358, 1989. 

[16] V. Belavkin. Quantum continual measurements and a posteriori collapse on CCR. 
Commun. Math. Phys., 146:611-635, 1992. 

[17] V. Belavkin. Quantum stochastic calculus and quantum nonlinear filtering. Journal 
of Multivariate Analysis, 42:171-201, 1992. 

[18] V. Belavkin. Measurement, filtering and control in quantum open dynamical systems. 
Rep. Math. Phys., 43:405-425, 1999. 

[19] V. Belavkin and O. Melsheimer. A stochastic Hamiltonian approach for quan- 
tum jumps, spontaneous localizations, and continuous trajectories. Quantum Semi- 
class. Opt, 8:167-187, 1996. 

[20] V. Belavkin and P. Staszewski. A quantum particle undergoing continuous observa- 
tion. Physics Letters A, 140:359-362, 1989. 

[21] J. BeU. On the Einstein-Podolsky-Rosen paradox. Physics, 1:195-200, 1965. 

[22] L. Bouten. Squeezing enhanced control, preprint. University of Nijmegen, 2004. 

[23] L. Bouten, M. Gu^a, and H. Maassen. Stochastic Schrodinger equations. J. Phys. A: 
Math. Gen., 37:3189-3209, 2004. 

[24] L. Bouten, H. Maassen, and B. Kiimmerer. Constructing the Davies process of 
resonance fluorescence with quantum stochastic calculus. Optics and Spectroscopy, 
94:911-919, 2003. 

[25] O. Bratteli and D. Robinson. Operator algebras and quantum statistical mechanics L 
Springer- Verlag, Berlin Heidelberg, 1979. 

[26] O. Bratteli and D. Robinson. Operator algebras and quantum statistical mechanics 2. 
Springer- Verlag, Berlin Heidelberg, 1981. 

[27] H. Carmichael. An Open Systems Approach to Quantum Optics. Springer- Verlag, 
Berlin Heidelberg New- York, 1993. 



100 



[28] H. Carmichael. Stochastic Schrodinger equations: what they mean and what they can 
do. In J. Eberly, L. Mandel, and E. Wolf, editors, Coherence and Quantum Optics, 
volume VII. Plenum Press, New York, 1996. 

[29] H. Carmichael, S. Singh, R. Vyas, and P. Rice. Photoelectron waiting times and 
atomic state reduction in resonance fluorescence. Phys. Rev. A, 39:1200-1218, 1989. 

[30] J. Dalibard, Y. Castin, and K. M0lmer. Wave-function approach to dissipative pro- 
cesses in quantum optics. Phys. Rev. Lett., 68:580-583, 1992. 

[31] E. Davies. Quantum stochastic processes. Commun. Math. Phys., 15:277-304, 1969. 

[32] E. Davies. Markovian master equations. Commun. Math. Phys., 39:91-110, 1974. 

[33] E. Davies. Quantum Theory of Open Systems. Academic Press, London New- York 
San Francisco, 1976. 

[34] E. Davies and J. Lewis. An operational approach to quantum probability. Com- 
mun. Math. Phys., 17:239-260, 1970. 

[35] L. Diosi, N. Gisin, J. Halliwell, and I. Percival. Decoherent histories and quantum 
state diffusion. Phys. Rev. Lett, 74:203-207, 1995. 

[36] A. Doherty, S. Habib, K. Jacobs, H. Mabuchi, and S. Tan. Quantum feedback and 
classical control theory. Phys. Rev. A, 62:012105, 2000. 

[37] D. Evans and J. Lewis. Dilations of irreversible evolutions in algebraic quantum 
theory. Communications of the Dublin institute for advanced studies. Series A (The- 
oretical Physics), Dublin institute for advanced studies, 1977. 

[38] R. Feynman. Quantum Electrodynamics. W.A. Benjamin Inc, New York, 1961. 

[39] R. Feynman. Simulating physics with computers. Int. J. Theor. Phys., 21:467-488, 
1982. 

[40] A. Frigerio. Construction of stationary quantum Markov processes through quantum 
stochastic calculus. In L. Accardi and W. von Waldenfels, editors, QP and Appli- 
cations II, volume 1136 of Lecture Notes in Mathematics, pages 207-222. Springer, 
Berlin, 1985. 

[41] A. Frigerio and H. Maassen. Quantum Poisson processes and dilations of dynamical 
semigroups. Prob. Th. Rel. Fields, 83:489-508, 1989. 

[42] C. Gardiner. Inhibition of atomic phase decays by squeezed light: a direct effect of 
squeezing. Phys. Rev. Lett, 56:1917-1920, 1986. 

[43] C. Gardiner and A. Parkins. Driving atoms with light of arbitrary statistics. Phys. 
Rev. A, 50:1792-1808, 1994. 



101 



[44] C. Gardiner, A. Parkins, and P. Zoller. Wave-function quantum stochastic differential 
equations and quantum-jump simulation methods. Phys. Rev. A, 46:4363-4381, 1992. 

[45] C. Gardiner and P. Zoller. Quantum Noise. Springer, Berlin, 2000. 

[46] G. Ghirardi, P. Pcarlc, and A. Rimini. Markov processes in Hilbcrt space and continu- 
ous spontaneous localization of systems of identical particles. Phys. Rev. A, 42:78-89, 
1990. 

[47] N. Gisin. Quantum measurements and stochastic processes. Phys. Rev. Lett., 52:1657- 
1660, 1984. 

[48] N. Gisin. Stochastic quantum dynamics and relativity. Helv. Phys. Acta, 62:363-371, 
1989. 

[49] M. Grcgoratti. The Hamiltonian operator associated to some quantum stochastic 
differential equations. PhD thesis. University Milan, 2000. 

[50] M. Gregoratti. The Hamiltonian operator associated with some quantum stochastic 
evolutions. Commun. Math. Phys., 222:181-200, 2001. 

[51] M. Gregoratti and R. Werner. Quantum lost and found. J. Mod. Opt., 50:915-933, 
2002. 

[52] A. Guichardet. Symmetric Hilhert Spaces and Related Topics, volume 261 of Lect. 
Notes Math. Springer, 1972. 

[53] J. Hellmich, R. Honegger, C. Kostler, B. Kiimmerer, and A. Rieckcrs. Couplings to 
classical and non-classical squeezed white noise as stationary Markov processes. Puhl. 
RIMS, 38:1-31, 2002. 

[54] R. Hudson and K. Parthasarathy. Quantum Ito's formula and stochastic evolutions. 
Commun. Math. Phys., 93:301-323, 1984. 

[55] R. Kadison and J. Ringrose. Fundamentals of the Theory of Operator Algebras, vol- 
ume I. Academic Press, San Diego, 1983. 

[56] R. Kadison and J. Ringrose. Fundamentals of the Theory of Operator Algebras, vol- 
ume II. Academic Press, San Diego, 1986. 

[57] G. Kallianpur. Stochastic Filtering Theory. Springer, Berlin, 1980. 

[58] P. Kochan and H. Carmichael. Photon-statistics dependence of single-atom absorp- 
tion. Phys. Rev. A, 50:1700-1709, 1994. 

[59] B. Kiimmerer. A dilation theory for completely positive operators. PhD thesis, 
Eberhard-Karls-Universitat, Tiibingen, 1982. 

[60] B. Kiimmerer. Markov dilations on VF*-algebras. J. Funct. Anal, 63:139-177, 1985. 



102 



[61] B. Kiimmcrcr. Survey on a theory of non-commutative stationary Markov processes. 
In L. Accardi and W. von Waldenfels, editors, Quantum Probability and Applica- 
tions III, volume 1303 of Lecture Notes in Mathematics, pages 154-182. Springer, 
Heidelberg, 1988. 

[62] B. Kiimmerer. Quantum Markov processes. In A. Buchleitner and K. Hornberger, 
editors. Coherent evolution in noisy environments, volume 611 of Lecture Notes in 
Physics, pages 139-198. Springer- Verlag, Berlin Heidelberg, 2002. 

[63] B. Kiimmerer and H. Maassen. The essentially commutative dilations of dynamical 
semigroups on M„. Commun. Math. Phys., 109:1-22, 1987. 

[64] B. Kiimmerer and H. Maassen. Elements of quantum probability. Quantum Proba- 
bility Communications, X:73-100, 1998. 

[65] H. Kushner. Stochastic stability and control. Academic Press, Inc, New York, 1967. 

[66] H. Kushner. Introduction to stochastic control. Holt, Rinehart and Winston, Inc, 
New York, 1971. 

[67] G. Lindblad. On the generators of quantum dynamical semigroups. Commun. Math. 
Phys., 48:119-130, 1976. 

[68] J. Lindsay and H. Maassen. The stochastic calculus of Bose noise. CWI syllabus 32, 
Centre for Mathematics and Computer science, Amsterdam, 1992. 

[69] H. Maassen. Quantum Markov processes on Fock space described by integral kernels. 
In L. Accardi and W. von Waldenfels, editors, QP and Applications II, volume 1136 
of Lecture Notes in Mathematics, pages 361-374. Springer, Berlin, 1985. 

[70] H. Maassen. Quantum probability applied to the damped harmonic oscillator. In 
S. Attal and J. Lindsay, editors. Quantum Probability Communications, volume XII, 
pages 23-58. World Scientific, Singapore, 2003. 

[71] H. Mabuchi and P. Zoller. Inversion of quantum jumps in quantum optical systems 
under continuous observation. Phys. Rev. Lett., 76:3108-3111, 1996. 

[72] P. Meyer. Quantum Probability for Probabilists. Springer, Berlin, 1993. 

[73] K. M0lmer, Y. Castin, and J. Dalibard. Monte Carlo wave-function method in quan- 
tum optics. J. Opt. Soc. Am. B, 10:524-538, 1993. 

[74] K. Parthasarathy. An Introduction to Quantum Stochastic Calculus. Birkhauser, 
Basel, 1992. 

[75] D. Petz. An Invitation to the Algebra of Canonical Commutation Relations. Leuven 
University Press, Leuven, 1990. 



103 



[76] P. Robinson and H. Maassen. Quantum stochastic calculus and the dynamical Stark 
effect. Rep. Math. Phys, 30:185-203, 1991. 

[77] L. Rogers and D. WiUiams. Diffusions, Markov Processes, and Martingales. Wiley 
Scries in Probability and Mathematical Statistics, 1987. 

[78] I. Segal. Tensor algebras over Hilbert spaces. /. Trans. Am. Math. Sac., 81:106-134, 
1956. 

[79] P. Shor. Polynomial-time algorithms for prime factorization and discrete logarithms 
on a quantum computer. SIAM J. Comp., 26:1484-1509, 1997. 

[80] B. Simon. The P{(j))2 Euclidean quantum field theory. Princeton University Press, 
Princeton, New Jersey, 1974. 

[81] F. Slawny. On factor representations and the C*-algebra of canonical commutation 
relations. Commun. Math. Phys., 24:151-170, 1971. 

[82] M. Srinivas and E. Davies. Photon counting probabilities in quantum optics. Optica 
Acta, 28:981-996, 1981. 

[83] R. Stratonovich. Conditional Markov processes and their application to the theory of 
optimal control. American Elsevier Publishing Company, Inc, New- York, 1968. 

[84] M. Takesaki. Conditional expectations in von Neumann algebras. J. Funct. Anal., 
9:306-321, 1971. 

[85] J. Tomiyama. On the projections of norm one in W^*-algebras. Proc. Japan Acad., 
33:608-612, 1957. 

[86] R. van Handel, J. Stockton, and H. Mabuchi. Feedback control of quantum state 
reduction. arXiv:quant-ph/0402136, Caltech, 2004. 

[87] J. von Neumann. Mathematische Grundlagen der Quantenmechanik. Springer, Berlin, 
1932. 

[88] D. Walls and G. Milburn. Quantum Optics. Springer Verlag, Berlin Heidelberg, 1994. 

[89] D. Williams. Probability with Martingales. Cambridge University Press, Cambridge, 
1991. 

[90] H. Wiseman. Quantum theory of continuous feedback. Phys. Rev. A, 49:2133-2150, 
1994. 

[91] H. Wiseman. Quantum trajectories and feedback. PhD thesis. University of Queens- 
land, 1994. 

[92] H. Wiseman and G. Milburn. Quantum theory of field-quadrature measurements. 
Phys. Rev. A, 47:642-662, 1993. 



104 



[93] H. Wiseman and G. Milburn. Quantum-theory of optical feedback via homodyne 
detection. Phys. Rev. Lett, 70:548-551, 1993. 

[94] H. Wiseman and G. Milburn. All-optical versus electro-optical quantum-limited feed- 
back. Phys. Rev. A, 49:4110-4125, 1994. 



105 



